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STATEMENT  OF  PROBLEM  STUDIED 


Antennas  are  being  used  presently,  or  have  been  proposed  for  use  on  systems  now  under 
development,  on  military  hardware  such  as  missiles,  projectiles,  aircraft,  satellites,  and  tanks. 
They  may  be  used  for  the  purpose  of  remote  guidance,  communications,  terrain  sensing,  fusing, 
navigation,  or  position  location.  In  recent  years  the  frequency  range  of  interest  for  military  systems 
has  gradually  progressed  upward.  For  several  applications  today  frequencies  beyond  those  of  the 
usual  microwave  band  are  required.  Of  special  importance  is  the  millimeter  wave  regions  from 
100  to  300  GHz  (with  a  resulting  wavelength  in  air  of  3.0  to  1 .0  mm). 

Many  of  the  same  properties  that  are  important  for  more  usual  microwave  antennas  are 
also  required  for  uses  at  millimeter  wave  frequencies.  Some  of  the  physical  properties  that  can  be 
important  include,  ruggedness,  ease  of  fabrication,  and  the  ability  to  be  mounted  conformally  to  the 
object  on  which  it  is  attached.  Electrical  properties  which  may  be  important  include  efficiency, 
bandwidth,  input  impedance,  and  radiation  pattern.  In  addition  to  these  usual  considerations,  a 
new  set  of  criteria  also  arises.  Since  the  millimeter  wave  energy  is  not  always  guided  by  standard 
coaxial  lines,  or  even  rectangular  waveguides,  additional  problems  arise  in  the  integration  of 
suitable  antennas  with  the  desired  waveguiding  system.  Thus  the  transition  region  between  the 
waveguide  and  the  antenna  and  the  particular  feeding  structure  is  an  important  part  of  the  antenna 
design  problem. 

The  properties  of  two  types  of  antennas  suitable  for  operation  at  millimeter  wave 
frequencies  (100  to  300  GHz)  have  been  investigated  both  theoretically  and  experimentally.  The 
type  of  radiators  considered  were  printed-circuit  or  microstrip  antennas  and  dielectric  resonator 
antennas.  Special  emphasis  was  placed  on  the  problem  of  integrating  the  antenna  with  the 
waveguiding  structure.  The  ultimate  objective  of  the  research  was  to  provide  aid  in  the  selection 
and  design  of  millimeter  wave  antennas  for  given  systems  applications  within  this  frequency  range. 


SUMMARY  OF  MOST  IMPORTANT  RESEARCH 

1)  Dielectric  cylinders  of  very  high  permittivity  have  been  used  in  the  past  as  resonant 
cavities,  but  since  the  structure  is  not  enclosed  by  metallic  walls,  electromagnetic  fields  do  exist 
beyond  the  geometrical  boundaries  of  the  structure  and  part  of  the  power  is  radiated.  Through  the 
proper  choice  of  geometry  and  permittivity  this  radiation  can  become  the  dominant  feature  of  the 
structure  and  become  an  efficient  antenna  for  use  at  millimeter  wave  frequencies  Both 
experimental  and  theoretical  investigations  of  a  variety  of  these  dielectric  resonator  antennas  have 
been  undertaken.  In  particular,  the  input  impedance  of  a  probe-fed  cylindrical  structure  was 
examined  in  detail  and  a  comparison  of  theoretical  and  experimental  results  was  made.  [4] 

2)  An  approximate  theory  was  derived  for  the  modes  in  a  cylindrical  dielectric  resonator 
antenna.  The  air-dielectric  interface  was  approximated  as  a  perfect  magnetic  conductor  and 
expressions  for  the  resonant  frequency  and  radiation  from  the  equivalent  surface  magnetic  currents 
were  derived.  A  broad  beam  with  a  maximum  normal  to  the  ground  plane  was  obtained  for  the 
lowest  mode.  A  more  directive  beam  or  a  beam  with  a  single  sharp  null  was  obtained  for  the  next 
mode  depending  on  the  cylinder  geometry.  Good  agreement  was  obtained  between  theory  and 
experiment. 

An  experimental  investigation  was  also  undertaken  to  demonstrate  that  the  coplanar 
waveguide  was  an  effective  excitation  mechanism  for  the  cylindrical  dielectric  resonator  antenna. 
Data  showing  the  coupling  behavior  and  the  radiation  patterns  for  various  cylinder  aspect  ratios, 
dielectric  constants,  and  feed  positions  were  presented.  The  results  show  that  efficient  and  practical 
antenna  systems  can  be  designed  using  such  a  feed  structure  which  is  much  more  practical  at 
millimeter  wave  frequencies. 

3)  A  systematic  investigation  of  the  excitation  of  a  cylindrical  dielectric  resonator  antenna 
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by  a  microstrip  transmission  line  was  undertaken  to  illustrate  the  possibility  of  such  a  feed  structure 
and  to  characterize  the  coupling  behavior  and  radiation  patterns  of  the  system. [1 8) 

4)  The  electromagnetic  properties  of  electrically  thick  rectangular  microstrip  antennas  were 
investigated  experimentally.  Antennas  were  fabricated  with  different  patch  sizes,  and  with  electrical 
thicknesses  ranging  from  0.03  to  0.23  wavelengths  in  the  dielectric  substrate.  The  resonant 
frequencies  were  measured  and  compared  to  existing  formulas.  The  bandwidth  was  calculated  as  a 
function  of  electrical  thickness  and  the  antenna  radiation  patterns  were  measured. 

5)  An  algebraic  formulation  for  the  resonant  frequency  of  a  rectangular  microstrip  antenna 
has  been  derived  which  is  valid  for  electrically  thick  substrates.  Predicted  values  of  resonance  are 
compared  with  previous  theories,  as  well  as  actual  experimental  measurements,  for  radiators 
approaching  one-quarter  wavelength  in  thickness,  which  may  correspond  to  physically  thin  printed 
circuit  boards  at  millimeter  wave  frequencies. 

6)  The  design  and  experimental  measurement  of  a  dual-band,  monolithic  microstrip 
antenna  is  presented.  The  structure  utilizes  a  short-circuited  length  of  microstrip  transmission  line 
to  provide  reactive  loading  and.  thereby,  retains  the  low-profile  characteristic  of  a  normal 
microstrip  patch  radiator  and  allows  accurate  fabrication  for  high  frequency  applications. 

7)  An  experimental  study  of  reactively  loaded  rectangular  microstrip  antennas  has  been 
undertaken.  The  reactive  loads  used  were  either  a  single  pin.  or  a  symmetrically  located  pair  of 
pins,  short  circuiting  points  on  the  patch  to  corresponding  points  on  the  ground  plane.  A  set  of 
points  forming  a  single  curve  for  the  case  of  one  shorting  pin  or  two  symmetrically  located  curves 
for  the  case  of  two  shorting  pins  was  found.  A  single,  or  a  pair,  of  short  circuits  placed  at  any 
point  on  these  curves  results  in  an  element  having  a  fixed  resonant  frequency  and  a  fixed  radiation 
pattern  independent  of  the  point  selected.  The  input  impedance,  however,  varies  over  a  very  large 
range  when  the  short  circuit  or  pair  of  short  circuits  is  moved  along  these  curves.  The  major 
application  of  this  work  is  the  study  of  the  feasibility  of  maintaining  an  impedance  match  for  the 
elements  of  a  microstrip  antenna  array  with  changing  scan  angle  by  moving  the  position  of  one  or 
more  short-circuiting  pins.  An  experimental  and  theoretical  study  of  the  mutual  impedance  for 
microstrip  array  elements  was  conducted  and  used  in  a  computer  program  simulating  a  finite  array. 
Though  the  model  used  for  mutual  impedance  was  somewhat  crude,  the  simulation  program 
predicted  a  variation  of  active  element  impedance  with  scan  angle.  For  a  number  of  different  scan 
angles,  the  position  of  shorting  pins  was  adjusted  to  obtain  a  match  for  all  of  the  array  elements 
thus  indicating  that  this  dynamic  matching  technique  may  be  feasible  for  scanned  arrays. 

8)  A  dual  mode  microstrip  antenna  element  was  investigated  which  had  two  independently 
excitable  modes  resonant  at  the  same  frequency.  This  element  has  been  shown  to  be  capable  of 
producing  a  broadside  maximum,  a  broadside  null,  or  an  end-fire  type  pattern  by  suitable  choices 
of  its  reactive  loads  and  suitable  excitation  of  its  degenerate  modes.  Appropriately  located  loads  can 
be  used  to  resonate  modes  normally  resonant  at  quite  different  frequencies,  at  a  single,  common 
frequency.  The  results  indicate  that  the  nodal  lines  of  the  loaded  element  are  accurately  predicted 
by  the  generalized  theory  of  loaded  microstrip  antennas,  and  that  two  modes  can  be  excited 
independently  of  each  other  by  feeding  each  mode  along  the  nodal  line  of  the  other  To  verify  the 
theoretical  predictions  an  actual  dual  mode  microstrip  element  was  fabricated  and  tested.  The 
results  of  this  experiment  correlate  well  with  the  theoretical  model  for  the  overall  characteristics  of 
the  radiator. 

9)  One  of  the  initial  studies  during  this  contract  concerned  the  properties  of  microstrip 
radiators  in  a  two-layered  substrate-superstrate  geometry.  It  was  observed  previously  that  certain 
phenomena  occur  when  a  printed  antenna  element  is  embedded  wi'hin  such  a  structure.  For 
example,  a  superstrate  layer  may  be  used  to  eliminate  surface  waves  excited  by  a  horizontal 
current  element,  resulting  in  100%  efficiency  for  such  an  element.  A  phenomenon  termed 
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"radiation  into  the  horizon"  may  also  he  produced,  in  which  the  far-field  radiation  pattern  of 
an  element  in  a  layered  microstrip  geometry  may  extend  down  to  the  horizon  with  a  nonzero 
value.  Using  this  effect,  a  substrate-superstrale  geometry  may  be  chosen  so  that  the  radiation 
pattern  front  an  embedded  element  will  be  nearly  omnidirectional.  Alternatively,  a  substrate- 
superstrate  geometry  may  be  used  to  produce  the  opposite  effect,  to  produce  patterns  which  have 
narrow  beams.  Some  of  these  effects  may  have  practical  uses  in  modifying  the  properties  of 
microstrip  elements,  as  well  as  the  radiation  patterns.  For  example,  in  [6]  it  is  noted  that  using 
a  superstrate  to  eliminate  the  excitation  of  surface  waves  may  significantly  reduce  the  mutual 
coupling  between  microstrip  elements.  The  effect  of  a  superstrate  layer  on  other  types  of 
microstrip  elements  such  as  the  electromagneticaliy  coupled  (EMC)  dipole  and  the  microstrip 
line  are  also  discussed  in  [6], 

During  the  course  of  this  investigation,  a  technique  for  efficiently  evaluating  the 
Sommerfeld  integrals  which  occur  in  microstrip  analysis  was  used,  which  is  explained  in  [5], 
This  technique  overcomes  most  of  the  difficulty  in  computing  the  Sommerfeld  integrals  for 
source  and  observation  points  which  are  widely  separated,  which  would  otherwise  be  a  serious 
limitation  in  the  computation  of  mutual  impedance  in  many  cases. 

10)  One  of  the  main  projects  which  was  initiated  during  this  contract  was  a  leaky-wave 
investigation  of  the  narrow-beam  effect  which  may  be  created  in  a  multiple-layer  dielectric 
geometry.  This  project,  started  with  an  investigation  of  the  two-layered  substrate-superstrate 
structure,  which  had  already  been  found  previously  to  produce  narrow  beams  as  the 
permittivity  of  the  superstrate  increases,  for  a  particular  choice  of  layer  thicknesses.  This  effect 
was  found  to  be  attributable  to  the  existence  of  weakly  attenuated  leaky  waves  which  are  excited 
on  the  structure  in  this  high-gain  condition.  The  layer  thicknesses  essentially  determine  the 
propagation  constant  of  the  waves,  and  hence  the  beam  scan  angle,  while  the  permittivity 
determines  the  attenuation  of  the  waves,  which  in  turn  determines  the  beamwidths.  An 
asymptotic  solution  of  the  characteristic  equations  for  the  TE  and  TM  leaky  waves  allowed  for 
closed-form  expressions  for  the  propagation  constants.  By  using  scalar  radiation  theory,  radiation 
patterns  for  the  case  of  a  two-dimensional  source  within  the  structure  were  obtained.  A  comparison 
with  exact  patterns  showed  excellent  agreement  near  the  beam  maximum,  where  the  leaky  waves 
contribute  significantly  to  the  total  pattern.  Details  of  the  analysis,  as  well  as  the  results,  are 
given  in  [9]. 

The  next  phase  in  the  leakv-wave  research  project  was  the  investigation  of  the 
multiple-layer  dielectric  structure,  which  consists  of  an  arbitrary  even  number  of  layers  having 
alternating  low  and  high  permittivities.  It  was  observed  in  previous  work  by  others  (H.  Y.  Yang 
and  N.  G.  Alexopoulos.  "Gain  Enhancement  Methods  for  Printed  Circuit  Antennas  Through 
Multiple  Superstates ".  IEEE  Trans.  Antennas  Propagat. .  vol.  AP-35.  pp.  860-863.  July  1987)  that 
very  narrow  beams  could  be  obtained  when  using  more  than  two  layers  in  the  structure.  Such 
a  structure  may  find  applications  at  millimeter-wave  frequencies,  where  the  thickness  would  not 
be  a  serious  limitation.  A  leaky-wave  antenna  constructed  from  only  dielectric  layers,  with 
a  simple  source  mechanism,  would  have  the  advantages  of  durability  and  simplicity.  A 
leaky-wave  analysis  of  this  structure  revealed  the  same  leaky  wave  behavior  as  for  the  two¬ 
layered  case.  A  more  general  analysis  technique  was  developed  to  calculate  the  propagation 
constants  of  the  leaky  waves,  however.  This  technique  utilized  the  similarity  in  field  structure 
between  the  leaky  wave  fields  and  the  fields  existing  inside  the  layered  structure  when  illuminated 
with  a  plane  wave  excitation.  A  power  calculation  inside  the  structure  was  then  used  to  derive 
asymptotic  formulas  for  the  leaky-wave  propagation  constants.  A  comparison  with  exact  values 
of  propagation  constants  showed  very  good  agreement,  which  improves  as  the  beamwidth 
decreases.  The  derivation  of  these  formulas  is  given  in  [15]. 

During  the  leaky-wave  investigation,  another  topic  which  was  explored  was  the 
extension  of  leaky-wave  radiation  formulas  to  the  case  of  cylindrically  propagating  leaky 
waves.  Previously,  leaky-wave  radiation  formulas  existed  only  for  the  case  of  a  one¬ 
dimensional  wave,  where  the  radiation  pattern  could  be  determined  from  a  Kirchoff-Huygens 
scalar  radiation  equation  (T.  Tamir  and  A.  A.  Oliner.  "Guided  Complex  Waves:  Part  2,  Relation 
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to  Radiation  Patterns'.  1'ioc  IF.E.  vol  110.  pp.  325-334.  Fed.  1063).  For  leaky-wave  antennas 
consisting  of  a  cylindrical  leaky  wave  propagating  on  a  planar  interface,  scalar  radiation 
theory  no  longer  suffices,  due  to  the  vector  nature  of  the  aperture  fields.  A  microstrip  source 
within  a  layered  structure  is  one  example  where  the  leaky  wave  is  cylindrical,  so  it  is  clearly 
important  from  a  practical  standpoint  to  have  radiation  formulas  available  for  cylindrical 
waves.  By  using  the  equivalence  principle  and  integrating  over  the  equivalent  sources,  closed- 
form  expressions  for  the  far-field  radiation  patterns  of  TE  and  TM  cylindrical  leaky  waves 
were  obtained.  One  practically  important  aspect  of  these  formulas  was  that  the  results  allowed  for 
arbitrary  aperture  size,  so  that  the  radiation  from  a  truncated  aperture  could  be  studied.  This  is 
practically  important,  since  all  planar  structures  must  be  truncated  at  some  finite  radius.  The 
formulas  which  were  developed  were  used  together  with  the  asymptotic  analysis  of  the  multiple¬ 
layered  structure  in  order  to  accurately  predict  the  radiation  patterns  from  sources  within 
this  structure.  The  radiation  formulas  are  derived  in  [14], 

11)  Another  project  which  was  investigated  under  the  contract  was  the  accurate 
computation  of  mutual  coupling  in  large  phased  arrays  of  patch  antenna  elements.  This 
investigation  led  to  a  formulation  of  mutual  impedance  from  which  series  expansions  for 
the  impedance  matrix  could  be  easily  obtained.  For  large  sparce  arrays,  the  first  few  terms  of 
this  matrix  series  usually  gives  accurate  results  for  the  patch  impedances,  while  requiring 
much  less  computation  time  than  the  matrix  inversion  required  in  a  moment-method 
formulation  of  the  phased  array  impedance  matrix.  A  derivation  of  these  series  expansions  may 
be  found  in  [II], 

12)  Also  initialed  during  the  contract  was  a  design  study  of  an  array  of  EMC  dipoles.  The 
EMC  dipole  is  a  very  promising  microstrip  element  due  to  its  small  size  and  simplicity.  At 
higher  frequencies  where  the  substrate  thickness  may  be  .05  wavelengths  or  more,  the  EMC 
dipole  also  exhibits  a  large  bandwidth,  ft  also  has  the  advantage  of  easy  matching  with  the 
feeding  transmission  line,  by  simply  adjusting  the  dipole  length  and  offset.  An  array  of  these 
elements  would  make  an  attractive  design  for  many  applications.  Previously,  a  design  procedure 
for  an  array  of  such  elements  was  developed  by  Elliott  and  Stern  (R.  S.  Elliott  and  G.  J.  Stern. 
"The  Design  of  Microstrip  Dipole  Arrays  including  mutual  Coupling.  Part  I:  Theory’’.  IEEE 
Trans.  Antennas  Propagal..  vol.  AP-29.  pp.  757-760.  Sept.  1981).  This  design  procedure  was 
more  directly  compatible  with  experimental  measurements  than  with  numerical  calculation, 
however.  During  the  present  investigation  a  design  procedure  was  developed  which  was  well- 
suited  to  a  direct  moment-method  analysis.  One  of  the  key  features  of  this  design  procedure  is 
the  introduction  of  coupling  coefficients  which  may  be  determined  numerically  by 
considering  only  two  EMC  dipole  pairs  at  a  time,  resulting  in  a  relatively  small  moment 
method  problem.  The  use  of  traveling-wave  basis  functions  on  (he  transmission  lines  further 
improved  the  efficiency  of  the  calculation.  Once  the  coupling  coefficients  are  determined,  a 
matrix  description  of  the  line  and  dipole  currents  is  used  to  formulate  a  matrix  set  of  design 
equations,  which  accounts  for  all  of  the  necessary  interactions  between  all  of  the  lines  and  dipoles 
in  the  array.  This  matrix  set  of  equations  allow-s  for  the  design  of  a  feed  network  which  will  excite 
the  required  dipole  currents  in  order  to  produce  a  specified  radiation  pattern.  In  order  to  test 
the  theory,  an  experimental  four  element  linear  array  was  constructed  and  tested,  with  good 
agreement  between  the  theoretical  and  measured  patterns.  Details  of  the  design  procedure, 
along  with  the  results,  may  be  found  in  [13], 

13)  An  investigation  of  the  scattering  properties  of  conductive  bodies  embedded  within 
lossy  media  was  undertaken  during  the  contract  period.  One  particular  example  of  such  a 
configuration  would  be  a  microstrip  patch  antenna  on  top  of  or  below  a  lossy  dielectric  layer.  A 
general  scattering  formula  was  derived  which  shows  the  fundamental  dependence  of  the  radar 
cross  section  of  the  body  on  the  parameters  of  the  lossy  meo.a.  One  important  result  of  this 
formula  is  that  the  body  radar  cross  section  (RCS)  is  directly  related  to  the  radiation  efficiency  of 
the  body,  so  that  the  RCS  may  be  directly  controlled  by  introducing  loss  into  a  surrounding  layer. 
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A  numerical  study  was  performed  for  a  rectangular  patch  antenna  embedded  within  a  lossy 
substrate-superstrate  geometry,  as  a  specific  case.  Results  were  obtained  for  three  particular 
configurations:  a  lossy  substrate,  a  lossless  substrate  with  a  lossless  superslrate  on  top.  and  a 
lossless  substrate  with  a  lossy  superstrate  on  top.  In  each  case,  the  results  confirmed  the  general 
trend  predicted  by  the  scattering  formula.  Details  and  results  may  be  found  in  [12]. 

14)  During  the  contract  period  a  feasibility  study  of  an  array  of 
electromagnetically  coupled  patch  antennas  was  performed.  In  this  array,  the  patch  antennas 
are  electromagnetically  coupled  between  adjacent  neighbors  by  short  sections  of  transmission 
line  which  are  embedded  within  the  substrate.  This  eliminates  the  need  for  a  corporate  feed 
network,  and  greatly  simplifies  the  array  construction.  The  feasibility  study  determined  that  it  is 
possible  to  design  such  an  array  to  produce  a  beam  scanned  to  any  specified  direction,  with 
reasonable  coupling  values  between  adjacent  patches  required.  A  moment  method  calculation 
was  used  to  determine  the  coupling  coefficients  between  adjacent  patches,  and  a  matrix  equation 
using  these  values  was  then  formulated  to  predict  the  patch  current  amplitudes  assuming  a 
given  feed  source  on  the  first  patch  and  a  specified  load  termination  on  the  last  patch.  A 
synthesis  procedure  was  then  developed,  whereby  the  coupling  coefficients  between  adjacent 
patches  may  be  determined  from  the  desired  set  of  patch  current  amplitudes,  which  determines 
the  radiation  pattern.  Hence,  the  dimensions  of  the  embedded  line  sections  between  the  patches  are 
determined  from  the  specified  pattern.  To  test  the  theory,  a  four  element  linear  array  of  patches  was 
built,  designed  for  a  45  beam.  The  measured  pattern  showed  good  agreement  with  the 
predicted  result,  and  demonstrated  the  feasibility  of  this  type  of  array.  The  details  of  the  design 
technique,  along  with  the  results,  are  given  in  (1 1). 

15)  An  investigation  of  the  properties  of  a  traveling- wave  array  of  monopoles 
embedded  within  a  substrate,  fed  by  a  meandering  stripline,  was  conducted  during  the  contract 
period.  Such  an  array  may  make  a  useful  scanned-beam  antenna  for  some  applications. 
Because  individual  monopoles  which  are  electrically  short  do  not  radiate  significantly,  it  is 
possible  to  construct  arrays  of  this  type  many  wavelengths  long,  with  control  of  the  monopole 
currents  achieved  by  variation  of  the  monopofe  heights  and  length  of  the  meandering  stripline 
between  monopoles.  A  quasi-periodic  structure  analysis  of  the  array  was  used  to  obtain  design 
equations  for  the  monopolc  heights  and  line  lengths  required  to  produce  a  specified  constant 
phase  shift  between  elements,  and  a  specified  amplitude  taper  in  the  monopole  currents  along 
the  length  of  the  array.  As  an  experimental  verification  of  the  theory,  a  ten-element  linear  array 
with  a  60°  beam  was  designed  and  constructed.  The  measured  pattern  showed  good 
agreement  with  the  theoretical  one.  A  detailed  discussion  of  the  theory  and  results  may  be 
found  in  [2], 
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[I]  S.  E.  Davidson.  S.  A.  Long  and  W.  F.  Richards.  "Dual-Band  Microstrip  Antennas  with 
Monolithic  Reactive  Loading".  Electronics  Letters,  vol.  21.  no.  20.  pp.  936-937.  September  26. 
1985. 


The  design  and  experimental  measurement  of  a  dual-band,  monolithic  microstrip  antenna  is 
presented.  The  structure  utilizes  a  short-circuited  length  of  microstrip  transmission  line  to  provide 
reactive  loading  and,  thereby,  retains  the  low-profile  characteristic  of  a  normal  microstrip  patch 
radiator. 
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[2]  F..  Chang.  S.  A.  Long,  and  W.  F.  Richards.  "An  Experimental  Investigation  of 
Electrically  Thick  Rectangular  Microstrip  Antennas”.  IEEE  Trans.  Aniennas  and  Propagation,  vol. 
AP-34,  June  1986. 

The  electromagnetic  properties  of  electrically  thick  rectangular  microstrip  antennas  were 
investigated  experimentally.  Antennas  were  fabricated  with  different  patch  sizes  and  with  electrical 
thicknesses  ranging  from  0.03  to  0.23  wavelengths  in  the  dielectric  substrate.  The  resonant 
frequencies  were  measured  and  compared  to  existing  formulas.  The  bandwidth  was  calculated  as  a 
function  of  electrical  thickness  and  the  antenna  radiation  patterns  were  measured. 

[3]  W  F.  Richards  and  S.  A.  Long.  “Reactively  Loaded  Microstrip  Antennas."  IEEE 
Antennas  and  Propagation  Society  Newsletter,  vol.  28.  No.  5.  pp.  11-17.  October  1986. 

Microstrip  elements  have  proven  to  be  very  versatile  and  have  found  uses  in  a  variety  of 
applications  from  satellite  arrays  to  automobile  anti-collision  warning  systems.  This  versatility 
stems  mainly  from  their  structural  simplicity,  ease  of  construction,  ruggedness,  conformability,  low 
profile,  and  low  cost.  In  addition  to  these  features,  microstrip  elements  are  capable  of  producing  a 
variety  of  different  patterns  and  polarization.  Most  research  efforts,  however,  have  been  limited  to 
basic  applications  using  only  a  single,  unperturbed  mode  of  a  rectangular  or  circular  element. 
Additional  applications  are  possible  through  the  used  of  reactive  loading.  These  include  multiband, 
frequency-agile  and  polarization-diverse,  adaptable-impedance,  and  multimode-pattern  antennas. 
Some  experimental  work  has  been  done  in  these  areas  but  very  little  theoretical  work  has  been 
done.  In  this  paper,  a  very  simple  theory  is  presented  which  will  guide  a  designer  who  wished  to 
take  advantage  of  the  added  flexibility  that  reactive  loading  provides.  In  addition,  several  application 
are  discussed  with  corresponding  theoretical  and  experimental  results. 

[4]  S.  A.  Long  and  M.  W.  McAllister.  "The  Input  Impedance  of  the  Dielectric  Resonator 
Antenna".  International  Journal  of  Infrared  and  Millimeter  H'ovcs.  vol.  7.  no. 4.  pp.  555-570. 
1986 


Dielectric  cylinders  of  very  high  permittivity  have  been  used  in  the  past  as  resonant 
cavities,  but  since  the  structure  is  not  enclosed  by  metallic  walls,  electromagnetic  fields  do  exist 
beyond  the  geometrical  boundaries  of  the  structure  and  part  of  the  power  is  radiated.  Through  the 
proper  choice  of  geometry  and  permittivity  this  radiation  can  become  the  dominant  feature  of  the 
structure  and  become  an  efficient  antenna  for  use  at  millimeter  wave  frequencies.  Both 
experimental  and  theoretical  investigations  of  a  variety  of  these  dielectric  resonator  antennas  have 
been  undertaken.  In  particular,  the  input  impedance  of  a  probe-fed  cylindrical  structure  was 
examined  in  detail  and  a  comparison  of  theoretical  and  experimental  results  was  made. 

[5]  D.  R.  Jackson  and  N.  G.  Alexopoulos.  "An  Asymptotic  Extraction  Technique 
for  Evaluating  Sommerfeld-Type  Integrals’’.  IEEE  Transactions  on  Antennas  and  Propagation. 
vol.  AP-34.  pp.  1 467-1470.  Dec.  1986. 

An  asymptotic  extraction  technique  is  developed,  which  leads  to  an  efficient  evaluation  of 
the  electric  field  from  a  Hertzian  dipole  in  a  layi  -ed  geometry.  This  technique  allows  for  a 
rapidly  converging  expression  for  the  electric  field,  which  remains  well  behaved  as  the  source 
and  observation  heights  above  the  ground  plane  coincide.  A  useful  application  of  this  method  is 
in  the  calculation  of  mutual  impedance  between  printed  dipoles. 


[6]  D.  R.  Jackson  and  N.  G.  Alexopoulos.  "Analysis  of  Planar  Strip  Geometries  in 
a  Substrate-Superstrate  Configuration”.  IEEE  Transactions  on  Antennas  and  Propagation,  vol.  AP- 
34.  pp.  1430-1438.  Dec.  1986. 
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A  moment  method  procedure  is  used  to  analyze  the  oehavior  of  several  different 
configurations  consisting  of  planar  strips  in  a  substrate-superstrate  geometry  These  include 
the  microstrip  transmission  line. 

center- fed  dipole,  the  mutual  impedance  between  two  dipoles,  and  the  transmission-line 
coupled  dipole.  In  each  case  some  of  the  basic  superstrate  effects  are  discussed. 


[7]  R.  Garg  and  S.  A.  Long.  "Resonant  Frequency  of  Electrically  Thick  Rectangular 
Microstrip  Antennas"  Electronics  Letters,  vol.  23.  no.  21.  pp.  1 149-1 151.  Octobei  8.  1987. 

An  algebraic  formulation  for  the  resonant  frequency  of  a  rectangular  microstrip  antennas 
has  been  derived  which  is  valid  for  electrically  thick  substrates.  Predicted  values  of  resonance  are 
compared  with  previous  theories,  as  well  as  actual  experimental  measurements,  for  radiators 
approaching  one-quarter  wavelength  in  thickness. 


[8]  R.  Garg  and  S.  A.  Long.  "An  Improved  Formula  for  the  Resonant  Frequencies  of  the 
Triangular  Microstrip  Patch  Antenna".  IEEE  Trans,  on  Antennas  and  Propagation.  vol.AP-  36. 
no.  4.  pp.  570-575,  April  1988. 

An  improved  formula  is  derived  to  predict  the  resonant  frequency  of  the  triangular 
microstrip  patch  antenna.  The  tecltni  ue  previously  used  for  the  case  of  a  circular  disk  is  modified 
to  produce  greater  accuracy  than  pre\  ious  investigations. 


(9]  D.  R.  Jackson  and  A.  A.  Oliner.  'A  Leaky-Wave  Analysis  of  the  High-Gain  Printed 
Antenna  Configuration".  IEEE  Transactions  on  Antennas  and  Propagation,  vol.  AP-36.  July 
1988. 

A  leaky-wave  analysis  is  used  to  explain  the  narrow-beam  "resonance  gain"  phenomenon, 
in  which  narrow  beams  may  be  produced  from  a  printed  antenna  element  in  a  substrate- 
superstrate  geometry.  The  leaky-wave  approach  furnishes  insight  into  the  physical  nature  of  this 
phenomenon,  as  well  as  providing  new  practical  information.  The  leakv-wave  approach  also 
permits  the  derivation  of  simple  asymptotic  formulas  for  the  leaky-wave  propagation  constants, 
which  determine  the  radiation  patterns. 

[10]  R.  A.  Kranenburg  and  S.  A.  Long  "Microstrip  Transmission  Line  Excitation  of 
Dielectric  Resonator  Antennas".  Electronics  Letters,  vol.  24.  no.  18.  pp.  1156-1157.  September  I. 
1988. 


A  systematic  investigation  of  the  excitation  of  a  cylindrical  dielectric  resonator  antenna  was 
undertaken  to  illustrate  the  possibility  of  such  a  feed  structure  and  to  characterize  the  coupling 
behavior  and  radiation  patterns  of  the  system. 


PAPERS  SUBMITTED 

[11]  D.  R.  Jackson.  W.  F.  Richards,  and  A.  Ali-Khan.  "Series  Expansions  for  the  Mutual 
Coupling  in  Microstrip  Patch  Arrays".  IEEE  Transactions  an  Antennas  and  Propagation  (accepted 
for  publication). 

A  dominant-mode  mutual  coupling  theory  is  developed  for  an  array  of  microstrip  patches.  A 
formulation  is  used  which  is  compatible  with  either  a  cavity  model  or  a  spectral-domain 
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calculation,  or  with  measurements.  Th's  formulation  allows  for  the  simple  derivation  of  series 
expansions  for  the  impedance,  admittance,  and  scattering  matrices,  which  do  not  require 
inversion  of  a  moment-method  matrix.  These  expansions  are  computationally  efficient  for  large, 
sparce  arrays. 


[12]  D.  R.  Jackson.  "The  RCS  of  a  Rectangular  Microstrip  Patch  in  a  Substrate- 
Superstrate  Geometry".  IEEE  Transactions  on  Antennas  and  Propagation  (accepted  for 
publication). 

A  general  scattering  formula  is  developed  for  the  radar  cross  section  of  a  general  resonant 
conductive  body  embedded  within  lossy  media.  This  formula  shows  that  the  RCS  of  the  body  is 
directly  related  to  the  radiation  efficiency  of  the  body.  The  specific  case  of  a  rectangular  microstrip 
patch  within  a  lossy  substrate-superstrate  geometry  is  analyzed,  with  results  presented  for  three 
different  cases:  a  lossy  substrate  only,  a  lossless  substrate  with  a  lossless  superstrate.  and  a 
lossless  substrate  with  a  lossy  superstrate.  It  is  concluded  that  the  lossy  superstrate  case  gives 
the  most  RCS  reduction  with  the  least  detrimental  effect  on  the  input  impedance  of  the  patch. 

[13]  D.  R.  Jackson.  A.  E.  Dinbergs.  and  S.  A.  Long.  "A  Moment  Method  Design  Procedure 
for  an  Array  of  EMC  Dipoles",  submitted  to  IEEE  Transactions  on  Antennas  and  Propagation. 

A  design  procedure  based  on  a  moment-method  analysis  is  given  for  an  array  of  microstrip 
dipoles  electromagnetically  coupled  (EMC)  to  microstrip  transmission  lines.  One  of  the  key  features 
of  the  design  procedure  is  that  only  one  pair  of  EMC  dipoles  needs  to  be  analyzed  at  a  time 
numerically,  in  order  to  develop  the  necessary  coupling  coefficients  which  account  for  the 
interactions  between  all  of  the  EMC  line-dipole  pairs.  A  convenient  set  of  matrix  design  equations 
are  obtained,  which  allows  for  the  design  of  a  feed  network  to  excite  the  desired  dipole  currents. 


[14]  A.  Ip  and  D.  R.  Jackson.  "Radiation  from  Cylindrical  Leaky  Waves" .  submitted  to 
IEEE  Transactions  on  Antennas  and  Propagation. 

Formulas  for  the  far-field  pattern  of  cylindrical  leaky  waves  propagating  on  a  planar 
surface  are  derived.  These  formulas  may  be  used  to  predict  the  radiation  pattern  of  a  general 
class  of  leaky-wave  antennas,  consisting  of  a  finite-size  source  which  excites  a  radially 
propagating  leaky  wave  on  some  planar  surface.  Leaky-wave  antennas  comprised  of 
microstrip  antenna  elements  fall  into  this  category.  Formulas  are  derived  for  both  TE  and  TM 
leaky  waves  of  arbitrary  aperture  size,  which  allows  the  practical  effects  of  a  truncated  aperture  to 
be  determined. 


[15]  D.  R.  Jackson.  A.  A.  Oliner.  and  A.  Ip.  "Leaky-Wave  Propagation  and  Radiation  for 
a  Multiple-Layer  Dielectric  Structure",  submitted  to  IEEE  Transactions  on  Antennas  and 
Propagation. 

A  leaky-wave  analysis  is  used  to  predict  the  beamwidth  of  a  narrow-beam  multiple  layer 
dielectric  structure,  consisting  of  an  alternating  arrangement  of  low>  and  high  permittivity 
dielectric  layers.  A  power-flow  analysis  is  used  to  derive  asymptotic  formulas  for  the  phase 
and  attenuation  constants  of  the  leaky  waves  which  propagate  on  this  structure,  which  directly 
determines  the  beamwidths  and  scan  angle.  The  asymptotic  formulas  show  that  the  beam  width 
decreases  very  rapidly  as  the  number  of  layers  increases,  so  that  very  narrow  beams  may  be 
obtained  with  only  a  modest  number  of  layeis. 

[16]  H.  Y.  Yang.  N.  G.  Alexopoulos.  and  D.  R.  Jackson.  ''Microstrip  Open-End  and  Gap 
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Discontinuities  in  a  Substrate-Superstrate  Structure",  submitted  to  IEEE  Transactions  on 
Microwave  Theory  and  Techniques. 

An  analysis  of  both  open-end  and  gap  discontinuities  in  niicrostrip  lines  within  a  substrate- 
superstrate  geometry  is  given.  The  analysis  uses  both  subdomain  and  semi-infinite  traveling-wave 
basis  functions  for  an  efficient  moment-method  calculation  of  the  scattering  parameters. 
Results  show  that  the  superstrafe  may  significantly  affect  the  discontinuity  parameters,  by 
increasing  the  space  wave  and  surface-wave  radiation  at  the  discontinuity. 


]I7]  IV  13.  Katchi.  1).  R.  Jackson,  and  N.  G.  Alcxopoulos.  "Microslrip  Dipoles", 
submitted  as  Chapter  5  of  The  Handbook  a/  Microslrip  Aniennas.  .1.  R.  James  and  P.  S.  Hall. 
Eds. 

A  general  overview  of  (he  analytic  and  numerical  aspects  of  modeling  microstrip  dipole 
configurations  is  presented,  along  with  results  showing  some  of  the  basic  properties  of  the 
niicrostrip  dipole.  Particular  attention  is  devoted  to  the  EMC  dipole,  both  as  a  single  element  and  as 
an  element  in  a  finite  array.  Practical  considerations  include  bandwidth  and  impedance  behavior, 
and  the  trade-off  between  dipole  bandwidth  and  line  radiation  for  the  EMC  dipole. 

[18]  R.  A.  Kranenburg.  S.  A.  Long  and  .1.  T.  Williams.  "  Coplanar  Waveguide  Excitation 
of  Dielectric  Resonator  Antennas".  IEEE  Transactions  on  Antennas  and  Propagation. 

An  approximate  theory  is  derived  for  the  modes  in  a  cylindrical  dielectric  resonator 
antenna.  The  air-dielectric  interface  is  approximated  as  a  perfect  magnetic  conductor  and 
expressions  for  the  resonant  frequency  and  radiation  from  the  equivalent  surface  magnetic  currents 
are  presented.  A  broad  beam  with  a  maximum  normal  to  the  ground  plane  is  obtained  for  the 
lowest  mode.  A  more  directive  beam  or  a  beam  with  a  single  sharp  null  is  obtained  for  the  next 
mode  depending  on  the  cylinder  geometry.  Good  agreement  is  obtained  between  theory  and 
experiment.  An  experimental  investigation  was  also  undertaken  to  demonstrate  the  coplanar 
waveguide  as  an  effective  excitation  mechanism  for  the  cylindrical  dielectric  resonator  antenna. 
Data  showing  the  coupling  behavior  and  the  radiation  patterns  for  various  cylinder  aspect  ratios, 
dielectric  constants,  and  feed  positions  are  presented.  The  results  show  that  efficient  and  practical 
antenna  systems  can  be  designed. 


PAPERS  PRESENTED 

[1]  S.  A.  Long,  Dielectric  Resonator  Antennas"  presented  at  U.  S.  Army  Research  Office 
Workshop  on  Millimeter  Wave  Technology.  New  York  Institute  of  Technology.  Old  Westbury. 
NY.  December  1984. 

[2]  S.  E.  Davidson  S.  A.  Long,  W.  F.  Richards.  "Monolithic  Design  of  Dual-Band  Microstrip 
Antennas  Using  Reactive  Loading”  presented  at  IEEE  International  AP-S  Symposium,  Vancouver. 
British  Columbia.  June  1985 

[3]  W.  F  Richards  and  S.  A.  Long,  "Pattern  Adaptation  Using  Loaded  Dual-Mode  Microstrip 
Antennas  presented  at  IEEE  International  AP-S  Svmposium.  Vancouver.  British  Columbia.  June 

1985 

[4]  M  W.  McAllister  anil  S.  A.  Long.  "The  Dielectric  Resonator  Antenna",  presented  at  Tenth 
International  Conference  on  Infrared  and  Millimeter  Waves.  Lake  3uena  Vista.  Florida.  December 
1985 
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[51  W.  F,  Richards  and  S.  A.  Long.  "Impedance  Matching  of  Microstrip  Antennas  with  Reactive 
Loads"  presented  at  National  Radio  Science  Meeting,  Boulder.  Colorado.  January  1986. 
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presented  at  IEEE  International  AP-S  Symposium.  Philadelphia.  PA.  June  1986 
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Meeting.  Philadelphia.  PA.  June  1986. 

[10]  S.  A.  Long.  "The  Cylindrical  Dielectric  Resonator  Antenna"  presented  at  the  U.  S.  Army 
Research  Office  Workshop  on  Fundamental  Issues  in  Millimeter  and  Suhmillimeter  Waves.  Los 
Angeles,  CA,  September  1986. 

[11]  W.  F.  Richards  and  S.  A.  Long.  "Impedance  Control  of  Microstrip  Antennas  Utilizing 
Reactive  Loading"  presented  at  International  Telemetering  Conference.  Las  Vegas.  NV.  October 
1986. 

[12]  W.  F.  Richards  and  S.  A.  Long.  "Adaptive  Pattern  Control  of  a  Reactively  Loaded  Dual-Mode 
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DEGREES  EARNED 

(I)  E.  Chang,  "An  Experimental  Study  of  Electrically  Thick  Rectangular  Microstrip 
Antennas."  M.S.E.E..  University  of  Houston.  December  1985. 

The  electromagnetic  properties  of  electrically  thick  rectangular  microstrip  antennas  are 
investigated.  Antennas  are  fabricated  with  different  patch  sizes  and  with  electrical  thicknesses 
ranging  from  0.03  to  0.25  wavelengths  in  the  dielectric  substrate.  The  resonant  frequency  and 
bandwidth  are  measured  and  compared  to  existing  theories.  The  radiation  patterns  are  also 
measured  and  the  3  dB  beamwidths  for  the  antennas  are  determined. 


(2)  S.  Toomsawasdi.  "Moment  Method  Analysis  of  a  Microstrip  Antenna  Using  a  Patch-Edge 
Integral  Equation".  Ph  D.  University  of  Houston.  May  1987. 

A  method  of  computing  the  input  impedance  of  a  microstrip  antenna  in  free  space  and  in  a 
modified  environment  is  presented.  A  pair  of  integral  equations  in  terms  of  the  patch-edge  electric 
and  magnetic  currents  are  solved  using  the  method  of  moments.  One  of  these  integral  equations  is 
based  on  an  approximate  exterior-region  Green's  function.  An  appropriately  chosen  magnetic 
current  distributed  on  the  side  of  a  conducting  mesa  in  the  presence  of  the  dielectric  is  capable  of 
reproducing  the  exact  field  of  a  patch  element.  This  mesa  is  formed  by  connecting  the  edge  of  the 
patch  to  the  ground  plane  with  conductor.  The  approximate  Green's  function  used  in  this  theory 
ignores  the  presence  of  both  the  dielectric  substrate  and  the  conducting  mesa.  It  is.  therefore,  just 
the  free-space  Green's  function  and  its  use  instead  of  an  exact  Green's  function  considerably 
reduces  the  computational  complexity  of  the  analysis.  Results  for  the  input  impedance  computed 
using  other  approximate  theories  and  measurements  are  presented  for  comparison  with  this  theory. 
Comparison  with  measured  results  shows  the  this  theory,  though  approximate,  very  accurately 
predicts  the  locus  of  the  input  impedance  of  microstrip  ar'ennas  of  practical  thicknesses. 
Surprisingly,  the  antenna  quality  factor  was  also  predicted  very  accurately  by  this  technique 
although  the  approximate  Green's  function  cannot  account  for  the  presence  of  a  surface  wave.  The 


11 


resonant  frequencies  predicted  by  this  theory  are  somewhat  higher  than  those  measured,  an  effect 
attributable  to  ignoring  the  presence  of  the  dielectric.  The  approximate  Green's  function  was  also 
modified  to  account  for  the  presence  of  an  elevated  dielectric  slab  for  applications  to  microstrip 
image  arrays.  Measurements  were  made  for  an  element  in  this  modified  environment; 
theoretically  predicted  impedances  agreed  well  with  the  measured  results. 


(.3)  A.  Ali  Khan.  "Study  of  Loaded  Microstrip  Antennas  and  Their  Application  To  Arrays." 
M.S.E.E.,  University  of  Houston.  December  1987. 

An  experimental  study  of  reactively  loaded  rectangular  microstrip  antennas  is  presented  in 
this  thesis.  The  reactive  loads  used  here  were  either  a  single  pin.  or  a  symmetrically  located  pair 
of  pins,  short  circuiting  points  on  the  patch  to  corresponding  points  on  the  ground  plane.  A  set  of 
points  forming  a  single  curve  for  the  case  of  one  shorting  pin  or  two  symmetrically  located  curves 
for  the  case  of  two  shorting  pins  was  found.  A  single  or  a  pair,  of  short  circuits  placed  at  any 
point  on  these  curves  results  in  an  element  having  a  fixed  resonant  frequency  and  a  fixed  radiation 
pattern  independent  of  the  point  selected.  The  input  impedance,  however,  varies  over  a  very  large 
range  when  a  short  circuit  or  pair  of  short  circuits  is  moved  along  these  curves.  The  major 
application  of  this  work  is  the  study  of  the  feasibility  of  maintaining  an  impedance  match  for  the 
elements  of  a  microstrip  antenna  array  with  changing  scan  angle  by  moving  the  position  of  one  or 
two  short-circuiting  pins.  An  experimental  and  theoretical  study  of  the  mutual  impedance  for 
microstrip  array  elements  was  conducted  and  used  in  a  computer  program  simulating  a  finite  array. 
Though  the  model  used  for  mutual  impedance  was  somewhat  crude,  the  simulation  program 
predicted  a  variation  of  active  element  impedance  with  scan  angle.  For  a  number  of  different  scan 
angles,  the  position  of  shorting  pins  was  adjusted  to  obtain  a  match  for  all  of  the  array  elements 
thus  indicating  that  this  dynamic  matching  technique  may  be  feasible  for  scanned  arrays. 


(4)  Antonio  Ip.  "Radiation  and  Propagation  of  Cylindrical  Leaky  Waves  on  a  Planar 
Dielectric  Multi-Layer  Structure".  M.S.E.E..  University  of  Houston.  May  1988. 

Radiation  formulas  for  cylindrical  leaky  waves  on  a  planar  surface  are  derived  These  formulas 
predict  the  far-field  patterns  of  either  TE  or  TM  waves,  with  an  arbitrary  finite-sized  aperture 
assumed,  so  that  the  effects  of  aperture  truncation  may  be  determined.  A  narrow-beam  antenna 
obtained  by  embedding  a  printed  antenna  element  within  a  multiple-layer  dielectric  structure  is 
then  explained  from  the  existence  of  weakly  attenuated  leaky  waves  which  exist  on  this  structure. 
Analytical  methods  for  asymptotically  calculating  the  leaky-wave  propagation  constants  for  small 
beamwidths  are  presented.  Radiation  patterns  based  on  the  leaky  waves  are  compared  with  the 
exact  patterns  for  horizontal  and  vertical  electric  and  magnetic  dipoles  in  the  structure,  with 
excellent  agreement  obtained  near  the  beam  maximum,  where  the  leaky  wave  contribution  is  most 
significant. 

(5)  Sai  C.  Kwok.  "Analysis  and  Design  of  a  Traveling-Wave  Array  of  Vertical 
Monopoles  in  a  Substrate".  M.S.E.E  .  University  of  Houston.  Aug.  1988. 

This  thesis  presents  an  analysis  and  design  technique  for  a  traveling-wave  array  of  monopoles 
embedded  within  a  substrate,  fed  by  a  meandering  stripline.  A  Block  analysis  is  used  to  analyze 
the  structure,  which  is  approximated  as  a  periodically  loaded  transmission  line,  with  the  loads 
being  the  active  impedances  of  the  monopoles.  An  iterative  design  procedure  is  established  which 
determines  the  necessary  monopole  heights  and  meandering  line  section  lengths  in  order  to 
produce  a  desired  inter-element  phase  shift  and  current  amplitude  taper  for  the  set  of  monopole 
currents,  and  hence  a  desired  radiation  pattern.  Special  attention  is  devoted  to  the  efficient 
evaluation  of  the  Sommerfeld  integrals  which  are  needed  in  the  self  and  mutual  impedance 
calculations.  As  an  experimental  verification  of  the  theory,  a  ten  element  linear  array  is  built 
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and  tested,  with  good  agreement  obtained  between  theory  and  experiment. 


(6)  Arturs  E.  Dinbergs,  'Analysis  and  Design  of  an  Array  of  Electromagnetically 
Coupled  Microstrip  Dipoles  ",  M.S.E.E..  University  of  Houston,  Aug.  1988. 

This  thesis  presents  an  analysis  and  design  technique  for  an  array  of  microstrip  dipoles 
electromagnetically  coupled  to  microstrip  lines  (EMC  dipoles).  A  general  analysis  technique 
using  a  combination  of  piecewise 

sinusoidal  and  traveling-wave  basis  functions  is  used  to  analyze  the  coupling  between  elements  in 
the  array.  A  design  approach  is  formulated  in  which  coupling  coefficients  between  only  two 
isolated  EMC  dipole  elements  are  needed.  These  coefficients  are  used  together  with  a  matrix 
description  of  the  line  and  dipole  currents  in  the  array  in  order  to  account  for  the  necessary 
mutual  interactions  between  all  of  the  elements  in  the  array.  A  four  element  linear  array  is  built 
and  tested,  with  good  agreement  obtained  lielween  the  predicted  and  measured  patterns. 


(7)  Prakash 
Electromagnetically 
1988. 


Manghnani.  "Design  and  Analysis  of  an  Array  of 
Coupled  Microstrip  Patches".  M.S.E.E..  University  of  Houston,  Aug. 


In  this  thesis,  an  analysis  and  design  of  an  array  of  microstrip  patches  that  are 
electromagnetically  coupled  by  transmission  lines  embedded  within  the  substrate  is  presented.  A 
method  of  moments  formulation  is  used  to  calculate  the  coupling  coefficients  between  elements, 
using  piecewise  sinusoidal  basis  functions  on  the  transmission  lines  and  dominant-mode  cavity- 
basis  functions  on  the  patches.  An  analysis  procedure  is  then  developed  in  which  the 
coupling  coefficients  between  adjacent  elements  are  used  in  a  matrix  equation  which  then 
determines  the  patch  current  amplitudes  on  all  of  the  elements  in  the  array.  A  synthesis 
procedure  is  then  developed,  in  which  the  coupling  coefficients  are  determined  in  order  to  produce 
a  desired  set  of  patch  currents,  which  determines  the  radiation  pattern.  A  four  element  linear  array 
having  a  45u  beam  is  built  and  tested,  with  the  measured  far-field  pattern  showing  good 
agreement  with  the  predicted  pattern. 


(8)  Roger  Kranenburg  "Planar  Transmission  Line  Excitation  of  Dielectric  Resonator 
Antennas",  M.S.E.E..  University  of  Houston.  November  1988. 

An  approximate  theory  is  derived  for  all  the  modes  in  a  cylindrical  dielectric  resonator 
antenna.  The  air-dielectric  interface  is  approximated  as  a  perfect  magnetic  conductor  and 
expressions  for  the  resonant  frequency  and  radiation  from  the  equivalent  surface  magnetic  currents 
are  presented.  A  broad  beam  with  a  maximum  normal  to  the  ground  plane  is  obtained  for  the 
lowest  mode.  A  more  directive  beam  or  a  beam  with  a  single  sharp  null  is  obtained  for  the  next 
mode  depending  on  the  cylinder  geometry.  Good  agreement  is  obtained  between  theory  and 
experiment. 

An  experimental  investigation  was  also  undertaken  demonstrating  the  microstrip 
transmission  line  and  the  coplanar  waveguide  as  effective  excitation  mechanisms  for  the  cylindrical 
dielectric  resonator  antenna.  Results  showing  the  coupling  behavior  and  the  radiation  patterns  for 
various  cylinder  aspect  ratios,  dielectric  constants,  and  feed  positions  are  presented  The  results 
show  these  feed  structures  may  result  in  efficient  and  practical  antenna  systems. 
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DUAL-BAND  MICROSTRIP  ANTENNAS 
WITH  MONOLITHIC  REACTIVE  LOADING 
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Introduction;  For  a  variety  of  applications  in  telecommunica¬ 
tions  it  is  often  necessary  to  utilise  two  separate  operating 
frequencies.  Previous  work1*  has  shown  that  a  single  micro- 
strip  radiating  element  could  be  used  in  this  fashion  by 
employing  reactive  loading.  A  general  theory  has  been  formu¬ 
lated  and  preliminary  experimental  measurements  have  been 
made  for  radiators  using  a  variable-length,  short-circuited 
coaxial  stub  as  the  mechanism  for  the  loading.1  These  investi¬ 
gations  have  demonstrated  the  general  techniques  involved 
and  have  verified  the  accuracy  of  the  analysis  in  predicting  the 
resonant  frequencies  of  the  structure. 

The  presence  of  the  coaxial  load  which  protrudes  into  the 
region  behind  the  low-profile,  printed-circuit  antenna  is 
usually  not  desirable  in  practical  systems.  For  this  reason 
alternative  forms  of  loading  which  are  monolithic  in  nature 
and  thus  retain  the  low-profile  characteristics  of  the  more 
usual  single  frequency  microstrip  antenna  have  been  investi¬ 
gated.  A  set  of  antennas  have  been  designed  and  tested  which 
use  microstrip  transmission  lines  etched  on  the  same  printed- 
circuit  board  as  the  patch  radiator  for  the  reactive  loading. 
Such  a  structure  has  been  shown  to  display  dual-band  behav¬ 
iour  while  retaining  the  usual,  quite  attractive  physical  charac¬ 
teristics  of  a  microstrip  radiator. 


Design:  To  model  most  simply  the  previously  investigated 
coaxially  loaded  antennas  one  can  straightforwardly  substi¬ 
tute  a  short-circuited  microstrip  line  attached  to  the  edge  of  a 
rectangular  patch  element.  Use  of  a  coaxial  load  was  consider¬ 
ably  more  flexible,  however,  since  the  separation  between  res¬ 
onances  as  well  as  the  resonant  frequencies  themselves  could 
also  be  adjusted  by  moving  the  load  position  away  from  the 
edge.  To  emulate  this  variation  the  microstrip  line  has  a/so 
been  designed  with  a  specified  ’inset’  dimension  S  as  shown  in 
Fig.  1.  Unfortunately,  positioning  the  load  in  the  interior  of 
the  patch  results  in  yet  another  parameter  G,  the  gap  spacing 
between  the  line  and  the  radiator.  This  parameter  seems  to 
have  little  effect  on  the  upper  resonance  but  tends  to  reduce 
the  lower  resonant  frequency  as  the  gap  spacing  is  increased. 
The  effect  of  the  inset  distance,  however,  is  precisely  the  same 
as  for  the  coaxial  loads.  As  the  inset  is  increased,  the  two 
resonant  frequencies  come  closer  together,  thereby  reducing 
the  overall  separation  between  resonances. 

Results:  Several  antennas  were  fabricated  with  varying  insets 
and  gap  spacings  with  stub  iengths  predicted  by  theoretical 
calculations  to  give  dual-band  behaviour.  Two  illustrative 


Fif.  2  Impedance  of  edge-loaded,  4x6  cm  patch  antenna 


cases  arc  presented.  First  an  edge-loaded  antenna  consisting 
of  a  6  cm  by  4  cm  patch  with  a  short-circuited  microstrip 
transmission  line  of  length  L  =  4-0  cm  and  width 
W  =-  0  33  cm  (inset:  S  =  0,  G  ~  0),  was  measured.  All 
antennas  were  etched  on  a  dielectric  substrate  with  a  relative 
permittivity  of  2-17  and  a  thickness  of  0  079  cm  and  coaxially 
fed  near  the  edge  and  at  the  centre  of  the  6  cm  side.  Good 
pattern  performance  was  observed  at  each  of  the  resonant 
frequencies  (2-275  GHz  and  2  666  GHz,  respectively)  with  the 
overall  impedance  behaviour  shown  in  Fig.  2. 


Fig.  3  Impedance  of  inset-loaded,  4  x  6  cm  patch  antenna 
L  =  4  0  cm,  W  d  0-33  cm,  0  -  0-3  cm,  S  *  I  S  cm 


Secondly,  a  similar  radiator  with  L  =  4  0  cm.  W  —  0-33  cm 
and  an  inset  or  S  =  15  cm  in  from  the  edge  with  a  gap 
spacing  6  =  0-3  cm  was  investigated.  Similar  patterns  were 
measured  and  the  resulting  impedance  is  shown  in  Fig.  3.  The 
dual-band  behaviour  is  still  apparent,  but  the  separation 
between  resonances  (2-47 1  GHz  and  2-514  GHz)  is  seen  to  be 
greatly  reduced. 


4  0  8  0  '  12 
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Fig.  4  Resonant  frequency  against  line  length  for  edge-loading 
4  0  x  6  0  cm  patch 

W  =  0  33  cm,  S  =  0,  G  =  0 

A  complete  listing  of  the  measured  upper  and  lower  reson¬ 
ant  frequencies  for  various  dual-band  radiators  is  shown  in 
tabular  form  in  Table  1  and  graphically  in  Figs.  4  and  5. 
Fig.  4  shows  the  decrease  in  spacing  between  the  two  reson¬ 
ances  as  the  length  of  the  loading  stub  is  increased  for  two 
otherwise  identical  edge-loaded  patches  (4  0  x  6-0  cm  patch 
and  W  =  0-33  cm).  Fig.  5  illustrates  the  effect  of  changing  the 
gap  size  for  radiators  with  inset,  5=1-5  cm,  L  =  4-0  cm  and 

Table  1  RESONANT  FREQUENCIES  OF 
MONOLITHIC  MICROSTRIP 
ELEMENTS 
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H-  =  0-33  cm.  The  upper  resonance  is  seen  to  remain  rela¬ 
tively  stable  but  the  lower  one  tends  to  decrease  somewhat  as 
the  gap  size  is  increased. 
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Fig.  5  Resonant  frequency  against  gap  size 

4-0  x  6-0  cm  patch,  W  =  0-33  cm,  S  —  1-5  cm,  L  *=  4  0  cm 

Conclusions:  The  utilisation  of  reactive  loading  to  produce 
dual-band  behaviour  in  microstrip  antennas  has  been  imple¬ 
mented  in  a  more  practical  manner  by  the  use  of  monolithic 
microstrip  transmission  lines.  This  technique  eliminates  the 
more  cumbersome  coaxial  stubs  used  in  previous  investiga¬ 
tions. 
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Abstract — The  electromagnetic  properties  of  electrically  thick  rectan¬ 
gular  microstrip  antennas  were  Investigated  experimentally.  Antennas 
were  fabricated  with  different  patch  sizes  and  with  electrical  thicknesses 
ranging  from  0.03  to  0.23  wavelengths  in  the  dielectric  substrate.  The 
resonant  frequencies  were  measured  and  compared  to  existing  formulas. 
The  bandwidth  was  calculated  as  a  function  of  electrical  thickness  and  the 
antenna  radiation  patterns  were  measured. 

I.  Introduction 

DURING  THE  PAST  TEN  YEARS,  microstrip  antennas 
experienced  a  great  gain  in  popularity  and  have  become  a 
major  research  topic  in  both  theoretical  and  applied  electro¬ 
magnetics.  They  are  well  known  for  their  highly  desirable 
physical  characteristics  such  as  low  profile,  light  weight,  low 
cost,  ruggedness,  and  conformability.  Numerous  researchers 
have  investigated  their  basic  characteristics  and  recently 
extensive  efforts  have  also  been  devoted  to  the  design  of 
"frequency  agile,"  “polarization  agile,"  or  dual-band  mi¬ 
crostrip  antennas  [1],  [2],  [3],  Most  of  the  previous  theoretical 
and  experimental  work  has  been  carried  out  only  with 
electrically  thin  microstrip  antennas.  Recent  interest  has 
developed  in  radiators  etched  on  electrically  thick  substrates. 
This  interest  is  primarily  due  to  two  major  reasons.  First,  as 
these  antennas  are  used  for  applications  with  increasingly 
higher  operating  frequencies,  and  consequently  shorter  wave¬ 
lengths,  even  antennas  with  physically  thin  substrates  become 
thick  when  compared  to  a  wavelength.  Second,  microstrip 
antennas  have  inherently  narrow  bandwidths  and  are  normally 
not  suitable  for  broad  bandwidth  applications.  Increasing  the 
bandwidth  is  possible,  but  the  methods  used  [4],  [5],  [6] 
invariably  increase  the  volume  of  the  antenna  by  either 
extending  the  radiating  surface  or  by  increasing  the  overall 
antenna  thickness.  To  aid  in  the  design  of  broader  band 
microstrip  antennas,  a  careful  experimental  study  of  the 
resonant  frequency,  bandwidth,  and  radiation  patterns  of 
rectangular  microstrip  antennas  as  a  function  of  electrical 
thickness  of  the  substrate  was  undertaken.  The  measured 
resonant  frequencies  were  compared  to  formulas  previously 
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developed  for  predicting  the  resonant  frequency  of  electrically 
thin  rectangular  microstrip  amennas. 

II.  Experimental  Procedures 

The  microstrip  antennas  investigated  are  rectangular 
patches  with  geometry  as  illustrated  in  Fig.  1.  They  are 
fabricated  on  3M  CuClad  233  and  on  Rogers  RT/duroid  5870 
microwave  substrates.  The  CuClad  material  is  made  of  a 
poiytetrafluoroethylene  (PTFE)  woven  glass  laminate  material 
while  the  RT/duroid  material  is  made  of  a  glass  microfiber 
reinforced  PTFE  composite.  Both  substrates  have  a  nominal 
dielectric  constant  (<,)  of  2.33,  and  all  the  antennas  are  fed 
using  an  SMA  coaxial  feed.  In  this  investigation  the  feed  is 
located  a!  the  midpoint  of  the  longer  side  (x'  =  a/2)  and  at  a 
distance  from  the  edge  (y'  =  0.15  cm).  In  each  case  the 
dimension  "a”  has  been  chosen  to  be  approximately  one  and 
one-half  times  the  dimension  “ b "  with  a  10  cm  x  10  cm 
ground  plane.  Two  sets  of  regular  microstrip  antennas  have 
been  fabricated.  The  ones  in  the  first  set  have  the  same 
substrate  thickness  “h"  but  have  nine  different  patch  sizes; 
the  ones  in  the  second  set  have  the  same  patch  size  but  have 
three  different  substrate  thicknesses.  In  addition  a  so-called 
"air-dielectric"  model  radiator  has  been  fabricated  to  allow  an 
even  more  detailed  study  of  the  resonant  frequency.  Its 
geometry  models  that  of  the  regular  rectangular  microstrip 
antenna  shown  in  Fig.  1  but  with  a  substrate  whose  height  can 
be  changed  by  placing  sheets  of  styrofoam  (e,  =  1 .05)  with 
varying  thicknesses  between  the  ground  plane  and  the  radiat¬ 
ing  patch.  The  1.78  cm  x  2.67  cm  aluminum  radiating  patch 
has  a  thickness  off  =  0. 16  cm  and  is  coaxial  fed  over  a  14  cm 
x  14  cm  aluminum  ground  plane.  This  fixture  allows  the 
resonant  frequency  to  be  measured  for  a  wide  range  of 
electrical  thicknesses  using  exactly  the  same  rectangular 
radiating  patch. 

The  resonant  frequency,  impedance  and  radiation  pattern 
measurements  were  all  performed  at  the  University  of  Hous¬ 
ton-University  Park,  Applied  Electromagnetics  Laboratories, 
using  an  automated  network  analyzer  system  and  dedicaled 
computer  programs.  Accuracy  enhancement  techniques  [7J 
have  been  used  to  patially  correct  for  effective  directivity, 
effective  source  match,  and  frequency  tracking  errors  when 
taking  impedance  measurements.  The  radiation  pattern  mea¬ 
surements  were  taken  with  the  antenna  under  test  placed  inside 
an  anechoic  chamber  and  mounted  on  a  one  meter  diameter 
circular  aluminum  ground  plane. 
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III.  Experimental  Results 

A.  Impedance 

During  the  course  of  this  research,  the  input  impedances  (Z 
=  R  +  JX)  and  the  radiation  patterns  of  each  of  the  antennas 
have  been  measured.  Before  the  impedance  data  were  used  to 
determine  the  resonant  frequencies  (/,)  and  the  bandwidths 
(BW),  they  were  smoothed  in  order  to  take  out  any  residual 
ripples  or  oscillations  that  are  due  to  reflections  internal  to  the 
measurement  equipment  and  have  not  been  corrected  by  the 
accuracy  enhancement  routines.  The  values  of  admittance,  Y 
-  G  +  jB,  as  a  function  of  frequency  (/)  were  computed 
from  the  measured  impedance  versus  frequency  data  through 
the  relation  Y  -  Z" Depending  on  the  values  of  admittance, 
either  a  cubic  or  quadratic  least  squares  regression  polynomial 
is  fitted  through  the  values  of  the  admittance  versus  frequency 
curve.  From  this  fitted  polynomial,  the  smoothed  admittance 
at  each  measured  value  of  frequency  is  computed  and  then  the 
reciprocal  is  taken  to  obtain  the  smoothed  impedance  The 
actual  smoothing  operation  is  carried  out  with  the  admittance 
data  since  both  the  real  and  imaginary  parts  of  the  admittance 
are  monotonic  functions  in  the  neighborhood  of  resonance  and 
thus  result  in  a  better  polynomial  fit.  The  smoothed  curves 
follow  the  general  form  of  the  measured  traces  very  closely 
and  allow  the  true  peak  of  the  resistance  curve  to  be 
determined  more  accurately  for  resonant  frequency  measure¬ 
ments.  Fig.  2  shows  a  comparison  of  typical  smoothed  and 
measured  impedance  versus  frequency  curves  with  data  points 
taken  every  10  MHz  for  an  antenna  with  h/\,i  =  0.1 10. 
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Comparison  of  measured  and  smoothed  impedance  (I. I  cm  x 
cm  radiating  patch,  0.3175  cm  substrate,  e,  =  2.33). 
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B.  Resonant  Frequency 

Generally,  the  resonant  frequency  of  a  microstrip  antenna  is 
defined  as  the  frequency  at  which  the  reactance  is  equal  to 
zero.  For  electrically  thin  antennas,  this  point  is  also  very 
close  to  the  frequency  where  the  resistance  reaches  a  maxi¬ 
mum.  However,  in  this  investigation  many  of  the  reactance 
curves  exhibit  an  inductive  shift  due  to  the  coaxial  feed  passing 
through  the  electrically  thick  substrate  [8],  [9].  In  fact,  for  the 
thicker  antennas,  the  reactance  curve  never  passes  through 
zero  at  all  (see  Fig.  2).  For  this  reason,  the  resonant  frequency 
has  been  redefined  as  the  point  at  which  the  resistance  reaches 
a  maximum,  independent  of  the  value  of  reactance.  Further¬ 
more,  since  the  bandwidth  of  an  electrically  thin  microstrip 
antenna  is  commonly  defined  in  terms  of  the  impedance  (and 
thus  is  dependent  on  the  reactance),  an  alternate  definition  of 
bandwidth  that  is  not  affected  by  the  inductive  shift  of  the 
reactance  is  used  to  obtain  the  antenna  bandwidths  of  this 
paper.  This  last  point  will  be  discussed  in  detail  in  a  later 
section. 

Since  the  main  concern  in  the  measurement  of  the  resonant 
frequency  is  the  effect  of  the  changing  electrical  thickness  of 
the  substrate,  a  normalized  resonant  frequency  is  defined 
where  /mTm  =  /r//rp,  and/,0  is  the  zeroth-ordcr  prediction  of 
the  resonant  frequency.  This  approximation  for  fr0  assumes 
that  the  antenna  thickness  is  infinitesimally  thin  and  that  b  is 
equal  to  \/2.  Then  knowing  that  \d  =  c/(fr0-JTr)  =  2b,  f, 0 
can  be  computed.  For  the  units  of  c  in  m/s  and  those  of  b  in 
cm. 


/,o=  15/(6  v/T,)  GHz.  (I) 

Table  I  shows  the  measured  resonant  frequency,  zeroth  order 
prediction,  physical  dimensions,  and  electrical  thickness  of 
each  antenna.  Fig.  3  shows  a  plot  of  the  normalized  resonant 
frequency  plotted  as  a  function  of  electrical  thickness  for  the 
nine  antennas  etched  on  the  same  thickness  of  substrate  and  for 
the  three  identically  sized  antennas  on  different  substrate 
thicknesses. 

Table  11  shows  the  measured  resonant  frequencies  arxl  the 
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TABLE  I 

MEASURED  AND  PREDICTED  ANTENNA  RESONANT  FREQUENCIES 


8 

b 

h 

flex'd 

James 

Hammerstad 

'v'Ni 

(cm) 

(cm) 

(cm) 

(OHz) 

(GHz) 

(GHz) 

5.7 

3.8 

0.3175 

2.31 

2.30 

2.30 

0.037 

4.55 

3.05 

0.3175 

2.89 

2  79 

2.90 

0.047 

2.95 

1.95 

0.3175 

4.24 

4.11 

4.34 

0.068 

1.95 

1.3 

0.3175 

5.84 

5.70 

6.12 

0.094 

1.7 

1.1 

0.3175 

6.80 

64? 

7.01 

0.110* 

1.4 

0.9 

0.3175 

7.70 

7.46 

3.19 

0125 

12 

0.8 

0.3175 

8.27 

0.13 

9.01 

0.141 

1.05 

0.7 

0.3175 

9.14 

8.09 

9.97 

0.148 

0.9 

0.6 

0.3175 

10.25 

9.02 

11.18 

0.166 

1.7 

1.1 

0.1524 

7.87 

7.46 

7.84 

0.061 

17 

1.1 

03175 

6.80 

6.47 

7.01 

o.no* 

1.7 

1.1 

09525 

4.73 

4  32 

5.2? 

0229 

1  hex  two  are  (Ns  same  antennas. 


Electrical  Thickness  (wavelengths) 

Fig.  3.  Normalized  antenna  resonant  frequency  versus  electrical  thickness. 


TABLE  II 

Resonant  frequency  versus  h  for  air-Dielectrjc  Fixture  n  .78 

cm  x  2.67  cm  RADIATING  PATCH) 


h 

Mess'd 

James 

Hammerstad 

h/\ 

Icm) 

15Hz) 

(GHz) 

(GHz) 

0.64 

5.14 

4  54 

575 

0.112 

0.79 

5.12 

4.18 

542 

0138 

0.99 

4.33 

3.78 

5.07 

0.146 

1.19 

4.2? 

346 

4.77 

0.174 

1.44 

3.32 

313 

446 

0.163 

1.64 

3.06 

2.91 

425 

0.171 

2.04 

2.56 

2.55 

3.91 

0.178 

2.34 

2.29 

234 

3.70 

0.®3 

corresponding  electrical  thicknesses  of  the  air  dielectric 
fixture,  while  the  circles  in  Fig.  4  represent  the  same  data  in 
graphical  form.  It  is  evident  from  Figs.  3  and  4  that  the 
resonant  frequency  indeed  decreases  as  the  antennas  become 
electrically  thicker  as  has  been  shown  in  previously  published 
results  [10],  [1 1],  It  is  perhaps  unexpected,  however,  that  this 
trend  continues  even  to  thicknesses  approaching  one  quarter 
wavelength. 

C.  Bandwidth 

The  percent  bandwidth  of  the  antennas  was  determined  from 
the  impedance  data.  For  ease  of  notation  the  term  bandwidth 


refers  to  percent  bandwidth  unless  otherwise  specified.  Band¬ 
width  is  normally  defined  as 

percent  BW  =  [(/^ -/,,)//,]  100  percent  (2) 

where  /,  is  the  resonant  frequency,  while  fr2  and  are  the 
frequencies  between  which  the  magnitude  of  the  reflection 
coefficient  of  the  antenna  is  less  than  or  equal  to  1/3  (which 
corresponds  to  a  voltage  standing-wave  ratio  (VSWR)  <  2.0). 
However,  this  definition  is  found  not  to  be  directly  applicable 


770 


IEEE  TRANSACTIONS  ON  ANTENNAS  AND  PROPAGATION.  VOL  AP-34.  NO.  6.  JUNE  1986 


a 

D 
C T 


o Moasurod 
A  Jamas  ot  al . 
v  Hammarstad 


1 


Substrata  Thickness  (cm) 


Fig.  4,  Air-dielectric  fixture  resonant  frequencies  versus  substrate  thickness 
(1.78  cm  x  2.67  cm  radiating  patch,  t,  *  1.05). 


to  the  experimental  data  because  of  the  inductive  shift.  Thus, 
an  alternate  definition  is  found  in  order  to  determine  the 
bandwidth  of  the  test  antennas. 

The  case  where  the  impedance  at  resonance  is  purely 
resistive  (Z,„  =  /?„,,)  may  be  represented  by  a  parallel  RLC 
circuit,  and  an  analytical  expression  for  the  input  impedance  in 
terms  of  the  antenna  2-factor,  Rm„  and  /,  may  be  written  as 


1 

r  Qf' 

- fQ  1 

n  j 

max 

*™,/J 

_i_+r 

Of 

of,  y 

^nux/r  . 

Using  (3),  it  was  found  that  R  and  the  magnitude  of  Z  had  a 
definite  relationship  to  that  is  independent  of  the  other 
parameters.  Namely, 

|Z(/ri)l  =  \Z(f,i)\  =  0.8l8Rm„,  (4) 


It  should  be  noted,  however,  that  no  attempt  has  been  made 
in  this  investigation  to  actually  match  the  antennas  to  a 
standard  50  Q  transmission  line  by  a  technique  such  as  moving 
the  feed  point  away  from  the  edge.  Thus  the  bandwidths 
reported  in  Table  III  and  Fig.  5  are,  in  effect,  projected  ones 
that  might  be  obtainable  under  the  more  usual  definition  of 
bandwidth.  These  values  are  most  useful  for  comparison 
purposes  to  characterize  the  dependence  of  the  bandwidth  on 
the  various  antenna  parameters. 

D.  Radiation  Pattern 

The  radiation  pattern  for  each  element  was  measured  at  its 
resonant  frequency  in  both  the  £-planc  and  the  //-plane  The 
experimental  results  show  that  the  radiation  patterns  of 
electrically  thick  antennas  are  very  similar  to  those  of 
electrically  thin  antennas.  The  //-plane  patterns  remain 
virtually  unchanged  while  the  £-plane  ones  begin  to  show 
some  asymmetries  only  for  the  larger  substrate  thicknesses. 
Fig.  6  shows  the  radiation  patterns  of  a  representative  thicker 
antenna  {h/\d  =  0.1 10,  taken  at  6.8  GHz  in  0.5’  steps)  and  it 
is  typical  of  the  remaining  antennas. 

IV.  Comparison  of  Measurfd  and  Predicted  Resonant 
Frequencies 

The  measured  resonant  frequencies  can  be  compared  to  the 
predicted  resonant  frequencies  as  a  mutual  check  of  the 
experimental  data  and  of  the  validity  of  the  theories.  The 
equations  for  the  resonant  frequency  proposed  by  Hammerstad 
[13]  and  by  James,  Hall,  and  Wood  [11]  are  used  for  these 
comparisons.  Both  methods  share  the  concept  of  an  effective 
dielectric  constant  (<tfr)  given  by  [14) 


,  (e,+  1)  (r,-  1)(1  +  \0h/w)-'n 

fett  M  =  — —  + - - - - 


O) 


where  w  is  a  variable  and  can  be  either  the  patch  dimension 
"a"  or  "b."  Hammerstad  gives  a  predicted 


and 


*(/„)  =  /?(/,  3)  =  0.670/?m„.  (5) 

For  the  special  case  of  a  prescribed  VSWR  =  2  .0,  this  method 
is  equivalent  to  the  previously  derived  analytical  expression 
[12] 

BW  =  (VSWR -  l)/(2\/VSWR).  (6) 

For  electrically  thick  radiators  with  the  associated  large 
inductive  shift  in  the  impedance,  the  VSWR  may  not  be  below 
2.0  for  any  range  of  frequencies.  Using  only  the  resistance 
data  a  projected  bandwidth  can  be  calculated,  however, 
following  this  resonant  circuit  model  by  locating  the  frequen¬ 
cies  where  the  resistance  is  equal  to  0.670  times  the  value  of 
/?„*,.  Table  111  shows  the  bandwidths  so  obtained  arranged  in 
order  of  increasing  electrical  thickness  along  with  the  corres¬ 
ponding  /,i  and  f, 2,  while  Fig.  5  shows  the  same  bandwidth 
data  in  graphical  form.  It  is  seen  that  values  of  bandwidth  on 
the  order  of  20  percent  may  be  achieved  using  electrically 
thick  antennas. 


2(b  +  2Ab)  -J7,  ' 

,  0.412/i(ftfir  (c)  +  0.3)(n//t  +  0.264) 

At>  = - 

(e,rr  (°)  ~  0.258)(a/h  +  0.8) 

James  et  al.  give  a  predicted 


(8) 

(9) 


fr  = 


frOtr 


where 

6  =  (/i/*)0.882h 


*^Fcrr  (ujtcfr  (b)(  1  +  5) 
0" 


(10) 


[ 


"0.164(fr 

_  7f 

((,+  1)[0.758  +  In  (b/h  +  1 .88)] 


on 


The  predicted  resonant  frequencies  obtained  using  these  two 
methods  are  shown  in  Tables  I  and  II  for  the  test  antennas  and 
for  the  air-dielectric  fixture. 
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TABLE  III 

Projected  Values  or  bandwidth  and  cor resi*on ding  f,,  and  /•;, 


h/A* 

frl 

fr2 

BUJ 

(GHzl 

(GHz) 

(X) 

0.037 

2.280 

2.352 

3.117 

0.047 

2.834 

2.952 

4.083 

0.061 

7.632 

8.152 

6.607 

0.068 

4.120 

4.396 

6.509 

0.094 

5.632 

6.140 

8.699 

0.110 

6.494 

7.272 

11.441 

0.125 

7.314 

8.468 

14.987 

0141 

7.848 

9.024 

14.220 

0.148 

8.380 

10.560 

23.850 

0.166 

9.442 

11.560 

20660 

0.229 

4.320 

5.180 

18180 

Electricol  Thickness  (wavelengths) 
Fig.  5.  Projected  anlenna  bandwidth  versus  electrical  thickness. 


Fig.  6.  Antenna  radiation  pattern  (1.1  cm  x  1.7  cm  radiating  patch.  0  3175 
cm  substrate,  «,  *  2.33). 
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Fig.  7.  Comparison  of  normalized  antenna  resonant  frequencies  versus 
electrical  thickness  ( h  *  0.3175  cm.  e,  =  2.33). 

The  predicted  resonant  frequencies  are  normalized  to  the 
zeroth  order  predictions  and  are  compared  to  the  measured, 
normalized  resonant  frequencies.  Fig.  7  shows  this  compari¬ 
son  for  the  nine  different  patches  on  the  substrates  with 
identical  physical  thickness,  while  Fig.  4  shows  the  predicted 
values  of  resonant  frequency  for  the  air-dielectric  fixture 
compared  to  the  actual  measured  values.  It  should  be  noted 
that  the  theoretical  data  in  Fig.  7  are  presented  as  discrete 
points  rather  than  a  continuous  curve  so  that  the  theory  can 
correspond  to  the  exact  cases  of  the  experimental  cases,  some 
of  which  have  slightly  varying  values  of  a/b  ratio  and  of 
dielectric  constant.  This  use  of  individual  points  causes  the 
data  to  no  longer  be  aligned  along  smooth  curves.  In  both 
cases  it  is  clear  that  the  theories  follow  the  trend  of  the 
experimental  data  quite  well  even  for  electrically  thick 
substrates.  In  fact,  based  on  the  information  presented  here 
and  on  additional  research  data  [15],  some  general  observa¬ 
tions  may  be  made  concerning  these  two  methods.  Both 
predict  the  resonant  frequency  very  closely  for  electrically  thin 
rectangular  microstrip  antennas,  but  as  h  becomes  greater  than 
0.1  James  el  al.  give  consistently  better  predictions  than 
the  method  by  Hammerstad.  Specifically,  James  el  al.  usually 
predict  values  approximately  4  percent  lower  than  the  mea¬ 
sured  resonant  frequency  while  Hammerstad  predicts  values 
around  8  percent  higher  than  the  measured  resonant  fre¬ 
quency.  It  should  be  noted,  however,  that  these  two  theories 
are  not  intended  for  use  with  electrically  thick  substrates.  For 
h  ^  0.1  X,;,  the  measured  resonant  frequency  is  very  nearly 
the  mean  of  the  predicted  resonant  frequencies  from  the  two 
different  formulas.  A  third  algebraic  formula  for  the  resonant 
frequency  has  been  proposed  by  Sengupta  [16],  but  as  the 
author  states,  it  only  applies  to  electrically  thin  structures  and 
does  not  predict  the  proper  behavior  for  the  thicker  substrates 
measured  in  this  investigation. 

V.  Conclusion 

The  effect  of  varying  the  electrical  thickness  for  rectangular 
microstrip  antennas  has  been  investigated  experimentally 
during  the  course  of  this  research.  In  addition,  an  air-dielcctric 
model  radiator  has  been  fabricated  with  a  single  patch  size  and 
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a  variable  substrate  thickness.  The  resonant  frequency,  band¬ 
width,  and  radiation  pattern  have  been  measured  over  a  range 
of  substrate  thickness  from  0.03  to  0.23  of  a  wavelength  in  the 
dielectric. 

The  resonant  frequency  of  a  rectangular  microstrip  antenna 
was  found  to  decrease  as  a  function  of  electrical  thickness.  The 
validity  of  this  finding  is  confirmed  by  previously  published 
results  and  by  the  existing  theories  of  Hammerstad  and  of 
James  et  af.  The  predicted  resonant  frequencies  are  all  very 
close  to  the  measured  resonant  frequency  for  electrically  thin 
substrates.  As  the  electrical  thickness  is  increased,  the  theories 
due  to  Hammerstad  and  to  James  et  al.  generally  predict 
values  that  are  approximately  8  percent  high  and  4  percent  low 
respectively  when  compared  to  the  actual  values.  The  impe¬ 
dance  of  the  thicker  antennas  is  characterized  by  an  inductive 
shift  in  the  reactance  away  from  zero  at  resonance. 

The  projected  bandwidth  of  the  unmatched  antennas  was 
calculated  to  determine  the  effect  of  substrate  thickness  on  this 
characteristic  as  well.  The  resulting  data  show  that  bandwidths 
as  high  as  20  percent  could  be  achieved  by  simply  using 
electrically  thick  substrates.  Finally,  the  radiation  patterns  of 
electrically  thick  antennas  were  seen  to  be  very  similar  to  those 
of  the  more  usual  thin  ones.  Overall  it  has  been  shown  that 
electrically  thick  rectangular  microstrip  antennas  retain  most 
of  the  desirable  electrical  characteristics  of  thinner  ones  and 
may  be  utilized  for  broad-band  applications  assuming  a 
reactive  network  is  used  for  impedance  matching. 
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1  Introduction 

Microstrip  elements  have  proven  to  be  very  versatile  and  have  found 
uses  in  a  variety  of  applications  from  satellite  arrays  to  automobile 
anti- collision  warning  systems.  This  versatility  stems  mainly  from 
their  structural  simplicity,  ease  of  construction,  ruggedness,  con- 
formability,  low  profile,  and  low  cost.  In  addition  to  these  features, 
microstrip  elements  are  capable  of  producing  a  variety  of  different 
patterns  and  polarizations  (1],  [2],  (3),  [4],  (5),  [6],  [7],  (8).  Most 
research  efforts,  however,  have  been  limited  to  basic  applications 
using  only  a  single,  unperturbed  mode  of  a  rectangular  or  circu¬ 
lar  element.  Additional  applications  are  possible  through  the  use 
of  reactive  loading.  These  include  multiband,  frequency-agile  and 
polarization-diverse,  adaptable-impedance,  and  multimode-pattern 
antennas.  Some  experimental  work  has  been  done  in  these  areas  (2], 
(3),  (4),  but  very  little  theoretical  work  with  the  exception  of  (3),  (5], 
(7),  [8],  [9),  [10]  has  been  done.  In  this  paper,  a  very  simple  theory  is 
presented  which  will  guide  a  designer  who  wishes  to  take  advantage 
of  the  added  flexibility  that  reactive  loading  provides.  In  addition, 
several  applications  are  discussed  with  corresponding  theoretical  and 
experimental  results. 

2  Theory 

A  simple  theory  that  makes  use  of  a  graphical  interpretation  of  the 
characteristic  equation  for  the  resonant  frequency  of  a  loaded  mi¬ 
crostrip  element  is  presented  in  this  section.  By  a  "reactively  loaded 
element”  we  mean  a  simple  microstrip  antenna  to  which  one  or  more, 
two-terminal  lumped,  lossless,  linear  load  impedances  are  attached — 
one  terminal  to  the  patch  and  the  other  to  the  ground  plane.  Also 
included  in  the  set  of  reactive  loads  are  sections  of  transmission  lines 
attached  at  one  or  more  points  on  the  patch.  The  simplest  form  of 
the  theory  results  when  the  original,  unloaded  antenna  has  one  of 
the  separable  shapes  which  fall  into  four  classes:  rectangular,  circu¬ 
lar,  elliptical,  and  parabolic  In  this  paper,  wc  will  deal  only  with 
rectangular  elements. 

Foundation — the  cavity  model —  The  foundation  upon  which 
the  theory  is  built  is  the  so-called  cavity  model  [11)  as  modified  in 
[12].  This  model  is  best  suited  for  understanding  loaded  elements 
because  it  can  reconstruct  the  localized  field  in  the  vicinity  of  the 
loads,  yet  it  is  simple  enough  that  basic  physical  principles  useful  in 
design  can  be  developed  and  is  also  surprisingly  accurate.  It  is  based 
on  two  observations. 

1.  Tin*  real  part  of  tin*  resonant  fn-qiiemy  «»f  a  microstrip  antenna 
is  approximately  the  same  as  that  of  a  cavity  formed  by  intro 
during  a  magnetic  conductor  between  the  patch  edge  and  the 
ground  plane 


2.  The  quality  factor  of  the  microstrip  element  can  be  determined 
with  good  accuracy  using  the  resonant  canity  field  rather  than 
the  exact  antenna  field.  Once  the  quality  factor  Q  is  deter¬ 
mined,  an  artificially  lossy  dielectric  with  a  loss  tangent  equal 
to  the  reciprocal  of  Q  can  be  imagined  to  be  filling  the  cavity. 
The  field  within  su:!.  a  cavity  approximates  the  field  between 
the  patch  and  ground  plane  of  the  original  antenna. 

Thus,  the  most  fundamental  step  in  determining  the  properties  of 
a  loaded  antenna  is  to  determine  the  resonant  frequency  and  mode 
distribution  within  the  corresponding  loaded  cai/ily. 

If  the  reactive  load  is  imagined  as  being  contained  within  a  small 
cylinder  extending  between  the  ground  plane  and  the  patch,  then 
Huygen’s  principle  can  be  used  to  replace  the  load  by  equivalent 
magnetic  and  electric  surface  currents  impressed  in  a  now  empty  cav¬ 
ity.  Since  microstrip  elements  are  thin,  the  magnetic  su  face  current 
is  primarily  a  horizontal  circulating  loop  of  current  and  the  electric 
surface  current  is  vertically  directed.  By  making  suitable  transforma¬ 
tions  [13],  one  can  approximate  the  field  due  to  the  magnetic  current 
as  a  constant  times  the  field  produced  by  the  electric  current  acting 
alone  [10].  The  effect  of  the  equivalent  magnetic  current  is  usually 
negligible  anyway  and  thus  the  effect  of  the  load  can  be  accounted 
for  by  impressing  an  equivalent  electric  current  on  the  surface  of 
the  small  cylinder  enclosing  the  load  in  an  empty  cavity.  The  total 
equivalent  current  is  just  the  total  load  current  /l- 

Characteristic  equation  for  the  resonant  frequency —  There 
are  two  expressions  for  the  load  voltage  that  one  can  write  which 
when  equated,  yield  the  characteristic  equation  for  the  resonant  fre¬ 
quency  of  the  loaded  cavity.  One  expression  is  that  the  load  voltage 
Vi,  is  the  driving-point  impedance  of  the  empty  cavity  jX\\  times 
the  load  current  7l,  where  X\\  is  the  driving-point  reactance.  The 
second  expression  is  the  defining  relationship  for  the  load  itself, 

Vi  =  -J'A'iA. 

where  A'l  is  the  load  reactance.  Combining  these  two  equations 
yields  the  characteristic  equation, 

A'l  +  A' u  =  0. 

A  very  useful  graphical  intci  pretation  of  this  characteristic  equation 
can  be  obtained  if  the  microstrip  cavity  is  modeled  as  seen  from 
the  load  terminal  pair  as  illustrated  in  Fig.  1.  The  resonant  circuit 
shown  in  the  figure  represents  the  contribution  to  the  driving-point 
impedance  from  the  resonant  mode  of  the  empty  cavity  whose  res¬ 
onant  frequency  is  nearel  that  of  the  corresponding  loaded  cavity 
mode  of  interest.  The  reactive  contribution  due  to  all  other  modes 
of  the  empty  cavity  in  the  Foster  expansion  of  the  driving-point 
impedance  is  what  we  call  the  "feed  reactance,”  A'r-  We  have  de¬ 
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Figure  I.  Circuit  model  for  the  input  impedance  of  a  microstrip  cavity. 


picted  this  as  a  simple  inductive  reactance  although  it  does  not  re¬ 
ally  correspond  to  a  single  inductor.  In  most  cases,  particularly  cases 
in  which  degenerate-  or  nearly  degcnoratc-modc  pairs  of  the  empty 
cavity  do  not  exist,  the  feed  reactance  is  indeed  positive  and  slowly 
varying  over  the  band  of  both  the  loaded-cavity  and  empty-cavity 
modes  of  interest.  Regardless  of  how  it  is  depicted,  it  is  computed 
using  the  full  modal  analysis  of  the  cavity  and  as  such,  it  is  precise. 
Using  this  model,  the  characteristic  equation  becomes 

-  (Al,  +  Af)  =  XT.  (1) 

An  even  simpler  way  to  obtain  this  characteristic  equation  is  to  apply 
Kirchoff’s  voltage  law  to  the  circuit  illustrated  in  Pig.  1  with  the  load 
reactance  attached.  Since  a  resonant  mode  of  the  loaded  element  is 
a  source-free  mode,  the  sum  of  the  voltage  drops  around  the  loop 
must  vanish  implying  that 

(Al  +  Af  +  Ar)/L  =  0. 

This  yields  the  characteristic  equation  (1)  since  Ii  must  not  be  zero. 

In  order  to  develop  important  design  concepts,  it  is  useful  to  solve 
equation  (1)  graphically.  Such  a  graphical  solution  is  illustrated  in 
Fig.  2(a).  The  simple  case  a  short-circuit  load  A'l  =■  0  is  consid¬ 
ered  first.  The  intersection  of  the  dashed  line  in  Fig.  2,  representing 
— Af  alone,  with  the  solid  line  representing  A'r  marks  the  resonant 
frequency  of  the  loaded  cavity. 

Several  design  principles  emerge  from  this  figure,  particularly  in  ref¬ 
erence  to  the  rectangular  element  illustrated  in  Fig.  3. 


1.  It  is  immediately  clear  that  the  introduction  of  a  short  circuit 
will  increase  the  resonant  frequency  of  the  cavity. 

2.  The  resonant  reactance  corresponding  say  to  the  m  =  0,  n  =  1 
mode  of  the  rectangular  clement  (the  mode  most  commonly 
used)  is 


Ar  =  ~J/it 


2c 7  cosz(ry/6) 
ah  u1  -  ’ 


where  the  patch  dimensions  arc  a  x  b  x  t,  c  is  the  speed  of  light 
in  fhe  Cfletfy  dielectric,  koi  ~  vfb  is  the  resonant  wavenum¬ 
ber  of  the  (0.  l)-mode,  and  y  is  the  y  coordinate  of  the  load 
position.  When  y  approaches  either  edge,  y  =  0  or  y  =  i, 
the  resonant  reactance  giv^n  in  equation  (2)  becomes  station¬ 
ary.  As  y  — •  bj 2,  the  resonant  reactance  approaches  zero  In 
graphical  terms,  the  reactance  curve  begins  as  a  very  sender 
departure,  from  zero  for  load  points  near  y  —  h/7  and  then  fat 
tens  out  as  y  moves  away  from  the  centerline.  For  load  points 
near  the  edge,  the  X,  curve  changes  little*  This  trend  is  illus¬ 
trated  in  Fig.  2(b)  by  the  sequence  cf  Ihrcf*  XT  curves  labeled 
1-3  corresponding  to  feed  points  1-3 


Frequency 


Figure  2.  (a)  Plot  of  Xr>  -  (A/  +  X]_)  versus  frequency;  (b)  Plot  of  Xr,  -  X f 
for  three  different  feed  locations. 


This  observation  illustrates  how  load  placement  affects  the  res¬ 
onant  frequency.  A  short  placed  near  the  centerline  changes 
the  resonant  frequency  little.  As  the  load  is  moved  away  from 
the  centerline  in  cither  direction  toward  an  edge  the  resonant 
frequency  increases  rapidly  at  first  and  then  more  slowly  as  the 
load  point  nears  the  edge. 

This  graphical  interpretation,  which  will  be  found  very  useful 
in  later  applications,  is  consistent  with  the  observation  that 
the  load  will  have  more  effect  on  the  resonant  frequency  when 
placed  where  the  empty-cavity  field  is  strongest  and  no  effect 
where  the  cavity  field  is  zero. 


Figure  3  Rectangular  element  with  three  different  load  locattons. 
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3.  The  preceding  observation  about  the  dependency  of  resonant 
frequency  on  load  position  made  the  tacit  assumption  that  the 
feed  reactance  is  independent  of  of  feed  position.  Nothing  could 
be  further  from  the  truth.  In  fact,  the  feed  reactance  can  vary 
by  a  factor  of  three  or  more  as  the  feed  point  is  moved  from 
place  to  place  on  the  patch  (13J.  The  feed  reactance  typically 
increases  as  the  feed  point  is  moved  near  an  edge.  The  increase 
is  rapid  when  the  feed  point  is  very  near  the  edge-  This  feed 
reactance  is  predicted  fairly  well  by  the  cavity  model,  particu¬ 
larly  for  thinner  elements,  except  for  feed  points  right  on  the 
edge.  There,  the  cavity  model  substantially  overestimates  the 
feed  reactance  because  the  equivalent  electric  currents  on  the 
magnetic  wall  that  would  make  the  cavity  model  precise  have 
been  neglected.  These  equivalent  currents  tend  to  oppose  the 
impressed  source  current  (in  this  application,  the  load  current) 
and  ignoring  them  results  in  an  increased  inductive  reactance. 
The  effect  of  this  variation  is  that  as  the  load  point  moves 
toward  the  edge,  the  -Xf  curve  moves  downward.  This  is  il¬ 
lustrated  by  the  sequence  of  three  curves  in  Fig.  2(b).  Thus, 
it  is  expected  that  the  resonant  frequency  should  begin  at  the 
resonant  frequency  of  the  empty  cavity  when  y  is  on  the  cen¬ 
terline,  it  should  increase  rapidly  at  first  due  to  a  widening  Xr 
curve,  and  finally  decrease  somewhat  due  to  a  rapidly  increas¬ 
ing  feed  reactance  for  load  points  near  the  edge.  This  is,  in 
fact,  exactly  what  is  observed  by  measurements. 


Determination  of  feed  reactance —  In  order  to  graph  the  -X\ 
curve,  the  input  reactance  of  a  microstrip  cavity  must  be  computed. 
This  can  be  obtained  by  using  the  Green’s  function, 


unnormalized  mode  function  ^(r)  of  the  loaded  cavity  in  terms  of 
the  modcr  of  the  empty  cavity  is 


'J'(r)  = 

A 


tfrmn(r)d>mn(r,) 
k*  - 


(3) 


where  k  is  the  resonant  wavenumber  of  the  resonant  mode  of  the 
loaded  cavity  and  rj  is  the  position  vector  for  locating  the  center  of 
the  load.  No  integration  is  used  in  this  expression  since  it  is  not  used 
to  compute  the  field  in  the  region  of  the  soutcc  (i.c.,  load)  current. 
The  evaluation  of  this  series  is  also  facilitated  by  the  use  of  an  SEM 
expansion  with  a  quasistatic  correction. 

It  is  apparent  from  equation  (3)  that  the  closer  that  the  resonant 
wavenumber  k  is  to  a  given  associated  with  a  particular  mode  of 
the  empty  cavity,  the  closer  the  loaded-cavity  mode  will  resemble  the 
V'mn (r)  empty- cavity  mode.  More  will  be  said  about  this  for  specific 
applications  discussed  in  the  next  section. 

Quality  factor,  pattern,  and  impedance —  Once  the  modal  dis¬ 
tribution  has  been  obtained,  the  pattern  can  be  approximated  as  the 
field  due  to  a  filamentary  magnetic  loop  of  current  coinciding  with 
the  edge  of  the  patch.  The  current  at  a  given  point  r  on  the  loop 
is  S'(r)  and  is  impressed  over  a  grounded  dielectric  slab  with  the 
patch  removed.  Once  the  pattern  is  obtained,  the  power  radiated  in 
the  space  and  surface  waves  can  be  computed.  The  average  energy 
stored  in  the  electric  field  is  computed  using 


G(rl*-')  =  M* 

A 


~  kln 


where  A  is  the  set  of  all  pairs  of  non-negative  integers,  (m,n),  and 
^mn(r)  is  the  mnlh  mode  of  the  cavity  satisfying 

V30mn(r)  =  n^mn(r),  re  S, 


n(r)  •  n  s  0,  r  e  dS. 


The  source  and  observation  position  vectors  are  r  =  xx  +  yy  and 
r>  -  x'x  +  y'y,  S  is  the  set  of  points  making  up  the  patch,  dS  is  the 
set  of  points  on  the  patch  edge,  and  n  is  the  outward-pointing  unit 
normal  to  the  patch  edge.  In  the  case  of  rectangular  elements, 


^nn(r) 


2 

(  mrz) 

1  COS  I 

(™y\ 

\A^  +  ^mo)(l  + 

v  a  y 

l  b  ) 

The  scries  for  G(r|r')  diverges  as  r  —  r'.  This  is  because  the  in¬ 
ductive  reactance  associated  with  a  filament  is  infinite.  In  order  to 
obtain  the  input  reactance  associated  with  the  actual  load  current, 
one  must  integrate  this  Green’s  function  over  the  load  current  distri¬ 
bution.  It  is  assumed  that  the  load  current  is  distributed  uniformly 
on  a  circular  cylinder  of  diameter  d.  The  “location”  of  the  load  is 
the  center  of  the  cylinder’s  cross  section.  The  required  integration 
needs  to  be  computed  efficiently  since  the  nonlinear  characteristic 
equation  for  the  resonant  frequency  must  be  solved  numerically,  or 
enough  points  plotted  to  solve  it  graphically.  A  very  efficient  SEM 
expansion  of  the  cavity  field  with  a  quasistatic  correction  has  been 
developed  [10j,  the  details  of  which  are  too  lengthy  for  presentation 
here.  Once  th^*  integration  has  been  carried  out,  the  result,  when 
multiplied  by  -jt  yields  the  total  reactance  Xu  at  the-  load  point 
from  which  the  resonant  reactance  can  be  subtracted  leaving  X\. 

Determination  of  the  modal  distribution —  Once  the  resonant 
frequency  is  determined,  the  task  turns  to  computing  the  modal  dis¬ 
tribution  of  the  loaded  cavity.  For  the  case  of  a  single  load,  the 


Where  Wl  is  the  average  electric  energy  stored  by  the  load.  From  the 
radiated  power  Pr  and  the  stored  electric  energy  Wt,  the  radiative 
quality  factor 


4r/0Wc 

Qr~  Pr 

is  computed  where  /o  is  the  resonant  frequency  of  the  loaded  an 
tenna.  In  addition  to  the  radiative  quality  factor,  the  quality  factors 
associated  with  the  power  consumed  in  heating  the  dielectric  and 
metal  cladding  are  Qd  -  Iff  and  Qm  =  tf A,  where  6  is  the  actual 
loss  tangent  of  antenna  dielectric,  and  A  is  the  skin  depth  in  the 
metal  cladding.  From  these  quality  factors,  the  quality  factor  of  the 
antenna  is  determined  using 


To  compute  the  input  impedance  (at  least  for  the  case  of  loaded 
elements  in  which  only  one  resonant  mode  contributes  significantly 
in  the  band  of  interest),  one  can  use  the  circuit  model  illustrated  in 
Fig.  4.  The  admittance  of  the  resonant  circuit  is 


4+*  ((£)’- Ol- 

The  conductance  G  can  be  determined  by  using  the  relation 

IKI’G  =  P, 


where  voltage  \\  across  the  resonant  circuit  is 
Vx  -  -<¥(rf). 


rf  being  the  location  of  the  feed  of  the  loaded  clement.  The  power 
P  is  the  total  power  absorbed  by  the  antenna  element  corresponding 


i_t  L 
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X(f) 


Figure  4.  Circuit  model  of  the  loaded  microstrip  antenna. 


the  resonant  mode  sEf (▼*)  computed  earlier  to  determine  the  antenna 
quality  factor.  Note  that  Vt  is  not  the  feed  voltage  since  'l'(rf)  does 
not  include  the  effect  of  the  feed  reactance.  (The  feed  reactance  of 
the  empty  cavity  at  the  load  point,  of  course,  does  affect  the  Vt.  The 
feed  reactance  referred  to  here  is  the  feed  reactance  at  the  actual  feed 
point  of  the  loaded  cavity). 

The  feed  reactance  can  be  modeled  as 


A'f  =  A'ro  +  (/  -  fo)Xi\. 


Let 


f+  (1  +  2C?)/o'  =  0  2  q)/o' 

The  coefficients  A'ro  and  An  can  be  approximated  as 

Xa  *jW+)  +  X(/.))f^, 

and 

Xn  *|  [*(/+)-*(/-) +  221], 

where  R  =  1  (G.  What  remains  is  to  determine  the  input  reactance 
of  the  ideal  cavity  at  the  upper  and  lower  frequencies  /+  and  /_ . 

To  do  this,  one  simply  views  the  input  reactance  at  the  feed  port  as 
the  input  reactance  of  a  loaded  two- port  network.  The  z  parameters 
of  the  two-port  network  arc  the  z  parameters  of  the  empty  cavity 
with  port  2  located  at  the  feed  and  port  1  located  at  the  load.  The 
mutual  reactance  between  the  ports  with  respective  position  vectors 
r i  and  r7  is  just 

X\i  =  JlG(r,|r,). 


The  input  reactance  is  then 


X  =  X77  - 


X\\  +  A‘l’ 


where  Xu  is  the  input  reactance  of  the  empty  cavity  at  r2  and  is 
computed  in  the  same  way  as  An- 

Extension  to  multiple  loads —  All  of  the  preceding  was  geared 
to  a  sinp/c  load.  This  was  done,  initially,  because  very  useful  design 
principles  could  be  developed.  Although  the  nice  graphic; I  iuterpre 
lation  of  the  characteristic  equation  will  be  lost,  the  results  can  be 
easily  modified  to  accommodate  the  case  of  multiple  loads  [10]  This 
extension  is  most  valuable  in  a  number  of  applications  since  symnirl 
ric  pairs  of  loads  can  then  be  used  to  suppress  large  cross- polarised 
fields  that  would  otherwise  be  generated  by  a  single  load. 


3  Applications 


A  number  of  applications  of  loaded  microstrip  antennas  are  discussed 
in  this  section.  Many  were  discovered  by  purely  experimental  inves¬ 
tigations,  but  all  benefit  from  the  simple  design  principles  and  the 
theory  developed  in  the  preceding  section. 

Frequency-agile  antennas —  One  of  the  first  applications  recog¬ 
nized  [3]  was  that  one  could  change  the  resonant  frequency  of  mi- 
ciostrip  antennas  by  shorting  the  patch  at  one  or  more  locations 
What  is  just  as  important  is  that  one  could  vary  the  resonant  fre¬ 
quency  over  a  wide  range  without  significantly  changing  the  radiation 
pattern  or  input  impedance.  This  means  that  these  elements  could 
be  loaded  with  PIN  diodes  at  a  number  of  points  and  the  center  fre 
queney  could  be  varied  over  a  wide  range  by  simply  biasing  different 
diodes  on  or  off. 

The  computed  voltage  distribution  of  the  loaded  mode  around  the 
edge  of  the  patch  for  a  mode  shorted  at  a  single  point  is  illustrated 
in  Fig.  5.  The  voltage  distribution  at  first  appears  to  be  severely 
perturbed  from  that  of  the  (0,  l)  mode  of  the  corresponding  unloaded 
antenna.  A  closer  look,  however,  will  show  that  voltage  distribution 
of  the  loaded  antenna  can  be  thought  of  as  that  of  the  unloaded 
clement  plus  an  oscillating  magnetic  current  which  has  little  effect 
on  the  radiation  pattern.  Thus,  much  the  character  of  what  might  be 
termed  the  “parent”  mode — the  unloaded  antenna  mode  from  which 
the  loaded  mode  was  adapted — is  carried  over  to  the  loaded  antenna 
mode.  The  predicted  results  are  consistent  wim  experimental  data 
as  shown  in  a  typical  case  illustrated  in  Fig  6. 

Polarization-diverse  antennas —  Microstrip  antennas  can  also  be 
designed  to  produce  a  field  along  the  broadside  direction  whose  po¬ 
larization  can  be  adjusted  by  turning  on  or  off  one  or  more  PIN 
diodes  (3),  or  by  adjusting  the  bias  on  one  or  more  varactor  diodes 
(12).  Two  nearly  degenerate  modes  can  be  excited  by  a  single  feed 
with  one  mode  producing  a  far-field  component  which  is  orthogonal 
to  that  of  the  other  mode.  By  aJtering  the  difference  in  the  resonant 
frequencies  of  the  two  modes,  one  can  adjust  the  relative  phasing  of 
the  modes.  Thus,  one  can  control  the  polarization  of  the  total  far 
field  produced  by  the  two  modes.  Elements  have  been  built  whose  po¬ 
larization  could  be  altered  continuously  between  left-  and  right-hand 
circular  polarization.  A  drawback  of  this  type  of  element  is  that  the 
band  over  which  the  axial  ratio  is  under  3dB  is  quite  narrow  (35/C? 
percent)  (12} 


Figure  5.  Magnetic  current  distribution  of  the  resonant  mode  of  a  frequency- 
agile  element  with  a  single  short  located  at  x *  5.72  cm  and  y  *  1 .52  cm.  The 
patch  dimensions  are  1 1  43  x  7.62  x  0. 1 52  cm  and  the  dielectric  constant  is 
2.43 
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o-o  Measured 
o  o  Computed 
5  MHz  increment 

Figure  6.  Input  impedance  of  a  frequency-agile  element  described  in  Fig.  5 
with  the  feed  point  is  at  (2.54,  1.27). 

Dual-band  elements — In  some  applications,  two  or  more  rather 
narrow,  discrete  frequency  band.»  of  operation  are  often  required  with 
band  separations  of  only  a  few  percent.  One  technique  that  has  been 
used  is  to  stack  elements  of  different  sizes  (14),  (15J.  One  experimental 
result  using  a  reactive  load  was  reported  in  (2).  The  present  theory 
shows  the  type  of  load  that  is  needed  and  how  to  control  the  band 
reparation. 

Figure  2(a)  illustrates  the  uses  of  a  load  reactance  which  has  a  pole 
lower  than  the  resonant  frequency  of  the  empty  cavity  resonant  par¬ 
ent  mode.  If  adjusted  properly,  two  intersections  between  the  XT 
curve  and  the  -(Aff  -f  A'l)  curve  occur  corresponding  to  two  daugh¬ 
ter  resonant  modes  of  the  loaded  cavity.  The  band  separation  can 
be  controlled  in  two  ways. 

1.  If  the  magnitude  of  the  slope  of  the  -(A'r  -1-  A'l)  curve  is  in¬ 
creased,  the  band  separation  is  reduced. 

2.  If  the  load  is  placed  closer  to  a  zero  of  the  parent  mode  function, 
then  the  Xr  curve  becomes  thinner  again  reducing  the  band 
separation. 

Of  these  two  methods,  the  latter  is  the  easiest  to  implement.  The 
resonant  load  can  be  simply  an  open-  or  a  short-circuited  transmis¬ 
sion  line,  which  can  be  etched  on  the  same  plane  as  the  patch  itself. 
In  this  case,  the  patch  may  need  to  be  slotted  and  the  line  recessed 
as  illustrated  in  Fig.  7  in  order  to  control  the  band  separation  and 
still  maintain  at  least  bilateral  symmetry  in  the  element  to  reduce 
cross  polarization. 

It  is  clear  from  Fig.  2(a)  that  there  is  a  continuous  set  of  possible 
pairs  of  dual  resonant  frequencies  obtainable  by  varying  the  resonant 
frequency  of  the  load.  However,  it  turns  out  that  there  is  really  only 
one  pair  of  resonant  frequencies  which  yields  resonant  modes  whose 
input  impedances  are  approximately  the  same.  For  all  other  pairs, 
the  maximum  input  resistance  of  one  is  larger  than  the  other.  In 
exlreiw  c;us«‘s,  one  inode  is  so  pooily  matched  that  it  is  useless. 
Figure  S  contains  the  measured  input  impedance  of  both  modes  of  a 
typical  dual-band  dement 

DC  mode  applications —  In  all  of  the  applications  discussed  above, 
the  parent  mode  used  almost  exclusively  was  the  (0,1)  mode  of  a 
rectangular  element.  This  mode  produces  a  pattern  with  a  beam 
maximum  in  a  direction  perpendicular  to  the  patch.  If  mounted  on 


Figure  7.  A  monolithic  duaJ-band  element  using  a  short-circuited  transmis¬ 
sion  line  recessed  into  the  patch. 


a  large  grounded  dielectric  slab,  the  patten.  Calls  off  to  zero  in  any 
direction  parallel  to  the  patch.  If  the  dielectric  is  truncated  at  the 
patch  edge,  the  gain  also  falls,  but  not  to  zero.  For  the  purposes 
of  the  present  discussion,  we  will  assume  that  the  dielectric  is  so 
truncated  Another  mode  of  a  nonloadcd  microstrip  antenna  is  the 
DC  mode — a  mode  resonant  at  zero  frequency  which  has  a  constant 
voltage  distribution  from  point  to  point  on  the  patch.  If  this  mode 
could  lx  excited  at  an  RF  frequency,  the  pattern  it  would  produce  is 
basically  the  same  as  that  of  an  electric  monopole  element  mounted 
perpendicular  to  the  ground  plane.  This  yields  a  pattern  which  has 
a  broadside  null  and  a  maximum  along  the  ground  plane.  The  first 
question  to  be  addressed  is  how  to  resonate  this  mode  at  a  useful 
frequency  instead  of  at  DC. 

To  answer  this  question,  we  consider  the  DC  mode  to  be  a  parent 
mode  and  plot  its  corresponding  X,  and  - A'r  curves  as  illustrated 
in  Fig.  9.  It  turns  out  that  if  Xl  =  0,  the  resonant  frequency  (oc¬ 
curring  where  the  dashed  curve  and  the  Xr  curve  intersect  in  Fig.  9) 
is  still  a  little  too  low  to  be  useful.  To  raise  the  frequency  even 
higher,  an  appropriate  amount  of  capacitance  can  be  added-  If  the 
load  capacitor  is  placed  in  the  center  of  the  square  patch,  then  the 
resonant  frequencies  of  the  (0,1)  and  (1,0)  modes  of  the  antenna  are 
unaffected  and  all  three  modes  can  be  made  to  resonate  at  the  same 
frequency.  An  alternative  to  using  a  capacitive  load  is  to  use  four 
short  circuits  as  illustrated  in  Fig.  10.  The  voltage  distribution  along 


Figure  8.  Measured  input  impedances  of  both  modes  of  a  typical  dual-band 
element  using  a  shorted  transmission  line  as  its  resonant  load. 
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Figure  9.  Plot  the  of  xfi  -Xf.  and  -( Xf+Xi ,)  curves  for  a  square  patch  ele¬ 
ment's  DC  mode.  The  X]^  corresponds  to  a  capacitive  load. 


one  edge  of  the  patch  is  also  illustrated — the  distributions  along  the 
remaining  edges  are  identical.  This  again  illustrates  the  principle 
that  the  daughter,  loaded  mode  has  the  features  of  its  parent  mode, 
although  it  is  not  identical.  The  magnetic  current  distribution  is 
basically  constant  around  the  edge  with  an  oscillating  perturbation 
having  little  significance  in  the  far  field.  At  least  three  applications 
of  this  DC  mode  can  be  imagined- 

End-fire  element.  By  combining  with  the  proper  magnitude  and 
phase,  the  magnetic  current  associated  with  the  DC  mode,  which  cir¬ 
culates  in  one  direction  around  the  patch,  and  the  magnetic  current 
of  a  (0,1)  mode,  which  is  directed  oppositely  on  opposite  radiating 
edges,  one  can  achieve  an  end-fire  pattern.  By  "end-fire”  we  mean 
a  pattern  with  a  beam  maximum  in  one  direction  parallel  to  the 
patch.  This  is  not  possible  using  a  single  mode  of  a  microstrip  an¬ 
tenna.  Moreover,  it  is  probably  not  possible  using  a  pair  of  naturally 
occurring  degenerate  modes  (such  as  the  (0,1)  and  (1,0)  modes  of  a 
square  patch)  since  such  inodes  typically  yield  differently  polarized 
fields — such  pairs  of  degenerate  modes  are  useful  for  creating  circular 
polarization,  however.  In  order  to  excite  both  the  loaded- DC  mode 
and  the  (0,1)  mode  independently ,  one  must  feed  each  mode  either 
at  a  point  where  the  other  mode  is  zero,  or  with  a  symmetry  opposite 
to  that  of  the  other  mode.  The  element  illustrated  in  Fig.  11  uses  a 
combination  of  these  two  methods.  The  dashed  circular  curve  in  the 


Magnetic  current 


figure  illustrates  a  line  along  which  the  electric  field  of  the  loaded-DC 
mode  is  zero.  In  principle,  one  should  be  able  to  feed  the  fO,  1)  mode 
at  that  point  without  coupling  to  the  DC  mode.  A  symmetrical  feed 
is  introduced  at  a  pair  of  points  at  which  the  field  of  the  J0,l)  mode 
is  zero  ih'.s  preventing  coupling  to  this  mode.  In  additioi.  the  (1,0) 
has  odd  symmetry  and  thus  is  also  not  excited  by  the  symmetrical 
feed.  The  only  mode  that  this  symmetrical  feed  can  couple  to  is  the 
loaded  DC  mode.  This  element  was  built  and  its  performance  was 
consistent  with  the  design  principles  outlined  here. 

Adaptive  element.  The  element  described  in  the  preceding  paragraph 
is  essentially  a  iwo-clcmcnt  array  within  a  single  patch  and,  as  such, 
lias  a  pattern  that  can  be  adapted  by  changing  the  excitation  of  the 
two  inodes.  Anolhci  type  of  adaptive  element  requiring  a  -ir.gle  futJ 
is  again  that  shown  in  Fig.  10,  but  with  the  short  circuits  imple¬ 
mented  as  PIN  diodes.  When  biased  on,  the  element  would  produce 
a  inonopole  pattern.  When  the)'  are  biased  off,  the  elemei*.  produces 
the  usual  (0, 1)  mode  pattern. 

Monopulse  element.  By  having  a  combination  of  four  fe~is  located 
at  the  center  of  each  edge  of  the  patch,  one  can  produce  a  sum  and 
difference  pattern  by  appropriate  even  or  odd  excitatioi  of  oppo* 


Figure  II  .  An  element  supporting  three  degenerate  modes:  tbe  loaded-DC 
mode,  the  (0,1)  mode,  and  the  (1,0)  mode 


site  pairs  of  feeds.  This  application  calls  for  the  capacitiv*ly  center- 
loaded  square  patch  supporting  degenerate  loaded-DC.  '0.1),  and 
(1,0)  modes.  If  a  pair  of  feeds  on  opposite  edges  are  excited  in 
phase,  then  only  the  DC  mode  is  excited  Exciting  the  ?air  out  of 
phase  would  excite  only  cither  the  (0.1)  mode  or  the  1 1 . 0)  mode 
depending  on  which  pair  of  feeds  were  used. 

Adaptive  impedance  matching —  Another  application  : i  to  change 
the  input  impedance  of  an  element  dynamically  in  response  to  a 
changing  environment  in  which  the  element  is  placed,  for  example  in 
a  phased  array.  In  this  application,  one  seeks  to  hold  the  resonant 
frequency  and  the  pattern  constant.  If  one  introduces  a  symmetri¬ 
cally  located  pair  of  short  circuits  that  lie  along  the  curve  illustrated 
in  Fig.  12.  the  maximum  input  resistance  varies  as  shown  in  Fig.  13. 


Patch 


l-igmc  10.  .Square  patch  element  with  four  short  circuits  protucing  a  loaded- 
DC  mode.  The  voltage  distributions  along  one  edge  and  alorg  a  cm  through 
a  centerline  of  the  element  are  also  shown 


16 


Figure  12.  The  locus  of  points  along  which  a  symmetrically  placed  pair  of 
shorts  will  result  in  a  constant  resonant  frequency  and  an  essentially  constant 
pattern. 
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i-component  of  short-circuit  location 


Figure  13.  Variation  of  the  maximum  input  resistance  as  the  load  pair  is 
moved  along  the  curves  illustrated  in  Fig.  12. 

Pattern  shaping —  There  is  not  a  lot  of  shaping  of  the  pattern  that 
one  can  do  with  a  microstrip  element  since  it  is  not  electrically  large 
and  the  magnetic  current  distribution  must  be  essentially  real-valued. 
However,  one  can,  io  principle,  obtain  a  more  nearly  isotropic  pattern 
in  the  £-plane  by  appropriately  loading  one  edge  of  the  element. 
(Again,  the  dielectric  must  be  truncated  at  the  edge  of  the  patch). 
The  concept  in  this  application  is  to  place  loads  at  one  edge  so  that 
the  net  magnetic  current  at  that  edge  is  zero  as  illustrated  in  Fig.  14. 
The  magnetic  current  at  the  opposite  edge  is  not  zero  resulting  in  a 
single  magnetic  dipole  for  this  pair  of  edges.  There  is  a  net  magnetic 
dipole  on  the  remaining  pair  of  edges,  but  they  do  not  contribute 
any  field  in  the  £-plane  and  their  maximum  field  strength  in  the  H- 
plane  has  been  computed  to  be  down  by  about  9dB  from  the  main 
beam.  This  yields  the  desired  isotropic  pattern.  Some  shaping  of 
this  pattern  is  possible  by  adjusting  the  net  magnetic  current  at  the 
loaded  edge. 

4  Conclusions 


The  basic  theory  of  loaded  microstrip  antennas  and  numerous  prac¬ 
tical  applications  have  been  discussed.  Proper  selection  and  location 


Figure  14.  Voltage  distribution  along  the  edge  of  a  rectangular  patch  loaded 
with  a  pair  of  capacitors.  The  net  magnetic  current  on  one  edge  is  zero,  but 
is  not  zero  on  the  other 


of  reactive  loads  can  be  used  to  control  the  resonant  frequency,  po¬ 
larization,  impedance,  and  radiation  pattern.  These  characteristics 
can  then  be  made  electronically  adaptable  to  various  systems  needs 
if  the  loads  are  PIN  or  varactor  diodes. 
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Dielectric  Resonator  Antenna 


agreement  is  seen  to  be  quite  good  considering  all  the  a'“  *  £r  -  13.2,  tan6  =  0.0064,  probe  length1 


;i-.euretical  and  experimental  impedance,  a  >  1.27  Fig.  0  Experimental  (tap)  and  theoretical  using  measured 

a/d  --  1.0,  Er  =  13.2,  tan<5  =  0.0004,  probe  ibottom)  impedance,  a  •  1.27  cm,  a/d  ■  0.5,  Er  - 
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Abstract— An  asymptotic  extraction  process  is  developed  which  leads 
to  an  efficient  method  foi  the  evaluation  of  the  eleclric  field  from  a 
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Hertzian  dipole  in  a  layered  geometry.  This  technique  allows  for  a  rapidly 
converging  expression  for  Ihe  electric  field  which  remains  well  behaved  as 
the  source  and  observation  heights  i' ,  z  coincide.  A  useful  application  of 
this  method  is  in  the  calculation  of  mutual  impedance  between  printed 
dipoles. 

I.  Introduction 

For  the  rigorous  analysis  of  microstrip  antenna  problems,  a  key 
step  is  usually  the  evaluation  of  Sommerfeld-type  integrals  to 
compute  the  reaction  between  current  functions.  This  is  often  the 
most  time-consuming  pan  of  the  total  calculation,  due  to  slow 
convergence.  Funhermore.  because  the  Sommerfeid  integral  for  the 
electric  field  from  a  Hertzian  dipole  becomes  nonconvergcnt  as  the 
source  and  observation  heights  z‘  and  z  coincide,  certain  techniques 
must  be  used  for  calculations  in  the  real-space  domain  (I)  An 
alternative  approach  is  the  plane-wave  spectrum  method  (2),  (3).  This 
method  eliminates  the  nonconvergent  behavior  when  z  «  z' ■  Self- 
term  problems  arc  also  eliminated,  making  the  method  very  attractive 
for  the  numerical  analysis  of  antenna  problems.  One  disadvantage  of 
this  method,  however,  is  that  the  convergence  of  the  resulting  integral 
expression  depends  on  the  Fourier  transforms  of  the  current  functions 
when  z  =*  z‘ ■  The  integrand  thus  oscillates  more  rapidly  as  the 
separation  between  currents  increases,  and  the  convergence  becomes 
worse  as  the  size  of  the  current  functions  decreases.  Hence  it  is 
difficult  to  compute  the  reaction  between  small  basis  functions  of 
current  which  arc  widely  separated.  One  technique  for  overcoming 
this  difficulty  has  been  introduced  by  Pozar  [4],  In  this  technique  a 
homogcncous-spacc  term  is  extracted  from  the  integrand  in  the 
spectral  domain,  resulting  in  a  more  rapid  convergence.  This 
technique  is  well-suited  for  an  efficient  analysis  of  patch  and  dipole 
antennas.  The  convergence  of  the  resulting  integrand,  is  still  slow  for 
small  current  functions  when  z  =  z',  however. 

In  this  communication  an  asymptotic  extraction  technique  for 
directly  evaluating  the  electric  field  of  a  Hertzian  dipole  is  intro¬ 
duced.  in  which  convergence  is  obtained  for  all  values  of  z  and  z' . 
Because  the  extraction  technique  is  performed  in  the  real-space 
domain,  the  method  efficiently  computes  the  reaction  between  basis 
functions  of  arbitrarily  small  size.  The  resulting  expression  is  also 
rapidly  convergent,  allowing  for  an  efficient  computation  even  for 
large  basis  function  separations. 

II.  Asymptotic  extraction  technique 

The  geometry  of  a  Hertzian  dipole  embedded  in  a  layered  structure 
over  a  ground  plane  is  shown  in  Fig.  1 .  Although  the  general  method 
to  be  discussed  is  applicable  to  an  arbitrary  number  of  layers,  an 
explicit  solution  will  only  be  given  for  the  case  of  a  dipole  within  a 
single  layer  of  thickness  B,  as  shown  in  the  figure.  Likewise,  the 
method  is  applicable  to  layers  of  arbitrary  loss,  although  results  will 
only  be  shown  for  the  lossless  case.  The  Hertzian  dipole  is  assumed 
oriented  in  the  x  direction,  at  a  distance  z  from  the  ground  plane. 
The  relative  permittivity  and  permeability  of  the  layer  arc  denoted  as 
er,  firy  with  n  -  Ve,pf  the  index  of  refraction. 

The  starting  point  of  the  analysis  is  the  Sommerfeid  solution  for  the 
Hertzian  potential  components  IT,,  n.  at  x,  y,  z  due  to  (he  dipole  at 
x' ,  y' ,  z' .  This  solution  may  be  written  in  the  form  (suppressing 

e*nis\ 

n,=  |" /(X)-/o(Xr)  rfX  (I) 

Jo 

n.=  [  g(\)Jt(\r)  COS  <t>  d\  (2) 

Jo 

in  which 

r=  ((*-*' )2  +  (v- ■y')’1)''1,  cos  <t>  = - — 


z 


Fig.  I .  Hertzian  dipole  embedded  in  a  single  layer. 


and /(X),  g(X)  are  complicated  functions  of  the  integration  variable  X 
which  also  depend  on  z  and  z\  although  this  dependence  is  not 
explicitly  shown  here.  For  the  case  of  a  single  layer  with  z,  z'  S  5, 
the  functions  /  and  g  arc  given  by 

Z  <  z'  <  B 


/(X)= 


2  X 
D,(X)  u, 


—  u ,  cosh  U| (/?-*') 

A. 


Z  <  B 


+  u  sinh  uj  (B  -  z ') 


sinh  u ,z  (3) 


/(X)  = 


2  X 

A-(X)  u, 


—  u,  cosh  Ui(B~z) 

Vr 


and  in  either  region, 


+  u  sinh  u, (fl- z) 


sinh  UiJ'  (4) 


£(X)  = 


0 -n1) 


Df(\)Dm(\)  nr 


X2  sinh  u,z'  cosh  utz  (5) 


where 


Dr(\)  =  u  sinh  u,Z?  +  —  uy  cosh  u,B  (6) 

Mz 

un1  cosh  U|5  +  U|  sinh  uyB  (7) 

Mr 


with 


w^CX2-*-2)1'2,  I/,  =  (X2  -  *2) l/2,  and  k\=k0n. 

In  (3)-(5)  a  normalization  factor  -jwnoHr/4*k2  has  been  sup¬ 
pressed.  An  important  feature  of  the  coefficient  functions  /(X)  and 
g(\)  is  that  they  arc  both  0(1)  as  X  -*  oo  when  z  =  Z ' ,  and  they  both 
decay  exponentially  as  X  -*  oo  otherwise.  This  general  property  holds 
true  for  any  arbitrary  number  of  layers.  These  functions  also  exhibit 
surface-wave  poles  at  the  roots  of  the  functions  D,(X)  and  D„(\)  in 
the  region  kQ  <  X  <  kon,  as  well  as  branch  points  at  X  =  ±k0.  The 
integration  in  (I)  and  (2)  is  along  a  Sommerfeid  contour  which  goes 
above  these  poles,  and  which  maintains  the  proper  branch  interpreta¬ 
tion  of  u  (1). 

To  illustrate  the  extraction  technique  the  Ex  component  of  the 
electric  field  is  considered,  given  by 


Er 


k2  nr+ 


d2nt 

dx2 


32n. 

dxdz 


(8) 


where  k  is  the  wavenumber  at  x,  y,  z. 
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To  calculate  Ex,  (2)  is  first  differentiated  as 

f"A(X)y,(Xr)  cos  4  X  d\ 
oz  Jo 

(9) 

where 

I  d* 

X  dz 

(10) 

The  function  A(X)  is  also  0(1)  as  X  -*  oo  when  z  -  z' ,  and 
exponentially  decaying  otherwise.  Hence  (9)  is  nonconvergent  when 
z  -  z' .  To  overcome  this  problem,  an  asymptotic  extraction  term 
may  be  formulated,  which  asymptotically  approximates  the  behavior 
of  the  integrands  in  (1)  and  (9).  Such  a  form  may  be  obtained  by 
approximating  the  hyperbolic  sinh  and  cosh  functions  appearing  in 
the  f  and  h  functions  with  exponential  functions  of  X.  In  particular, 
the  following  relations  are  used: 

sinh  (if;#)  —  j  6ae  +  Xtf 

cosh  (w,a)  — ^  (2-60)e*  *‘a 

as  X  -*  oo,  with 

fu  o>0 
^"(O,  0  =  0. 

Using  these,  the  coefficient  functions  are  seen  to  become  asymptotic 
to  the  following  forms  as  X  -»  oo; 

fW~t  Te'Mn 

(ID 

«  A.(X), 

(12) 

and 


The  number  N  depends  on  the  number  of  layers  in  the  problem,  with 
N  =  I  for  the  single  layer  case.  The  constants  a„,  b„  depend  only  on 
the  layer  thicknesses  and  material  parameters,  while  the  constants 
depend  on  z  and  z'  as  well  as  the  layer  thicknesses.  For  example, 
for  the  single  layer  case,  the  following  results  arc  obtained  from  (3)- 
(7)  for  the  case  z  £  z ' : 


Hr+l 


[2-6g-z>  +  Hr&B-t'  ) 


?i  =*z‘-z 


and 


2(1 -rr3) 


(p,+  l)(«r+l) 

Cl  =  2  B-z-z’. 


In  (13).  (15)  a  factor  -30 juJn1  has  been  suppressed. 

The  asymptotic  terms /„(X)  and  A,(X)  are  now  subtracted  and  added 
from  the  integrands  in  (I)  and  (9).  respectively.  The  resulting 
expressions  front  the  added  terms  may  Ik  evaluated  in  closed  foun. 
using  the  relations 


d\  = 


1 


T-M^r)=  cos  <#>). 

dx 


(18) 


Substituting  Ihe  results  irdo  (S)  and  differentiating  then  results  in 

£>  f'  £(X)  d\  +  k>  2  T^+  2  T  T»  (19) 

J0  Kq  Ag  i  *0 


/r(X)  =  [/(X)-/<r(\)]f,(X)  +  [/I(A)-Afl(X)]£2(X)  (20) 


with 


£,(X)  =  (Ar3-X3  cos!  d)70(Xr)- X!(l  -2  cos1  4)  (21) 

£j(X)  =  X3  cos!  4To(Xr)  +  X 2 (1  - 2  cos3  4)^4^  (22) 


Xr 


and 


r'”“(r3  +  f3)''3 


r2„  = 


r3  +  f„- Xx-x’)2 
(r3  +  {3) s,! 

rJ  +  tl-Xx-x')1 


(c!  +  f?)s 


(23) 

(24) 

(25) 


The  functions /(X)  and  A(X)  approach  their  corresponding  asymp¬ 
totic  forms  fairly  quickly  in  general,  as  shown  in  Fig.  2  for  the  case  of 
a  single  teflon  layer  with  z  =  z'  -  B  (the  asymptotic  functions  are 
constants  in  this  case,  since  and  fi  are  both  zero).  The  terms /(X) 
-  /#(X)  and  A(X)  -  /r,(X)  arc  0(1  fXs)  as  X  -*  ®.  so  the  integration  in 
(19)  converges  for  all  z  and  z" 

Because  of  the  surface-wave  poles  the  integration  in  (19)  may  be 
broken  into  iwo  regions:  (0,  A)  and  ( A ,  oo),  where  A  is  larger  than 
all  surface-wave  poles.  A  pole  extraction  technique  may  then  be  used 
to  evaluate  the  integration  over  the  first  region  [I],  The  integration 
over  the  second  region  may  be  improved  by  noting  that 


and 


(13) 

(H) 

(15) 

(16) 


lim  X![/(X)-/„(X)]  =  c, 


lim  X3[A(X)-A„(X)]  =  c2 


(26) 


(27) 


(r3  +  03)'' 


(17) 


where  c( ,  c2  are  constants  which  are  nonzero  when  z  =  z'. 

Assuming  that  the  limiting  results  of  (26),  (27)  are  valid  when  X  > 
AT  for  some  sufficiently  large  value  AT,  the  integration  from  AT  tc  os 
may  be  evaluated  analytically  as 

f"  £(X)  rfX  =  ^^3  [d  -cos3  4)7(ATr)  +  (1  - 2  cos3  4)7, (ATr)] 
r 

-  C\n2r[f(Mr)  +  J0(Mr)/ Mr)  (28) 

where 

l(x)  s  f  Ml)  dl  =  -y,(AT)(l  -X'1) 

+  J-^1  (I  -3j--!)  +  0(x--m)  (29) 

X 

a  s  x  —  oo . 
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X/kQ 

Fig.  2.  A  plot  of  the  coefficient  functions /(X),  Jt(X),  and  the  corresponding 
asymptotic  functions  /,(X),  Jt4(X)  (in  all  functions  a  factor  j/k0  has  been 
suppressed). 

Choosing  M/k0  -  25  will  usually  insure  that  the  limiting 
approximations  of  (26),  (27)  as  well  as  that  of  (29)  are  valid.  Hence 
the  integration  in  X  need  not  usually  extend  past  this  value  of  M  in 
order  to  obtain  accurate  results  for  all  values  of  z  and  z' . 

As  an  illustration  of  results  obtained  with  this  method,  the  real  and 
imaginary  pans  of  the  £>fieid  are  plotted  versus  x-axi$  separation  in 
Fig.  3  (Ax  =  x  -  x')  for  the  case  <t>  =  0  with  a  teflon  layer.  The  Ex- 
field  decays  slowly  with  r  in  this  case  due  to  the  endfire  nature  of  the 
dominant  TMi  surface-wave  mode. 

One  practical  use  of  the  £-field  computation  by  this  method  is  in 
the  calculation  of  mutual  impedance  between  dipoles  [8].  When  using 
the  method  of  moments  with  subdomain  basis  functions,  the  reaction 
between  widely  separated  basis  functions  on  different  dipoles  may  be 
computed  efficiently  by  approximating  the  currents  as  point  sources, 
which  reduces  the  computation  to  a  single  £*-field  evaluation.  The 
plane-wave  spectrum  method  may  be  used  for  reactions  along  the 
same  dipole,  thus  using  both  methods  in  a  complimentary  way. 

III.  Conclusion 

An  asymptotic  extraction  technique  for  the  efficient  real-space 
evaluation  of  the  electric  field  from  a  Hertzian  dipole  in  a  layered 
media  has  been  presented.  The  method  involves  the  extraction  of  z- 
dependent  terms  from  the  integrands  in  the  Sommerfeld  expressions 
for  the  Hertzian  potential  components,  resulting  in  an  expression 
which  converges  rapidly  for  all  values  of  z  and  z' .  This  method  is 
well  suited  for  finding  the  reaction  between  basis  functions  of  current 
which  are  small  and  widely  separated,  for  which  a  point  dipole 
approximation  may  be  used. 
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Analysis  of  Planar  Strip  Geometries  in  a 
Substrate-Superstrate  Configuration 
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Abstract — A  momenl  method  procedure  is  used  to  analyze  the  behavior 
of  several  different  configurations  consisting  of  planar  strips  in  a 
substrate-superstrate  geometry.  These  include  the  microstrip  transmission 
line,  center-fed  dipole,  the  mutual  impedance  between  two  dipoles,  and 
the  transmission-line  coupled  dipole.  In  each  case  some  of  the  basic 
superstate  effects  are  discussed. 

I.  Introduction 

ECENTLY,  some  of  the  basic  effects  of  a  superstate 
(cover)  layer  on  printed  circuit  antenna  performance  were 
investigated,  based  on  the  analysis  of  a  Hertzian  (infinitesimal) 
dipole.  By  properly  choosing  the  layer  thicknesses  and 
material  constants,  it  was  observed  that  several  different 
phenomena  could  occur,  each  of  which  could  be  used  in  a 
different  way  to  improve  the  radiation  characteristics  of  the 
dipole  [l]-[4].  These  results  extend  directly  to  apply  to  the 
radiation  characteristics  of  full-size  printed  antennas.  How¬ 
ever,  for  the  full-size  antenna,  input  impedance  and  mutual 
coupling  problems  are  also  of  concern,  and  these  require 
calculation  of  the  near-field  together  with  a  suitable  numerical 
technique.  The  purpose  of  this  investigation  is  the  numerical 
analysis  of  four  representative  classes  of  full-size  geometries 
in  a  substrate-superstrate  configuration,  each  composed  of 
planar  strips.  These  four  geometries,  shown  in  Fig.  1,  include 
the  microstrip  transmission  line,  the  center-fed  dipole,  the 
mutual  impedance  between  two  dipoles,  and  the  transmission¬ 
line  coupled  dipole.  In  each  case  the  conductive  strips  of  width 
w  are  embedded  within  a  two-layered  structure  consisting  of  a 
grounded  substrate  of  thickness  B  with  relative  permittivity 
and  permeability  U ;  with  a  superstrate  of  thickness  f  on  top, 
having  parameters  e;,  p2.  The  strip  axes  are  taken  to  be  in  the  x 
direction  in  all  cases.  An  illustration  of  the  layer  geometry  is 
shown  in  Fig.  2  for  the  case  of  the  microstrip  line.  A 
convenient  analysis  for  each  of  these  cases  can  be  provided  by 
using  a  plane-wave  spectrum  formulation  for  the  reaction 
between  two  arbitrary  x-directed  current  sources,  together 
with  a  Galerkin  moment  method  formulation. 

After  the  formulation  of  each  problem  is  given,  results  arc 
shown  and  the  effects  of  a  superstrate  layer  are  discussed.  For 
the  microstrip  transmission  line,  the  effect  of  the  superstrate 
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Microstrip  Line  Center-Fed  Dipole 


Mutual  Impedance  Between  Iwo  Dipoles 


TOP  VIEW 

Coupling  Between  Transmission  Line  and  Dipole 
Fig.  1  Planar-strip  geometries  in  a  substrate-superstrate  configuration. 


Fig.  2.  Substrate-superstrate  geometry. 


on  the  line  dispersion  is  demonstrated.  For  the  center-fed 
dipole,  the  superstrate  effect  on  the  resonant  length  and  input 
resistance  is  seen,  and  simple  formulas  are  given  to  approxi¬ 
mately  calculate  these  quantities.  The  effect  of  a  superstrate 
layer  on  the  mutual  impedance  between  two  dipoles  is  then 
discussed,  and  it  is  demonstrated  that  a  superstrate  may  be 
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used  lo  signiflcamly  reduce  ihe  mutual  impedance  for  dipoles 
in  an  endfirc  configuration.  Finally,  the  effect  of  a  superstrate 
layer  on  the  coupling  between  a  dipole  and  transmission  line  is 
examined,  and  the  capacitive  nature  of  the  coupling  is 
demonstrated. 

!!.  Formulation  for  Current  Reaction 
A  foundation  for  the  numerical  analysis  of  full-size  planar 
strip  structures  is  the  calculation  of  the  reaction  between  two 
x-directed  currents.  Denoting  J ,(x,  y)  and  J2(x,  y)  as  two  such 
currents  located  at  z  =  z,  and  z  =  z2l  respectively,  the 
reaction  is  defined  as 

ri2  =  (Ji,  Ji)=  J  j  E,,(x,  y)J2(x,  y)  dx  dy  (I) 

where  Ex,  is  the  £, -field  from  current  J,(x.  y)  and  the 
integration  is  performed  in  the  plane  z2.  A  starting  point  for 
this  calculation  is  the  Sommerfeld  solution  for  the  magnetic 
Hertzian  vector  potential  at  (x,  y,  z)  due  to  an  electric  Hertzian 
dipole  source  at  (x',  y' ,  z')  where  z'  =  Zi.  This  can  be 
written  in  cylindrical  coordinates  as 


Both  (2)  and  (5)  arc  then  in  the  form  of  cylindrical  (Fourier- 
Bessel)  transforms,  and  a  conversion  to  rectangular  (Fourier) 
form  may  be  accomplished.  Substituting  (4)  into  (I)  then 
results  in  the  plane-wave  spectrum  form 

rn=Tx  n  n*  i 

MK.  X2) — Xx,  —  Xy)  d\x  d\y  (8) 

where 

r.  Xj,)=J  j  y(x' ,  y ')e-/IVr'  ♦yr'J  dx'  dy' 

and 

xWx’  +  x;. 

Defining  new  origins  (x, ,  and  (x2,  y2)  for  currents  J,,J2, 
and  assuming  the  currents  arc  even  functions  in  x  and  y  about 
the  new  origins.  (8)  can  be  written  as 


>1=  r/(X)70(Xr)  d\ 

•Ml 

(2) 

n;  =  cos  0  (  g(X)y,(Xr)  d\ 

(3) 

where 

r=[(x-x')!  +  (y-y')JJl/! 

and  /(X),  8(X)  are  complicated  functions  of  X  which  also 
depend  on  z,  z',  although  this  dependence  is  not  explicitly 
shown  here  (the  time  dependence  e*'"'  is  also  being  sup¬ 
pressed  throughout  the  analysis).  These  functions  are  given 
explicitly  in  (2).  The  £r-field  in  the  plane  z2  is  then  found  from 

£.,i(*,  y,  Zi)=  J  _  |  -/|(x',  y ') 


2  rT/1  r 

'•12  =  -  [(A-J-^)/(X)  +  X^(X))/,(X„  X,) 

T  J0  Jr 

•  /2(X„  Xj.)  •  cos  (X,Ax)  cos  (X^Ay)  d\  di  (9) 

where  Ax  =  x,  -  x2,  Ay  =  y,  -  y2,  and  4>  =  tan'1  ). 

The  contour  C  in  the  Sommerfeld  integration  is  from  0  to 
o°,  above  the  surface-wave  poles.  A  root  extraction  technique 
is  used  to  remove  the  singular  behavior  of  the  integrand  [5], 
(6),  and  other  numerical  techniques  are  used  to  improve 
computational  efficiency.  A  matrix  storage  of  the  complicated 
part  of  the  integrand  in  brackets  is  one  such  technique,  and  the 
use  of  a  Filon-type  integration  to  account  for  the  oscillatory 
nature  of  the  cosine  terms  is  another.  In  the  application  of 
Galerkin’s  method  to  the  antenna  geometries,  a  convenient 
choice  of  basis  functions  is 


3!n,.  3  TLX 


AJn,  +  — f + 


3xJ  3x3z/;»:2 


dx'  dy'  (4) 


where 


b(x,  y)  =  ?Wl(f) 


(10) 


where  the  integration  is  performed  in  the  plane  z,,  and  k  is  the 
wavenumber  at  z  =  z2,  given  as  A  =  Ayi*  when  z  is  in  the  ith 
layer,  with  n,  being  the  index  of  refraction  of  the  rth  layer. 
Because  of  the  differentiation  appearing  in  (4),  it  is  convenient 
to  write 

3n.  r« 

— -  =  cos  <t  /i(X)y,(Xr)XrfX  (5) 

3z  Jo 

where 

I  3g 

h(X)  =  --2.  (6) 

X  3z 

Equation  (5)  may  be  put  into  a  similar  form  as  (2)  by  using 
3 

—  y0(Xr)=-X  cos  07((Xr).  (7) 

dx 


sin 

?(*}  = — 

A  ,{d- |x|) 
sin  k,d 

(ii) 

4  , 

n(>)  =  — 

5  w 

'.  +  | 
V 

2 

-y, 

w 

O' 

W 

\y\*j 

(12) 

[7),  [8J.  An  advantage  of  using  (12)  over  the  Maxwell  current 
is  faster  convergence  of  the  transform  [8],  The  planar  strips  of 
width  w  are  then  divided  into  N  subsections  of  length  d  in  the  x 
direction,  and  the  representation 

y)  =  rj(y)  £  7„£(x-  nd  +  Z./2)  (13) 

n  -  1 

is  used  for  each  strip,  with  L  the  strip  length  and  x  measured 
from  the  strip  center. 

One  of  the  difficulties  in  using  (9)  is  that  the  cosine  terms 
oscillate  faster  for  increasing  separation  A*,  A y  between  basis 
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functions.  Furthermore,  convergence  is  provided  by  [he 
Fourier  transforms  J\,  /2  when  z,  =  z2  since /(X)  and  A( X)  are 
0(1)  as  X  -»  oo,  so  the  integration  in  X  must  extend  to  further 
values  as  the  size  of  the  basis  functions  decrease.  This  makes 
the  necessary  computation  for  mutual  impedance  very  time 
consuming  when  computing  the  reaction  between  basis  func¬ 
tions  which  are  small  and  widely  separated.  In  this  case  it  is 
convenient  to  formulate  the  reaction  in  the  real-space  domain. 
If  the  basis  functions  are  approximated  as  point  dipoles  of 
current  this  is  equivalent  to  simply  evaluating  the  Ex- field  of  a 
Hertzian  dipole,  and  this  may  be  formulated  in  an  efficient 
manner  [9]. 


III.  Microstrip  Transmission  Line 
Using  (8),  an  analysis  of  the  propagation  on  a  microstrip 
line  is  straightforward.  A  current  expansion  function  of  the 
form 


J i(->f,  y)  =  i){y)e-‘Sx  (14) 

is  assumed  on  the  line,  and  a  testing  function 

MX.  y)  =  t>(x)Hy)  (15) 

is  used,  which  forces  the  £,-fleld  at  the  strip  center  to  be  zero. 
This  is  equivalent  in  principle  to  the  original  formulation  of 
Denlinger  [8),  although  the  approach  is  based  here  on  the 
three-dimensional  Green's  function.  Using  the  integral  repre¬ 
sentation  of  the  5-function,  the  result 


f"|7(X,)K*J-i3V(X)  +  i3J/i(X)]  d\-0  (16) 

J0  A 


is  obtained,  where 


X  =  V/32  +  X} 

and  v(\y)  is  the  one-dimensional  transform  of  rj(y).  As  pointed 
out  in  [1 1)  the  value  of  0  must  be  greater  than  all  surface-wave 
poles,  since  a  residue  contribution  would  result  in  a  complex 
value  for  the  integration,  and  a  real-valued  solution  for 
would  not  be  possible. 

Equation  (16)  allows  for  a  convenient  solution  for  the 
propagation  on  an  infinitely  thin  strip.  As  w  -»  0,  -*  1 

and  does  not  provide  any  convergence,  so  the  term  in  brackets 
must  vanish  as  X  -»  oo,  resulting  in 


02~k2 


/(°°) 

/(oo)-h(co) 


(17) 


as  w  -+  0.  From  this  formula  the  classical  results  for  an 
infinitely  thin  wire  may  be  obtained  (12J,  [13],  where  the 
effective  permittivity  and  permeability  are 

<*  =  -(«,  +  <;).  /*,=  j^2  (Ff'  +  F/'lj  (18) 

when  the  strip  is  at  the  interface  between  the  / th  and y'th  layers, 
and  f e  =  e,,  iLe  =  /x,  if  the  line  is  embedded  within  the  /th 
layer.  Equation  (18)  may  be  used  as  a  basis  for  choosing  kr  in 
(II),  with  k,  =  k0sf^r  [6], 


Results  for  a  line  at  the  interface  (z  =  B)  are  shown  in  Fig. 
3(a)  for  the  case  of  a  GaAs  superstrate  (e2  =  12.5)  over  a 
teflon  substrate  (e i  =  2. 1).  As  the  frequency  increases,  r,  -*  e. 
since  e2  >  f|.  The  dispersion  is  less  for  the  smaller  line  width 
at  lower  frequencies,  but  regardless  of  line  width  the  solution 
is  above  the  dominant  TM,  surface-wave  mode  here.  As  the 
superstrate  thickness  increases  the  amount  of  dispersion  does 
as  well,  seen  from  Fig.  3(b).  This  is  due  to  the  fact  that  <}  is 
significantly  larger  than  If  a  superstrate  with  e2  =  e,  is  used 
instead,  the  dispersion  would  be  less  than  without  a  super¬ 
strate. 


IV.  Center-Fed  Dipole 

For  the  three  antenna  geometries  shown  in  Fig.  1 ,  a  moment 
solution  using  Galerkin's  method  is  used.  For  the  center-fed 
dipole  the  current  distribution  (13)  is  assumed,  and  a  5-gap 
excitation  is  taken  at  jr  =  0.  The  Z,„„  coefficients  of  the 
Galerkin  matrix  then  represent  the  reaction  between  basis 
functions  centered  at  jr  =  md  -  £/2  and  x  =  nd  -  L/2 
where  L  is  the  length  of  the  dipole.  The  Galerkin  set  of  matrix 
equations 


[Z, „„][/„]  =  [£,,]  (19) 

is  then  formulated,  where  the  excitation  column  vector  £  has 
components  E„  =  6„p  for  N  even,  where  p  -  N/2.  A 
variational  expression  for  the  input  impedance  is  [14] 

=  <20) 

where  /(0)  is  the  current  at  x  =  0.  Using  (19),  this  formula 
immediately  reduces  to  the  usual  result  Zin  =  1/7(0).  Once 
input  impedance  is  found  in  this  way,  the  resonant  length  may 
easily  be  determined.  If  the  resonant  length  Lr  is  small 
compared  to  a  wavelength,  an  approximate  formula  for  the 
input  resistance  is  [15) 


£,=  120 


2 


Pt 

I5tVX2 


(0) 


(21) 


where  pT  is  the  total  (radiated  +  surface  wave)  power  (in 
watts)  delivered  by  a  Hertzian  dipole  of  strength  1  A-m  at  the 
same  z  location. 

In  Fig.  4(a)  the  exact  resonant  length  of  a  dipole  at  the 
interface  is  shown,  and  compared  with  the  approximate  value 


£,/Xo=l/V2(f,  +  e2)  (22) 

as  predicted  by  (18)  when  fringing  fields  are  neglected  in  the 
infinitely  narrow  nonmagnetic  case.  When  the  thickness  t 
becomes  comparable  to  w  the  resonant  length  changes  rapidly, 
becoming  asymptotic  to  a  limiting  half-space  value  close  to  the 
value  predicted  by  (22).  In  Fig.  4(b)  the  resonant  input 
resistance  is  shown  for  the  same  case  of  Fig.  4(a),  and 
compared  with  (21).  From  (21)  it  is  seen  that  a  superstrate 
affects  resonant  input  resistance  in  two  ways;  by  shortening 
the  resonant  length,  and  by  changing  Pr.  The  latter  effect  is 
discussed  in  (1). 
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V.  Mutual  Impedance  Between  Two  Dipoles 
The  mutual  impedance  between  two  center-fed  dipoles  has 
been  analyzed  in  two  ways:  by  using  the  method  of  moments, 
and  by  using  the  electromotive  force  (emf)  method.  In  the 
moment  method  approach  basis  functions  are  assumed  on  each 
dipole,  with  dipole  1  excited  by  a  delta-gap  source  and  dipole  2 
short  circuited.  A  partitioned  Gaterkjn  matrix  set  of  equations 
is  then  solved  to  find  the  dipole  currents.  Using  [Z]  =  [  K) ' 1 
where  [  Y\  is  the  admittance  matrix  results  in  the  formula 


-  1 


/a(0)//,(0) 


(23) 


Other  formulas  for  Z[2  may  be  derived,  but  an  advantage  of 
this  form  is  the  avoidance  of  numerical  difficulties  when 
|/2(0)|  <  |/,(0)|,  which  occurs  for  larger  dipole  separations. 

In  the  emf  method  the  mutual  impedance  is  given  by 

Z^lMW)\\E*'{X'y)JAX'y)dXdy  (24) 

h 

where  Ex\  is  the  field  from  the  current  Jt  on  dipole  1  when 
excited  in  the  absence  of  dipole  2,  and  J2  is  the  current  on 
dipole  2  (region  denoted  S2)  with  dipole  l  open-circuited 
Currents  J ,  and  J2  are  usually  assumed,  making  the  method 
approximate.  In  the  case  of  infinitely  narrow  dipoles  the 
currents  become  sinusoidal  with  an  effective  k,  as  discussed  in 
Section  III  however,  allowing  for  an  exact  calculation.  Since 
no  self-term  problems  are  encountered  for  separated  dipoles, 
an  efficient  calculation  of  Z,2  may  be  achieved  using  integra¬ 
tion  in  real  space.  For  narrow  dipoles  this  approach  is 
computationally  more  efficient,  and  has  the  advantage  of  being 
numerically  stable  as  w  0  since  there  are  no  self-term 
problems. 

For  simplicity,  results  will  be  confined  to  the  special  case  of 
both  dipoles  at  z  =  Zo-  First,  a  comparison  of  the  moment 
solution  with  the  emf  method  is  shown  in  Fig.  5,  for  dipoles  in 
an  endfire  configuration  on  a  single  teflon  layer  (zo  =  B).  The 
real  and  imaginary  parts  of  the  mutual  impedance,  f?>2  and 
Xi2,  are  plotted  for  several  different  strip  widths.  The  moment 
method  solution  is  used  for  the  nonzero  strip  widths,  and  the 
emf  method  is  used  for  the  zero  strip  width.  As  can  be  seen 
from  the  figure,  the  moment  solution  converges  toward  the 
emf  solution  for  narrower  strip  widths,  supporting  the  validity 
of  the  assumed  current  using  (18).  Because  all  of  the  basic 
superstrate  effects  on  mutual  impedance  can  be  seen  from  the 
infinitely  narrow  dipole  case,  the  remainder  of  the  results 
shown  in  this  section  will  be  for  this  case. 

The  basic  properties  of  mutual  impedance  for  dipoles  in 
broadside  or  endfire  configuration  can  be  seen  from  Fig.  6,  for 
dipoles  on  a  single  layer.  The  mutual  impedance  decays  fairly 
rapidly  with  separation  in  the  broadside  case  (Fig.  6(a)),  but 
much  slower  in  the  endfire  case  (Fig.  6(b))  This  is  due  to  the 
facl  that  the  far-ficld  surface  wave  pattern  of  a  transverse 
magnetic  (TM)-mode  wave  is  strongest  in  the  endfire  direc¬ 
tion,  having  a  null  in  the  broadside  direction  [16j.  The 


AX  /A, 

Fig.  5.  Mutual  impedance  versus  separation  showing  a  comparison  of  the 
moment  method  solution  for  two  different  strip  widths  and  the  emf  solution 
for  zero  strip  width. 

opposite  is  true  for  a  transverse  electric  (TH)-mode  wave,  but 
for  the  layer  thickness  used  here  only  the  dominant  TM  n  ode 
exists.  This  strong  endfire  coupling  is  even  more  pronou  iced 
for  dipoles  on  a  magnetic  layer,  since  the  TM,  mode  is  nuch 
more  strongly  excited.  Such  a  case  is  shown  in  Fig.  7,  where 
Ut  =  1 0.0.  This  strong  endfire  coupling  can  be  reduced  by 
using  a  superstrate  layer  to  reduce  or  eliminate  the  sutface 
waves.  An  example  of  this  is  shown  in  Fig.  8,  in  which  a 
magnetic  layer  with  p2  =  10.0  is  used  as  a  superstrate  over  a 
teflon  substrate,  with  Zo  =  B.  The  superstrate  thickness  is 
chosen  from 


ftjlc  _  ”2 

\>  2n'ln\  — 

which  results  in  an  elimination  of  all  surface  waves  in  this 
case.  As  discussed  in  (1),  the  elimination  of  surface  waves  is 
often  difficult  to  achieve  when  using  a  nonmagnetic  super¬ 
strate.  because  of  the  thin  substrate  limitation.  However 
another  way  to  reduce  -urface  waves  is  by  embedding  the 
dipoles  within  the  substrate  layer  and  using  a  superstrate  with  a 
higher  index  of  refraction.  This  reduces  surface  wave  excita¬ 
tion  due  to  the  exponentially  decaying  nature  of  the  surface- 
wave  fields  inside  the  substrate.  An  example  showing  mutual 
impedance  for  this  configuration  is  presented  in  Fig.  9,  where 
Zp  =  B/2-  Once  the  separation  becomes  sufficiently  large  in 
this  case,  the  mutual  impedance  is  very  small.  Also  of  interest 


1 

’ Jn\  —  n ? 


(25) 
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Fig,  6.  (a)  Mutual  impedance  versus  separation  for  dipoles  in  broadside 

configuration  on  a  single  layer,  (b)  Mutual  impedance  versus  separalion  Tor  Fig.  8 .  Mutual  impedance  versus  separation  for  dipoles  in  endfire  configura- 
dipoles  in  endfire  configuration  on  a  single  layer  wi,h  a  superstate  used  to  eliminate  surface  waves. 
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Fig.  9.  Mutual  impedance  versus  separation  for  dipoles  in  endfire  configura¬ 
tion  embedded  in  the  middle  of  the  substrate. 

is  the  almost  monotonic  nature  of  the  curves  in  this  case, 
which  can  be  explained  by  assuming  the  superstrate  acts  as  a 
conductor  for  the  space- wave  fields,  and  using  image  theory. 
The  images  will  add  in  phase  for  z0  and  B  -  z0  an  odd  number 
of  electrical  quarter-wavelengths  in  the  substrate. 

VI.  Dipole  Excitation  by  Transmission-Line  Coupling 

Another  configuration  which  has  been  analyzed  using  the 
moment  approach  is  the  strip  dipole  excited  by  electromag¬ 
netic  coupling  to  a  microstrip  transmission  line,  shown  in  Fig. 
10.  This  problem  has  been  successfully  analyzed  with  the 
moment  method  in  the  past  for  the  line  and  dipole  within  a 
single  layer  (17J.  A  6-gap  source  is  placed  near  the  end  of  the 
line  as  shown,  and  the  same  procedure  as  for  the  mutual 
impedance  is  used  to  find  the  line  and  dipole  currents.  From 
this  information  the  normalized  admittance  Y/Y0  may  be 
found,  referred  to  the  dipole-end  of  the  line.  For  fixed  line  and 
dipole  heights  Zi  and  Zd,  an  input  match  may  be  achieved  by 
varying  the  dipole  length  L  and  offset  in  either  the  x  or  y 
directions,  provided  Zl  and  Zd  are  sufficiently  close.  An  offset 
variation  in  the  x  direction  was  chosen  here  in  order  to  take 
advantage  of  an  efficient  matrix  storage  procedure  for  the 
reaction  values.  The  offset  is  then  taken  in  integral  number  of 
subsections.  Ax  =  Nd. 

One  of  the  main  objectives  in  feeding  a  dipole  in  this  manner 
is  to  minimize  the  radiation  from  the  line  while  maintaining  a 
significant  radiation  from  the  dipole.  This  implies  keeping  Zl 
small,  while  simultaneously  keeping  a  prescribed  level  of 
coupling  between  the  two.  Also,  to  have  a  significant 
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Fig.  10.  Geometry  of  coupling  betweeh  dipole  and  transmission  line,  (a) 
Dipole  on  lop  of  supersiralc.  (b)  Dipole  av  Ihe  interface. 

bandwidth,  Zd  should  not  be  too  small.  However,  decreasing 
Zi  lowers  the  coupling  for  a  fixed  Zd ,  so  it  is  important  to 
investigate  ways  of  increasing  the  coupling.  One  approach  is 
to  stack  multiple  dipoles  above  the  line  [1 8] .  Another  approach 
is  to  choose  a  different  material  between  the  line  and  dipole. 
Qualitatively,  the  coupling  between  the  line  and  dipole  may  be 
regarded  as  largely  capacitive,  which  suggests  an  improved 
coupling  using  e2  >  «i  in  the  geometry  of  Fig.  10(a).  To  verify 
this,  the  end  admittance  is  plotted  versus  offset  for  different 
dipole  lengths  in  Figure  1 1(a)  for  Ihe  case  of  e*  =  r2  =  2.1. 
As  can  be  seen,  the  coupling  is  not  sufficient  to  allow  for  an 
input  match  with  the  layer  thicknesses  used.  In  Fig.  11(b)  the 
same  layer  electrical  thicknesses  are  used,  but  with  e2  =  10.0. 
In  this  case  the  coupling  is  sufficiently  strong  to  allow  for  an 
input  match,  for  L/\o  =  0.186  and  Ax/Xo  =  0.123.  In  Fig. 
11(c)  a  magnetic  superstrate  is  used,  with  p2  =  1 0.0.  The 
coupling  in  this  case  is  not  significantly  different  than  the  case 
of  Fig.  11(a),  supporting  the  conclusion  of  a  capacitive-like 
coupling  mechanism.  A  similar  improvement  in  coupling  is 
possible  when  using  stacked  dipoles.  The  stacked  dipole 
approach  allows  fo  an  even  greater  improvement  in  band¬ 
width,  but  the  total  physical  layer  thickness  will  be  larger 
when  using  e2  =  e,. 

The  dipole  can  also  be  placed  at  the  interface  and  excited  by 
a  line  in  the  subsiraie,  as  shown  in  Fig.  10(b).  In  this  case  the 
superstrate  is  not  improving  coupling,  but  may  be  used  to 
increase  the  efficiency  of  the  dipole  or  achieve  other  desired 
effects  [1). 
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Fig.  11.  (a)  Normalized  end  admittance  versus  offset  for  different  dipole  lengths  with  t2  =  <,  =  2.1.  The  offset  along  each  curve  is 
specified  by  Ax  =  Nd  with  N  =  A,  8,  12,  16,  20  corresponding  to  points  O,  □  ,  A,  X,  V,  respectively,  (b)  Normalized  end 
admittance  versus  offset  for  different  dipole  lengths  with  <j  =  10.0.  =  2.1.  The  dashed  curve  shows  a  dipole  length  which  allows 

for  an  input  match.  The  offset  along  each  curve  is  specified  by  Ax  =  N-d  with  N  -  3,5,  7,  10,  1 1,  12  corresponding  to  points  O, 
D,  A ,  X,  V ,  +  ,  respectively,  (c)  Normalized  end  admituncc  versus  offset  for  different  dipole  lengths  with  n2  =  (0.0,  f ,  =  2.1. 
The  offset  along  each  curve  is  specified  by  Ax  =  Nd  with  N  =  8,  12,  14,  16,  18.  20  corresponding  to  points  O,  □.  A,  X.  +  , 
respectively. 
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VII.  Conclusion 

Four  geometries  consisting  of  planar  strips  in  a  substratc- 
supcrstratc  configuration  have  been  analyzed,  including  ihc 
microstrip  transmission  line,  the  center-fed  dipole,  the  mutual 
impedance  between  two  dipoles,  and  the  transmission-line 
coupled  dipole.  The  starting  point  in  the  analysis  of  each  of 
these  was  the  plane-wave  spectrum  formulation  for  the 
reaction  between  two  current  sources.  A  real-space  evaluation 
technique  for  finding  the  electric  field  was  also  used,  which 
improved  the  efficiency  of  the  moment-method  calculation  for 
mutual  impedance,  and  also  allowed  for  an  easy  determination 
of  the  mutual  impedance  between  infinitely  narrow  dipoles. 

For  the  microstrip  transmission  line,  the  amount  of  disper¬ 
sion  was  seen  to  significantly  increase  as  a  superstrate  is  added 
when  62  is  much  larger  than  6, .  The  dispersion  decreases  as  the 
line  width  decreases,  although  must  always  be  above  the 
largest  surface-wave  pole.  In  the  case  of  the  center-fed  dipole, 
the  resonant  length  and  input  resistance  were  found  from  the 
moment  method  and  compared  with  simple  approximate 
formulas.  These  approximate  formulas  allow  for  an  easy 
extension  of  Hertzian  dipole  results  to  the  resonant  strip 
dipole.  The  mutual  impedance  between  dipoles  was  then 
investigated  using  both  the  moment  method  and  the  emf 
method.  The  emf  method  allows  for  an  efficient  computation 
for  the  case  of  very  narrow  dipoles,  and  gives  accurate  results 
when  the  proper  value  of  ke  is  used  in  the  assumed  current.  It 
was  demonstrated  that  a  superstrate  may  be  used  to  signifi¬ 
cantly  decrease  the  mutual  impedance  between  dipoles, 
especially  for  dipoles  in  an  endfire  configuration,  by  reducing 
or  eliminating  surface  waves.  Finally,  results  were  shown  for 
the  case  of  a  dipole  excited  by  coupling  to  a  transmission  line, 
and  it  was  verified  that  the  coupling  mechanism  is  predomi¬ 
nantly  capacitive.  Hence  the  coupling  is  improved  by  using  a 
material  with  a  higher  dielectric  constant  between  the  line  and 
dipole. 
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RESONANT  FREQUENCY  OF  ELECTRIcAuLV 
THICK  RECTANGULAR  MICROSTRIP 
ANTENNAS 


't*hc  function/(£r)  describes  the  influence  of  the  dielectric  inho- 
mogeneity  on  A  W.  It  is  taken  to  be  the  same  as  previously 
found  for  the  open  end  of  the  microstrip  line: 


Indexing  terms;  Antennas,  Microstrip 

An  algebraic  formulation  for  the  resonant  frequency  of  a 
rectangular  microstrip  antenna  has  been  derived  which  is 
valid  for  electrically  thick  substrates.  Predicted  values  of  res¬ 
onance  are  compared  with  previous  theories,  as  well  as 
actual  experimental  measurements,  for  radiators  approaching 
one-quarter  wavelength  in  thickness 


/(£,)  =  +  0-3 


2  t„{!)  -  0  258 


(4) 


In  cqn.  4,  e,f(/)  means  that  e.rr  is  the  effective  dielectric  constant 
of  a  microstrip  line  of  width  l.  Using' eqns.  2  and  3  in  the 
expression  for  Weff,  we  obtain 


Introduction:  Rectangular  microstrip  patch  antennas  are 
among  the  most  widely  used  printed-circuit  antennas.  A  com¬ 
prehensive  analysis  of  the  configuration  has  been  carried  out 
by  a  number  of  authors,  resulting  in  simple  design  equations 
for  the  resonant  frequency,  input  impedance,  beamwidth  and 
bandwidth.12  The  bandwidth  of  the  rectangular  patch  is  typi¬ 
cally  very  small  for  low-profile,  electrically  thin  configurations. 
One  of  the  techniques  to  increase  the  bandwidth  is  to  increase 
the  thickness  proportionately.  However,  most  of  the  design 
equations  for  the  resonant  frequency  are  only  accurate  for 
thin  patches,  normally  of  the  order  of  h/X  £  0  02,  where  h  is 
the  thickness  of  the  dielectric  substrate  and  X  is  the  wave¬ 
length  in  the  substrate.  In  this  letter  we  present  a  simple  semi- 
empirical  design  equation  for  the  resonant  frequency  of  an 
electrically  thick  rectangular  patch. 

Analysis:  The  design  equation  is  based  on  the  simple  observa¬ 
tion  that  the  open-end  extension  of  a  rectangular  patch  is 
related  to  the  edge  extension  of  a  similar  transmission  line.  In 
the  case  of  striplines  or  triplate  lines  these  extensions  are  given 
by  b  In  (2/tc),  independent  of  the  dielectric  constant  of  the 
substrate  and  frequency  for  wide  strips  (W/b  P  0T0). 
However,  in  the  case  of  microstrip  such  a  simple  behaviour  is 
not  expected  to  exist  because  of  the  inhomogeneity  of  the 
medium  above  the  ground  plane. 


h -  - Hra 

Fig.  1  Geometry  of  rectangular  microstrip  antenna 

The  resonant  frequency  fmn  of  a  rectangular  patch  of  width 
V/  and  length  /,  both  comparable  to  2/2,  is  given  by 


fmn  = 


*te)i 


(l) 


where  is  the  dynamic  dielectric  constant  for  the  patch 
defined  by  Wolff  and  Knoppik,1  and  Wtff  and  leff  are  the 
effective  dimensions  taking  into  account  the  energy  stored  in 
the  fringing  fields.  The  physical  and  effective  dimensions  of  the 
patch  are  shown  in  Fig.  I,  To  determine  Wrff  we  consider 
W 'tff  =  W  +  A W,  where  AW  is  the  open-end  effect,  which 
according  to  the  observation  above  is  related  to  the  edge  effect 
of  a  microstrip  line  of  width  /.  This  is  given  by 


A  tv  -  /) 


(2) 


”  l  2  /  -  0-2 

Similarly,  \ef{  is  given  by 
\,.ff  =  I  + 


3_ 

258 


K,  ~  ,rr(W]  +  0  3 
t:,r{W)  -  0-258 


(5) 


(6) 


Results:  For  evaluating  r.tt  the  simple  expression  due  to 
Schneider5  was  used: 


10  1~U2 
W/h_ 


O) 


For  calculating  Z0  for  a  given  value  of  W/h,  the  equations 
provided  by  Wheeler  were  used.4  The  effect  of  dispersion  was 
ignored.  The  resonant  frequency /,0  obtained  using  the  design 
equation  given  here  is  compared  with  that  due  to  James  el  al 2 


Tabic  1  CALCULATED  AND  MEASURED 


NORMALISED  RESONANT  FREQUENCIES 
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IV 

w. 

fo 

fmn 

fi» 

fhn 

5-70 

3-80 

0037 

2  586 

0893 

0918 

0-889 

0-920 

4  55 

3  05 

0-047 

3  222 

0897 

0899 

0  866 

0-900 

2-95 

1  95 

0-068 

5  039 

0  841 

0  848 

0816 

0-861 

1  95 

1  30 

0-094 

7  559 

0-773 

0-786 

0-754 

0-810 

MO 

M0 

o-no 

8  934 

0-761 

0-754 

0-724 

0  785 

1  40 

0-90 

0  125 

10-919 

0  705 

0-715 

0-683 

0-750 

1  20 

0-80 

0  141 

12-284 

0  673 

0-693 

0-662 

0  733 

1  05 

0-70 

0-148 

14-039 

0651 

0664 

0  633 

0  710 

0-90 

0-50 

0-166 

16-378 

0-626 

0  629 

0-600 

0  683 

1  70 

l  to 

0-229 

8  934 

0-529 

0  513 

0484 

0-590 

r.,  -  2  33.  h  =  0-31 75cm  (except  last  entry,  where  h  =  0-9525 cm) 


and  Hammerstad15  in  Table  I.  The  set  of  parameters  for  the 
rectangular  patch  arc  those  given  by  Chang  et  al6  In  each 
case  a  normalised  resonant  frequency  is  found  by  dividing  the 
actual  value  of  resonant  frequency  by  /0,  the  zero-order  pre¬ 
diction  assuming  the  substrate  is  infinitesimally  thin  and  that 
W  is  one-half  wavelength  in  the  dielectric.  For  units  of  c  in 

i  0 


0  9 


£0  8 

£ 

3 

CT 

i07 

•o 


where  is  the  equivalent  dimension  obtained  from  the  planar 
waveguide  model  :J 


1 20n/i 

=  zmJTJd 


(3) 


0  4  I _ 4 _ , _ , _ 1 _ , 

0  0  I  0  ? 

[53523  electrical  thickness,  wavelengths 

Fig.  2  Calculated  and  measured  normalised  resonant  frequencies 
r.r  -  2-33,  h  =  0  3175cm  (except  last  entry,  where  h  =  0-9525 cm) 
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inctrc.s/socond  and  ilio.sc  of  W  in  cculimehcs.  Jt)  — 
\5}{W Jer)GHz.  The  entries  for  JmH,  f„,  fJlt  and  fhm  thus  rep¬ 
resent,  respectively,  the  values  measured,  calculated  by  this 
method,  calculated  by  Jaimes  at  al.,  and  calculated  by  Ham- 
merstad  divided  by  the  appropriate  f0.  The  graphical  repre¬ 
sentation  of  these  data  in  Fig.  2  shows  that  the  resonant 
frequencies  predicted  here  lie  in  between  those  calculated  by 
the  other  two  methods,  and  are  closer  to  measured  values  for 
most  of  the  cases.  An  overall  accuracy  of  approximately  3% 
was  found  for  values  of  electrical  thickness  hfX  <  0-229. 

Acknowledfimenix;  This  work  was  supported  by  the  US  Army 
Research  Office  through  contract  DAAG29-H4-K-0I66  and  by 
a  grant  from  the  Texas  Advanced  Technology  Research  Pro¬ 
gramme. 

R.  GARG  Iflth  Auqust  /P«V7 

S.  A.  LONG 

Department  of  Electrical  Enginceriiuj 
U  nwersity  of  Houston 
Houston,  TX  77004,  USA 


References 


1  bahl,  ).  and  bmarTia,  »>..  'Microstrip  antennas'  (Artech  House, 
1980) 

2  James.  J.  R..  hall,  h.  s„  and  wood.  C.  ‘Microstrip  antennas — 
theory  and  design’  (Peter  Pcregrinus,  1981) 

3  wolh-,  t..  and  KNOi*P»K.  N..  ’Rectangular  and  circular  nucrosvnp 
disk  capacitors  and  resonators’.  IEEE  Trans.,  1974,  MTT-22,  pp. 
857-864 

4  GOTTA.  K.  <:..  gakg.  R-,  and  bahl,  i.  ‘Microstrip  lines  and  slot- 
lines’  (Arlcch  House,  1979) 

5  SCHNEIDER,  M.  v.  ‘Microstrip  lines  for  microwave  integrated  cir¬ 
cuits'.  Belt  Syst.  Tech.  J..  1969,  48,  pp.  1422-1444 

6  CHANG.  L\.  LONG.  s.  a.,  and  RICHARDS,  w.  F. :  'An  experimental 
investigation  of  electrically  thick  rectangular  microstrip  antennas'. 
IEEE  Trans..  1986.  AP-34,  pp.  767-772 


APPENDIX  H 


570  IEEE  TRANSACTIONS  ON  ANTENNAS  AND  PROPAGATION.  VOL  36.  NO.  4.  APRJL  1988 

Communications _ 


An  Improved  Formula  for  the  Resonant  Frequencies  of 
the  Triangular  Microstrip  Patch  Antenna 

RAMESH  GARG  AND  STUART  A.  LONG.  SENIOR  MEMBER,  IEEE 

Abstract — An  improved  formula  is  derived  to  predict  the  resonant 
frequency  of  the  triangular  microstrip  patch  antenna.  The  technique 
previously  used  for  the  case  of  a  circular  disk  is  modified  to  produce 
greater  accuracy  than  previous  investigations. 

Introduction 

Recently,  Daheie  and  Lee  [1]  have  reported  on  the  resonant 
frequencies  of  an  equilateral  triangular  microstrip  patch.  The 
calculated  values  of  these  frequencies,  based  on  an  effective  length 
aeff  and  an  effective  relative  permittivity  ceff ,  have  been  compared  with 
measured  values.  A  comparison  with  measured  values,  however, 
shows  some  reasonable  difference  between  the  two.  Here,  we  show 
that  if  a  different  effective  length  and  just  the  actual  substrate 
permittivity  arc  used,  the  computed  values  arc  much  closer  to  the 
measured  values. 

Calculations 

To  compute  the  resonant  frequency  of  an  equilateral  triangular 
patch,  the  effect  of  the  fringing  field  can  bo  assumed  the  same  as  that 
of  an  equivalent  circular  disk  which  has  the  same  metalized  area  as 
the  given  triangular  patch.  Previously,  the  effective  radius  at  of  the 
circular  disk  was  given  [2J: 


Resonant  Frequencies  of  an  Equilateral  Triangular  Patch 
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ae  =  a 


In— +  1.7726 
2  h 
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where  a  is  the  radius  of  the  disk  and  h  is  the  substrate  thickness.  This 
formula  can  then  be  modified  to  provide  a  reasonably  accurate 
prediction  of  the  resonant  frequency  of  the  triangular  patch  by 
substituting  an  equivalent  sidelength,  into  (1)  in  the  place  of  a  to 
then  produce  a  new  effective  sidelength.  This  can  be  done  by  defining 
at q  such  that  S  =  *fleq  where  S  is  the  area  of  the  original  triangular 
patch.  Then  this  value  of  af  can  be  used  in  the  formula  for  the 
resonant  frequency  (1J: 


2c 


'3  ae{tr)'/2 


(m2  +  mn  +  n2)vl 


C) 


where  c  is  the  velocity  of  light  and  the  integers  m,  n,  I  satisfy  the 
condition  m  +  n  +  l  =  0. 

This  approach  was  first  suggested  by  Suzuki  and  Chiba  [3]  for  the 
case  of  a  pentagonal  patch.  The  results  are  shown  in  Tabic  I,  along 
with  the  previous  theories  [1]  and  experimental  data.  As  can  be  seen 
the  new  method  consistently  predicts  resonant  frequencies  more 
accurately  through  the  first  five  resonances. 
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A  Leaky-Wave  Analysis  of  the  High-Gain  Printed 
Antenna  Configuration 

DAVID  R.  JACKSON,  member,  ieee,  and  ARTHUR  A.  OLINER,  life  fellow,  if.ee 


Abstract — A  leaky-wave  analysis  is  used  to  explain  the  narrow-beam 
"resonance  gain"  phenomenon,  in  which  narrow  beams  may  be  produced 
from  a  printed  antenna  element  in  a  substrale-superslrate  geometry.  The 
leaky-wave  approach  furnishes  insight  into  the  physical  processes  in- 
volved,  provides  new  practical  information,  and  permits  the  derivation  of 
simple  asymptotic  formulas  for  the  leaky-wave  properties  and  for  the 
radiation  behavior.  Results  arc  presented  as  a  function  of  frequency,  the 
scan  angle,  and  (he  e  of  the  superstate. 


I.  Introduction 

RECENTLY,  a  method  for  increasing  the  gain  of  a 
printed-circuit  antenna  involving  the  use  of  a  superstrate 
layer  was  discussed  [1],  In  this  method  the  antenna  element  is 
embedded  within  a  two-layered  geometry  consisting  of  a 
bottom  (substrate)  layer  of  thickness  b  having  relative  permit¬ 
tivity  and  permeability  e t .  mi,  and  a  top  (superstrate)  layer 
having  parameters  e2,  M2,  as  shown  in  Fig.  1.  As  explained  in 
[1],  a  narrow  beam  about  any  desired  angle  8P  may  be  obtained 
for  e2  >  1  or  M2  **  1  when  the  layer  thicknesses  are  chosen 
properly.  In  particular,  the  following  dual  types  of  "resonance 
conditions"  are  discussed  in  [1], 


Case  I: 


Xq 


hi 

Xq 


Vl  -  sin2  0  //rf  =  — 

1  2 

(2/P-  1) 


V 1  -  sin2  0p/n\ 


for  e2  >  1. 
Case  2: 


n,b 

Xo 

n2t 


'/l  -sin2  8p/n] 


—  V 1  -  sin2  8p/n\  = 
Xo 


(2m  -  1) 
4 

(2p~  I) 


U) 

(2) 


(3) 

(4) 


for  m  >  1,  where  rt|  =  Ve  t M t .  n2  =  Vr2M2  and  in,  n  are 
positive  integers. 

When  (1),  (2)  or  (3),  (4)  are  satisfied,  the  gain  at  8P  becomes 
increasingly  large  as  r2  -»  o=  or  m2  -*  »,  respectively,  with  a 
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Substrate  -  Superstrate  Geometry 


Fig.  1.  Two-layered  geometry  with  embedded  dtpolc  antenna  element  from 
which  narrow  beam  can  be  radiated.  e3  (or  ^2)  of  superstrate  (upper  layer) 
must  be  greater  than  e,  (or  n,)  of  substrate.  Horizontal  antenna  element  is 
located  at  height  Zo  and  peak  of  radiated  beam  occurs  at  angle  6P. 

resultant  decrease  in  beamwidth.  Asymptotic  formulas  for 
gain  and  beamwidth  are  presented  in  [1).  Although  providing 
for  convenient  formulas,  this  derivation  does  not  by  itself 
provide  any  insight  into  the  fundamental  cause  of  the  narrow- 
beam  phenomenon. 

In  this  article  it  is  shown  that  the  resonance  gain  effect  may 
be  elegantly  described  in  terms  of  leaky  waves.  It  is 
demonstrated  that  the  phenomenon  is  attributable  to  the 
presence  of  both  transverse  electric  (TE)-  and  transverse 
magnetic  (TM)-mode  leaky  waves  which  are  excited  on  the 
structure.  Asymptotic  formulas  for  the  leaky-wave  properties 
are  derived,  which  demonstrate  that  the  attenuation  rate  of  the 
leaky  waves  decreases  as  the  gain  of  the  pattern  increases.  The 
migration  with  frequency  of  the  leaky-wave  poles  in  the 
steepest-descent  plane  is  also  discussed,  and  it  is  shown  that 
several  important  physical  effects  are  confirmed  by  that 
investigation.  Radiation  patterns  based  on  the  relevant  leaky 
wave  are  then  compared  with  the  exact  patterns  to  demonstrate 
the  dominant  role  of  the  leaky  waves  in  determining  the 
pattern.  One  practical  result  of  this  investiga'ion  is  a  simple 
rule  for  determining  the  size  of  the  structure  necessary  to 
produce  the  narrow-beam  patterns,  a  result  which  would  be 
difficult  to  obtain  otherwise. 

Due  to  the  similarity  of  results,  only  case  1  (f  2  1)  will  be 

discussed.  This  is  the  more  practical  of  the  two  cases,  since 
low-loss  high-permeability  materials  are  difficult  to  obtain. 

II.  Leaky-Wave  Poles  and  the  Steepest-Descent 
Representation 

The  theory  of  leaky  waves  and  their  role  in  determining  the 
radiation  behavior  of  certain  antennas  is  well  established  [2]- 
[7],  To  emphasize  the  main  points,  the  analysis  for  the  case  of 
a  horizontal  Hertzian  electric  dipole  in  the  stratified  medium 
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(b) 


Fig.  2.  (a)  A,-plane,  where  k,  is  wavenumber  transverse  to  z  direction. 
Surface-wave  poles  lie  along  real  axis  (denoted  by  at’s)  and  branch  cuts 
(wiggly  lines)  emanate  from  branch  points  ±A'0.  Integration  path  T  is 
deformed  around  surface-wave  poles  in  accordance  with  conventional 
assumption  of  infinitesimal  loss  (and  e /l"  time  dependence).  Branch  cuts  are 
chosen  so  that  upper  Riemann  sheet  is  proper,  whereas  lower  Ricmann 
sheet  is  improper.  Leaky-wave  poles  occur  on  improper  (lower)  sheet,  (b) 
Steepest-descent  plane  corresponding  to  time  dependence  c*1.  Saddle-point 
is  present  at  f  =  0  on  real  axis,  and  example  is  given  of  leaky-wave  pole 
and  surface-wave  pole  captured  by  deformation  from  original  path  T  to 
stccpcst-desccnt  path. 

of  Fig.  1  will  be  considered.  Starting  with  the  classical 
Sommerfeld  solution  [8],  the  expressions  for  the  fields  are 
written  in  terms  of  the  Hertzian  magnetic  vector  potential  ft  as 

ff=ju(VxTl  (5) 

£=ft2fi  +  V(V  ■  n).  (6) 

After  solving  the  boundary-value  problem  for  fl,  the  two 
components  of  the  potential  nv  and  IT  may  be  written  in  the 
region  z  >  ft  as  (suppressing  e*'"') 

r  f(k,r) 

n,(r,  z>  T^Tre-Jk‘{:-h)H%'(krr)  dkr  (7) 

J-o>Ve{Kr) 

n;(r,  z)  =  cos  <t>  \  - ^  '  - - e~ik:t:~h)Ht2){krr)  dk, 

‘  J-„  Dr(kr)D„,(kr)  ' 

(8) 

where 

r  =  \fx2  +  y2 
kt={k2-k2)</2. 

In  these  expressions  /,  g,  Dc,  and  D,„  arc  complicated 


functions  j9]  of  the  integration  variable  kr.  Pe  and  D„t  have 
zeros  on  the  real  axis  of  the  k,  plane  in  the  range  ka  <  kr  < 
max  ( k where  k,  is  the  wavenumber  in  the  fth  layer.  These 
zeros  correspond  to  TE-  and  TM-mode  surface  waves, 
respectively.  The  integration  path  goes  around  these  poles,  as 
shown  in  Fig.  2(a).  Branch  cuts  from  the  branch  points  at  ±  k0 
are  also  necessary  to  define  the  proper  branch  interpretation  of 
k.,  which  is  Im  (ft,)  <  0,  corresponding  to  fields  that  decay  as 
z  increases. 

The  steepest-dcscent  transformation 

k,  =  k0  sin  f  (9) 

k.-kac os  f  (10) 

then  results  in  an  integration  along  the  contour  T  in  the 
steepest-descent  f-plane,  shown  in  Fig.  2(b).  This  contour 
may  be  deformed  into  the  steepest-descent  path  (SDP)  which 
passes  through  the  saddle  point  at  f  =  0,  where  R,  6,  <t> 
describe  the  observation  point  in  spherical  coordinates  with  r 
=  R  sin  0,  z  -  ft  =  R  cos  6.  The  far-field  evaluations  of  (7), 
(8)  may  be  obtained  via  the  saddle  point  method  in  the  usual 
way  [7],  It  is  important  to  recognize  that  the  single-sheeted 
steepest-descent  plane  as  a  whole  is  neither  proper  nor 
improper,  but  in  fact  portions  of  it  are  proper  and  other 
portions  improper.  The  semi-infinite  strip  for  which  0  <  (V  < 
tt/2  and  f,  <  0  is  an  improper  region  (Im  (ft£)  >  0).  Complex 
zeros  of  £>,(A',),  Dm(k,)  may  exist  in  this  region  at  f,,  +  j^„, 
which  result  in  TE-  and  TM-mode  “leaky-wave"  poles 
(LWP).  Complex  poles  can  only  occur  in  the  improper  region 
for  the  lossless  case  when  e,  >  1 ,  n,  >  1  [2] .  Hence  the  only 
other  poles  in  the  f-plane  are  the  surface  wave  poles  (SWP), 
which  lie  along  the  lines  f,  =  ±  *72,  as  shown. 

If  8  is  sufficiently  large,  some  leaky-wave  poles  will  in 
general  be  captured  during  the  path  deformation,  as  shown  in 
Fig.  2(b).  The  leaky-wave  fields  increase  exponentially  with 
z,  but  decrease  exponentially  with  R  within  the  region  in 
which  they  are  defined,  since  8  >  when  they  are  captured 
(2],  The  leaky-wave  phase  and  attenuation  constants  are  given 
by 

0  =  Re  ( krp )  =  Re  ( k0  sin  fp)  =  k0  sin  frj,  cosh  (1 1) 

a=  -Im  (k,p)=  -Im  (k0  sin  fp)=  - k0  cos  X,P  sinh  f,,,, 

(12) 

respectively,  where  (3  —  jot  is  the  pole  location  in  the  A,-plane 
(on  the  improper  sheetV  If  a  leaky-wave  pole  is  close  to  the 
real  axis  so  that  a  is  sufficiently  small,  the  leaky-wave  field 
will  dominate  the  total  field  at  the  air-dielectric  interface, 
provided  the  residue  is  not  too  small  [31-  In  this  case  the 
radiation  pattern  exhibits  a  sharply  peaked  behavior  at  6  =  8P 
where 

0  =  A'o  sin  8„.  (13) 

A  more  detailed  discussion  of  leaky-wave  properties  is  given 
in  [2],  [3). 

III.  Asymptotic  Evaluation  of  thf.  Leaky-Wavi:  Properties 

In  view  of  the  previous  discussion,  it  may  be  expected  that 
leaky  waves  could  be  responsible  for  the  high-gain  effect.  This 
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presumption  may  be  verified  by  a  straightforward,  but 
lengthy,  asymptotic  solution  for  the  leaky-wave  pole  locations 
of  the  characteristic  equations,  which  are  as  follows. 


TE-Mode: 


D,  =  0  =  cos  {kzlt)hilUkl0b)  +  (.kzlb)  cot  (A-;,i>)] 


sin  ( kgl ) 
knt 


-  (kz2t)(kz2b) 

Lw 


-jpi(ktib)(kz0t)  cot  (ktlb) 


TM-Mode: 


(14) 


Dm  =  0  =  cos  (,kz2t)[j(kz 0b)ti-(kzlb)  tan  (kz,b)] 

-  (kz2t)(kz2b) 

+y«j(*,i b)(.kz0t)  tan  (kzl 6)1  (15) 


sin  {kz2t) 
kz2l 


where 

kza=(kl~  k))ul 

kzl-{k\-ky* 

kz2=(k\-k))'n. 

The  sign  taken  for  the  square  root  of  kzl,  kz2  is  arbitrary  since 
these  signs  occur  in  even  function  fashion  in  (14)  and  (15),  but 
Im  (,kz o)  >  0  is  used  to  ensure  the  improper  pole  interpreta¬ 
tion. 

To  solve  for  complex  pole  locations  asymptotically,  k,  is 
written  as  a  perturbation  about  the  expected  high-gain  value  as 


k,  =  k0  sin  6„  + k0(b,+j&:).  (16) 

The  arguments  of  the  tan  and  cot  functions  may  then  be 
approximated  for  the  cases  8„  >  0  and  0P  =  0  and  then 
substituted  into  (14),  (15),  with  only  the  dominant  terms  kept 
as  e2  -»  oo.  After  straightforward  but  lengthy  simplification, 
the  results  summarized  below  are  obtained. 


where 


a,  =  —  —  (nj-sin!  8pyn  cot  6P  (19a) 

rmc  n  | 

b,  - - C|jj2("2_sin!  0p)'n  (sin  0.,  cos  0„)  ~ 1  (19b) 

nn r 


(19c) 


Several  important  conclusions  follow  from  these  results. 
First,  there  exists  a  TE-mode  and  a  TM-mode  leaky-wave  pole 
which  both  become  dominant  ( a  -*  0)  for  e  2  >  1 .  For  these 
dominant  poles  0/ko  —  sin  6P  when  8P  >  0,  in  agreement  with 
the  conclusion  of  a  peaked  beam  at  6P  due  to  the  leaky  waves. 
For  Sp  >  0  the  attenuation  constant  ct  behaves  as  0(l/e2).  For 
the  case  0P  =  0,  0/ko  tends  to  zero,  but  relatively  slowly, 
behaving  as  0(1  /Vel).  In  this  case,  a/k0  has  the  same 
asymptotic  form  as  0/ko,  and  hence  the  leaky-wave  pole 
asymptotically  approaches  the  origin  along  the  diagonal  line, 
in  either  the  A,-plane  or  the  f-plane.  Furthermore,  in  this  case 
the  form  is  the  same  for  both  the  TE-  and  TM-mode  poles. 

To  verify  the  validity  of  the  asymptotic  formulas,  results 
showing  a  comparison  of  the  exact  and  asymptotic  leaky-wave 
propagation  wavenumbers  for  the  TE-mode  case  are  presented 
in  Fig.  3  for  a  broadside  beam,  and  in  Fig.  4  for  a  scanned 
beam  ( 8P  =  45"),  for  the  lowest  mode  (m  =  1,  p  =  1).  The 
exact  values  are  found  by  a  secant-method  search  in  the 
complex  plane.  The  agreement  is  good  in  all  cases,  and 
improves  as  r,  increases. 

One  of  the  practical  uses  of  these  formulas  for  the  leaky- 
wave  properties  is  in  determining  the  required  structure  size  in 
the  x,  y  directions  to  realize  the  narrow  beam  patterns.  The 
size  required  is  such  that  the  leaky  waves  emanating  from  the 
dipole  are  sufficiently  attenuated  at  the  perimeter.  If  the 
structure  radius  is  rs,  then  a r,  s>  1  is  required. 


IV.  Leaky-Wave  Radiation  Patterns 


TE-Mode: 

0P  >  0: 

iS/Ao-sin  dp  (1 7a) 

a/k0~  at/e2.  (17b) 

6P  =  0: 

(S/k0~a/k0~Ci/'J72.  (17c) 

TM-Mode: 

Sp  >  0: 

0/ko~sin8p  (18a) 

a/ko-b  |/r2.  (18b) 

8„  =  0: 

g/k0~a/k0~c,/'ft2  (18c) 


The  TE-mode  leaky  wave  primarily  determines  the  Et 
component  of  the  far  field,  while  the  TM-mode  wave 
primarily  determines  Et.  The  TM-mode  leaky  wave  thus 
determines  beamwidth  in  the  £'-plane,  while  the  TE-mode 
wave  determines  beamwidth  in  the  //-plane.  Because  the 
asymptotic  form  for  a/kr,  is  the  same  for  the  TE-  and  TM- 
modes  when  Bp  -  0,  a  broadside  pencil  beam  is  produced 
which  is  circular  in  shape.  For  a  scanned  beam  the  forms  are 
different,  so  a  conical  beam  having  different  beamwidths  in 
the  two  planes  is  created. 

For  the  //-plane,  it  is  possible  to  predict  easily  the  leaky- 
wave  radiation  pattern  using  scalar  diffraction  theory.  This  is 
because  the  //-plane  pattern  of  the  dipole  is  the  same  as  that 
for  an  infinite  line  source  in  the  x  direction.  Because  of  the 
two-dimensional  nature  of  the  line-source  problem,  a  scalar 
Kirchhoff-Huygens  formulation  may  be  used  to  determine  the 
radiation  field  from  the  amplitude  distribution  of  the  leaky- 
wave  on  the  layer  surface.  Derived  in  [3],  the  formula  for  the 
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(a) 


(b) 

Fig.  3.  Comparisons  between  exact  and  asymptotic  TE-modc  leaky-wave 
propagation  wavenumbers  as  function  of  dielectric  constant  t2  of  super¬ 
state,  for  lowest  mode  with  Bp  =  0*  (broadside  beam),  (a)  Phase  constant 
(3/kn  versus  (b)  Attenuation  constant  ct/k0  versus  <5. 


(b) 

Fig.  4.  Same  as  Fig.  3,  but  for  6P  =  45’  (scanned  beam).  Note  th^ 
4(a),  for  0/A'o,  ordinate  scale  is  greatly  enlarged. 


radiated  power  pattern  is 


cos  6  2 

sin2  sin2  6 


To  provide  a  convincing  demonstration  of  the  don 
nature  of  the  leak)  waves  in  determining  the  pattern,  the 
//-plane  radiation  pattern  (from  a  reciprocity  method  ( 
compared  with  (20)  for  the  lowest  mode  (m  =  p  =  1)  ir 
5(a)  and  (b)  for  =  0'  and  $„  =  45",  respectively 
patterns  are  normalized  to  1.0  at  the  maximum), 
agreement  is  excellent  in  both  cases. 

Since  the  asymptotic  formulas  for  a/k0  do  not  depend  o 
dipole  position  Zq  but  only  on  the  layer  parameters, 
beamwidth  and  gain  of  the  pattern  are  likewise  independe 
z0  in  the  asymptotic  limit.  The  amount  of  radiated  powc 
however,  dependent  on  the  dipole  position.  The  optii 


•O/Vii'  t'KiiN  t  in j  /\i>  i  chl'iA  v.OnnoUhA 


vuv 


Fig.  5.  Comparisons  of  exact  and  leaky-wave  //-plane  radiation  patterns  for 
lowest  mode  with  e2  =  10.  (a)  Qp  =  0*.  (b)  d„  =  45*.  Agreements  are  seen 
to  be  excellent. 

dipole  height  to  maximize  radiated  power  is 

^  Vl  —  sin2  9p/n\  =  \  (2/i  —  1),  n=  1,2,  •••  (21) 

Aq  4 

where  z<>  <  b.  This  places  Zo  at  a  voltage  maximum  position  in 
a  transmission-line  model  of  the  layered  structure  (1). 

V.  Pole  Migration  With  Frequency 

Some  interesting  features  can  be  observed  as  the  leaky-wave 
pole  positions  are  traced  as  a  function  of  frequency  in  the 
steepest-descent  plane.  In  Fig.  6,  such  a  plot  is  obtained  by 


Fig.  6.  Migration  of  leaky-wave  poles  in  sleepesl-descent  plane  with 
normalized  frequency/,,  for  lowest  mode  with  =  10.  Interval  between 
data  points  is  0.10. 

introducing  a  frequency  ratio  /,,  and  then  using  layer 
thicknesses  according  to 


n,b 

——  =  0.5/, 

Ao 

(22a) 

nit 

~=  0.25/, 

Ap 

(22b) 

with  r,  =  2.1,  e2  =  10,  and  p,  =  p2  =  1. 

At  low  frequency  the  radiated  beam  has  a  maximum  at 
broadside  with  a  relatively  low  gain.  For  low  frequency  both 
the  TE-  and  TM-mode  poles  lie  significantly  below  the  real 
axis,  as  seen  in  Fig.  6.  As  f,  increases  toward  1.0  the 
broadside  beam  narrows,  with  a  maximum  broadside  gain 
obtained  for  /,  =  1 .0.  At  this  point  conditions  (22a)  and  (22b) 
then  correspond  to  the  resonance  condition  for  0P  =  0  for  the 
lowest  mode  (m  =  1,  p  =  1)  as  given  by  (1)  and  (2).  In  the 
steepest-descent  plane  both  poles  are  seen  to  come  near  to  the 
origin  at  fr  =  1.  At  this  point  -  f,  =  f,,  in  accordance  with 
the  prediction,  made  in  association  with  (1 7c)  and  (18c),  that 
for  large  e2  the  poles  would  lie  near  to  the  origin  along  the 
diagonal  line. 

As  /,  increases  past  1.0  the  beam  scans,  since  (1)  is  then 
satisfied  for  larger  values  of  Bp\  (2)  will  then  not  be  satisfied 
exactly,  but  this  feature  is  not  too  significant  in  determining 
the  resonance  condition.  It  is  primarily  the  substrate  thickness 
that  determines  the  beam  scan  angle.  The  scanned  beam  is 
conical  in  shape,  having  a  narrower  beamwidth  in  the  //-plane 
than  in  the  E-plane.  These  results  are  in  agreement  with  the 
pole  behavior  for  /,  >  1.0  in  Fig.  6.  As  /,  increases,  f, 
increases  for  both  poles.  However,  |f,|  decreases  steadily 
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for  the  TE-mode  pole  but  begins  increasing  for  the  TM-mode 
pole.  The  attenuation  constant  a  is  thus  smaller  for  the  TE- 
mode  pole,  in  accordance  with  (17b),  (18b)  which  predict  that 
the  TE-mode  pole  is  less  attenuated  when  c2  >  1. 

Also  shown  in  Fig.  6  is  a  portion  of  the  steepest  descent  path 
(dashed  line)  for  the  case  of  a  saddle  point  at  f  =  9  =  ir/2, 
with  the  path  described  by  sin  f,  cosh  f,  =  1 .  All  leaky-wave 
poles  to  the  left  of  this  line  will  be  captured  for  some  angle  0, 
while  those  to  the  right  are  never  captured.  This  path  also 
serves  as  a  dividing  line  between  fast  waves  (/J  <  kD)  to  the 
left  and  slow  waves  (/3  >  k0)  to  the  right,  seen  from  (11). 
Hence  all  leaky-waves  poles  which  are  captured,  and  thus 
directly  contribute  to  the  fields,  are  fast  waves.  From  Fig.  6  it 
is  seen  that  both  the  TE-mode  and  TM-mode  poles  cross  the 
7t/2  steepest-descent  curve,  but  at  different  frequencies.  The 
poles  continue  to  proceed  into  the  slow  wave  region,  intersect¬ 
ing  the  line  f,  =  t/2.  In  the  lossless  case,  poles  in  the  f-plane 
are  symmetrically  located  about  the  line  =  rr/2  [2].  Hence, 
at  the  intersection  frequencies  the  poles  merge  with  their 
symmetrical  counterparts,  which  are  also  improper  poles.  As 
fr  increases  further,  the  pole  pairs  split.  For  each  pair,  one 
pole  proceeds  downward  along  the  line  f,  =  x/2,  and 
represents  an  improper  surface  wave  pole,  which  is  never 
captured  for  any  angle  6.  The  other  pole  moves  upward  from 
the  intersection  point  along  the  line  f,  =  ir/2,  starting  as  an 
improper  surface-wave  pole.  This  pole  crosses  the  f  =  ir/2 
point  (real  axis)  at  the  frequency  corresponding  to  a  TE  or  TM 
surface-wave  cutoff  frequency,  for  the  respective  pole  type. 
The  pole  continues  to  move  upward  as  /,  increases,  becoming 
a  proper  surface-wave  pole.  The  two  upward-moving  poles 
may  then  be  identified  as  the  usual  TEj  and  TM2  surface-wave 
poles,  respectively  (for  the  indexing  used  here,  the  TM0  and 
TE|  modes  have  the  lowest  cutoff  frequencies  of  their 
respective  mode  types,  with  the  TM0  mode  having  zero  cutoff 
frequency). 

Although  only  the  TEi  and  TM2  poles  are  shown  in  Fig.  6, 
there  are  an  infinite  number  of  both  TE-  and  TM-mode  leaky- 
wave  poles  [6] ,  although  only  a  finite  number  are  captured  by 
the  steepest-descent  path  (and  hence  contribute  directly  to  the 
field)  at  any  given  frequency.  Only  the  TE2  and  TM2  poles  are 
dominant,  however,  with  the  others  being  further  away  from 
the  real  axis.  All  of  the  leaky  waves  exhibit  a  similar  behavior 
regarding  the  intersection  witii  a  symmetrical  image  pole  on 
the  7t/2  line,  however,  with  the  upward-moving  pole  crossing 
the  t/2  point  at  some  frequency,  corresponding  to  a  surface- 
wave  cutoff  frequency.  In  this  way  all  leaky-wave  poles  may 
be  indexed  according  to  their  corresponding  surface  wave  type 
and  number. 

VI.  Conclusion 

The  narrow-beam  resonance  gain  phenomenon  discussed  in 
(1]  is  shown  to  be  attributable  to  the  excitation  of  weakly 
attenuated  leaky  waves  on  the  structure.  Both  a  TE-  and  a  TM- 
mode  leaky  wave  become  dominant  for  f2  t>  1,  where 
subscript  2  refers  to  the  superstrate.  The  attenuation  constant  a 


of  the  leaky  waves  tends  to  zero  as  t2  increases,  beha- 
0(l/£2)  for  a  scanned  beam  ( 0P  >  0),  and  as  0(1  ft/e. 
broadside  beam  ( 6P  =  0).  The  phase  constant  /3  of  tfu 
waves  tends  to  k0  sin  8P  for  8P  >  0  and  tends  i< 
asymptotically  in  the  same  form  as  a  when  6P 
Asymptotic  formulas  for  the  explicit  behavior  have 
presented.  One  practical  application  of  these  formula 
simple  rule  for  determining  how  large  the  structure  rai 
needs  to  be  to  realize  the  narrow  beam:  ars  >  1. 

The  TE-mode  leaky  wave  determines  the  H- plane  pi 
while  the  TM-mode  leaky  wave  determines  the  E 
pattern.  In  the  //-plane,  the  pattern  due  to  the  leaky  wa 
itself  may  be  computed  from  Huygen’s  principle,  ar. 
results  show  excellent  agreement  with  the  total  rad 
pattern,  furnishing  additional  proof  that  the  fields  at  the 
interface  are  indeed  dominated  by  the  leaky  waves.  The  1 
wave  approach  discussed  in  this  paper  provides  both  si 
cant  insight  into  the  physical  processes  involved  in  the  na 
beam  resonance  gain  phenomenon,  and  methods  to  deri\ 
relevant  expressions  in  much  simpler  ways. 
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I  mm  to  25-4  mm. 
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MICROSTRIP  TRANSMISSION  LINE 
EXCITATION  OF  DIELECTRIC  RESONATOR 
ANTENNAS 


Indexing  terms :  Dielectric  resonators.  Microstrip,  Antennas, 
Dielectrics  and  dielectric  devices 


A  systematic  investigation  of  the  excitation  of  a  cylindrical 
dielectric  resonator  antenna  was  undertaken  to  illustrate  the 
possibility  of  such  a  feed  structure  and  to  characterise  the 
coupling  behaviour  and  radiation  patterns  of  the  system. 


Introduction :  Previous  work  by  McAllister1-3  has  shown  the 
dielectric  resonator  antenna  to  be  an  efficient  and  practical 
radiator,  but  this  work  concentrated  on  the  coaxial  probe  feed 
which  is  not  normally  acceptable  for  higher  frequency  applica¬ 
tions.  This  investigation  concentrates  on  the  use  of  a  micro¬ 
strip  feed  structure  which  may  be  applicable  to  frequencies 
well  within  the  millimetre  wave  bands. 

Geometry  and  structure :  The  dielectric  resonator  antenna 
simply  consists  of  a  block  of  relatively  high  permittivity 
material  (en  «  6  to  30)  situated  on  top  of  a  conducting  ground 
plane.  For  the  case  considered  in  this  work,  the  radiator  was 
chosen  to  be  cylindrical  in  shape  and  placed  on  top  of  a 
common  microstrip  substrate  for  its  excitation  as  shown  in 
Fig.  1.  The  microstrip  line  was  3-2  mm  wide  and  on  a  1-19  mm 
thick  Teflon -fibreglass  substrate  with  a  relative  dielectric  con¬ 
stant,  tti  «*  2-5,  to  give  a  characteristic  impedance  of  50 fi.  The 
cylinders  used  were  of  radius  12  7  mm,  with  heights  from 
2-54  mm  to  12-7  mm  and  relative  dielectric  constants  from  81 


Fig.  1  Geometry  of  dielectric  resonator  antenna  fed  by  microstrip  trans¬ 
mission  line 


Results:  Typical  impedance  and  insertion  loss  plots  arc  shown 
in  Fig.  2  and  Fig.  3  for  the  case  =  20  8,  =  8  5  mm, 

a  =  12-7  mm  and  x  =  5  mm.  The  resonant  frequency  is  seen  to 


Fig.  2  Impedance 

~  20-8,  d  *=  8-5 mm,  a  =  12-7 mm  and  x*5mm  (point 
1  =  2  00 GHz;  point  2  =  2-72 GHz;  point  3  =  3  00  GHz) 


Fig.  3  Insertion  loss 

£„  =  208,  d  *=  8-5  mm,  a  =  12  7  mm  and  x  **  5  mm 

be  2-72  GHz  while  the  predicted  theoretical  resonance  should 
occur  at  2  46  GHz  (a  resulting  error  of  10-6%).  Radiation  pat¬ 
terns  for  both  the  E-plane  (the  plane  going  through  the  centre 
of  the  microstrip  line  and  perpendicular  to  the  substrate 
surface)  and  H-plane  are  shown  in  Fig.  4  and  are  typical  of  tbe 
performance  of  the  antenna  as  long  as  the  ratio  d/a  does  not 
exceed  1-5. 


The  critical  parameter  that  determines  the  amount  of  coup¬ 
ling  and  the  particular  mode  excited  by  the  microstrip  trans¬ 
mission  line  is  the  overlap  distance  x.  Its  effect  on  the  coupling 
is  illustrated  in  Fig.  5  which  shows  the  insertion  loss  against 
overlap  distance  for  the  cylinder  mentioned  in  the  previous 
pa/agraph.  It  can  be  seen  that  the  strongest  coupling  occurs 


x.  mm  CSF5 

Fig.  5  /  nsertion  loss  against  overlap  distance  x 
t„  =  20  8,  d  =  8-5  mm,  and  a  =  12-7  mm 


around  x  —  5  mm  and  x  =  21  mm,  which  are  the  points  where 
the  termination  of  the  line  is  slightly  less  than  one  quarter  of  a 
dielectric  wavelength  from  the  edge  of  the  cylinder.  It  is  also 
noted  that  the  radiation  patterns  deteriorated  around 


x  =  15  mm,  where  the  end  of  the  line  was  near  the  centre  of 
the  cylinder  and  thus  not  able  to  excite  the  proper  lowcst- 
order  mode.  More  detailed  results  arc  available  in  Reference  4. 

The  investigation  has  shown  that  a  microstrip  line  can  be 
strongly  coupled  to  the  dielectric  resonator  antenna.  This 
coupling  results  in  an  efficient  and  practical  antenna  system 
for  frequencies  well  within  the  millimetre  waveband. 
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